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nearly  Identical  beyond  about  6  Hz.  Spectral  ratios  are  used  to  briefly  examine  seismic 
source  properties  and  the  partitioning  of  energy  between  Pn  and  Pg. 

Freqeuncy-wavenumber  analysis  at  the  12-element  array  is  used  to  obtain  estimates  of 
signal  gain,  phase  velocity,  and  source  azimuth.  This  analysis  reveals  frequency-dependent 
biases  in  velocity  and  azimuth  of  the  coherent  Pn  and  Pg  arrivals.  Incoherent  scattering 
occurs  after  a  little  over  one  second  into  the  Pn  wavefield  and  throughout  the  Pg  wavefield. 
The  lack  of  signal  correlation  on  the  horizontal  components  severely  limits  the  utility 
of  3-comporient  processing. 

The  principal  factor  governing  array  performance  is  signal  correlation,  and  it  is 
examined  here  in  terms  of  spatial  coherence  estimates.  The  coherence  is  found  to  vary 
between  the  three  sites.  In  all  cases  the  coherence  of  Pn  is  greater  than  that  for  Pg. 

The  coherence  estimates  are  used  to  construct  spatially-continuous ,  frequency-dependent 
models  of  cross-spectra,  which  can  be  used  to  simulate  array  processing  performance  for 
arbitrary  sensor  configurations.  Simulations  are  used  to  rank  the  three  loctions  in 
terms  of  their  potential  as  regional  monitoring  array  sites. 
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Chapter  1 


Introduction 

This  dissertation  addresses  seismic  monitoring  of  underground  nuclear  weapons  testing. 
Specifically,  the  focus  here  is  on  array  analysis  of  high-frequency  regional  seismic  data  from 
nuclear  weapons  tests  conducted  at  the  Nevada  Test  Site.  Before  introducing  the  technical 
aspects  of  this  research,  some  background  information  is  provided  on  nuclear  testing.  More 
complete  background  can  be  found  in  Bolt  (1976),  Dahlman  &  Israclson  (1977),  and  U.S. 
Congress  (1988). 

1.1  Background  on  Nuclear  Testing 

There  are  currently  five  countries  in  the  world  with  acknowledged  nuclear  weapons  capa¬ 
bility.  They  are  the  United  States,  the  Soviet  Union,  Great  Britain,  China,  and  France.  Each  of 
these  countries  maintains  a  nuclear  testing  program.  In  the  United  States  nuclear  weapons  test¬ 
ing  takes  place  almost  exclusively  at  the  Nevada  Test  Site  (NTS)  in  southwest  Nevada.  The 
majority  of  these  tests  are  conducted  for  the  purposes  of  weapons  development:  75%  to  80%  of 
the  nuclear  tests  conducted  each  year  contribute  to  the  engineering  of  specific  new  warheads  or 
new  weapons  systems.  The  remaining  tests  are  conducted  to  determine  llie  survivability  of 
military  systems  under  the  effects  of  nuclear  explosions,  to  improve  the  understanding  of  the 
physical  phenomena  associated  with  nuclear  explosions,  and  to  ensure  that  existing  weapons 
systems  are  working  correctly.  It  has  been  argued  that  continued  nuclear  testing  is  required  for 
national  security  and  to  enhance  the  effcctivcnass  of  nuclear  deterrence.  However,  these  argu¬ 
ments  have  been  in  question  for  more  than  thirty  years,  since  the  beginning  of  negotiations 
toward  a  Comprehensive  Test  Ban  Treaty  (CTBT),  banning  all  nuclear  testing.  There  was  at 
that  time,  and  still  is,  a  great  concern  that  the  continued  development  of  more  powerful  and 
sophisticated  nuclear  weapons  svstems  through  continued  testing  can  increase  both  the 


1 


likelihood  and  the  destructive  power  of  a  nuclear  war.  Today,  despite  extensive  negotiations 
and  a  great  deal  of  progress,  a  Comprehensive  Test  Ban  Treaty  has  not  yet  been  achieved. 

Efforts  were  made  by  the  United  States,  Britain,  and  the  Soviet  Union  between  1958  and 
1963  to  achieve  a  CTBT,  however,  ostensibly  because  of  disagreements  on  verification  pro¬ 
cedures,  no  such  treaty  was  achieved.  These  parties  instead  produced  the  1963  Limited  Test 
Ban  Treaty  G-TBT),  which  banned  testing  in  the  atmosphere,  in  outer  space,  and  under  water, 
but  placed  no  limit  on  the  size,  or  yield,  of  nuclear  tests.  This  treaty  greatly  reduced  the 
danger  of  radioactive  fallout  and  has  since  been  signed  by  over  HX)  countries,  though  China 
and  Fiance  arc  not  currently  among  them.  The  Nuclear  Non-Proliferation  Treaty  (NPT)  of 
1970  has  served  as  a  deterrent  to  conducting  nuclear  weapons  testing  of  any  kind.  Non¬ 
nuclear  weapons  countnes  who  are  party  to  this  treaty  (there  are  currently  138)  pledge  not  to 
use  nuclear  energy  for  the  purpose  of  weapons  development.  The  treaty  further  states  that  the 
three  weapon  states  that  arc  party  to  the  treaty,  the  United  States,  Britain,  and  the  Soviet 
Union,  are  to  work  toward  nuclear  disarmament  and  a  discontinuance  of  nuclear  testing.  The 
two  other  declared  nuclear-weapon  states  of  China  and  France  have  not  signed  tJiis  treaty.  In 
1995  the  NFT  will  be  examined  to  decide  to  what  extent  to  reform,  strengthen,  or  expand  it. 
Yield  limits  were  finally  placed  on  nuclear  tests  by  the  1974  Threshold  Test  Ban  'freaty 
(TTBT),  which  bans  underground  tests  by  the  United  States  and  Soviet  Union  having  an  explo¬ 
sive  yield  greater  than  150  kilotons  (kt).  The  1976  Peaceful  Nuclear  Explosions  Treaty 
(PNET)  restricts  individual  peaceful  nuclear  explosions  by  the  United  States  and  the  Soviet 
Union  (for  engineenng  and  demolition  purpo.ses)  to  yields  also  no  greater  that  150  ki.  (One  kt 
represents  l(XX)  tons  of  TNT,  or  about  a  million  .sticks  of  dynamite.  The  yield  of  the  tximbs 
dropped  on  Japan  in  World  War  11  was  about  13  kt.)  To  date,  the  Threshold  Test  Ban  and 
Peaceful  Nuclear  Explosion  Treaties  have  not  been  ratified,  though  both  parlies  have  apparently 
abided  by  the  set  regulations. 

In  the  late  1970’s,  during  the  Carter  admini.stration,  a  good  deal  of  progress  was  made  in 
the  tri-lateral  negotiations  for  a  comprehensive  test  ban  between  the  United  Slates,  Britain,  and 
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the  Soviet  Union.  However  a  setback  to  the  furthering  of  these  negotiations,  and  to  U.S.- 
Soviet  relations  in  general,  came  in  July  1982,  when  the  Reagan  Administration  announced 
that  the  United  States  was  withdrawing  from  efforts  to  ban  all  nuclear  tests.  The  principal  rea¬ 
son  given  was  that  further  development  of  U.S.  nuclear  weapons  systems  was  required  to 
ensure  the  stability  of  nuclear  deterrence.  An  additional  reason  was  the  feeling  that  the  Soviets 
had  cheated  on  the  1974  TTBT  by  exceeding  150  kt  with  a  number  of  explosions.  This  has 
been  countered  by  many  seismologists,  who  find  that  this  claim  was  based  on  a  misinterpreta¬ 
tion  of  seismic  data  (e.g.,  Sykes  and  Davis,  1987).  Also  significant  was  the  Administration’s 
lack  of  faith  in  the  verifiability  of  a  CTBT.  An  important  point  to  make  here  is  tliat,  practi¬ 
cally  spettking,  technical  verification  of  a  CTBT  cannot  be  achieved;  verification  measurements 
from  very,  very  small  explosions,  will  either  be  within  the  noise  of  the  measurements,  or  the 
uncertainty  in  identification  and  yield  estimation  will  be  so  large  as  to  make  the  measurements 
meaningless.  Instead,  one  should  think  of  a  CTBT  as  a  low-yield  TTBT  with  the  maximum 
allowable  yield  set  below  that  required  for  militarily  significant  weapons  tests,  but  not  so  low 
as  to  present  significant  verification  uncertainties.  Notice,  then,  that  the  two  primary'  factors 
controlling  the  threshold  level,  the  yields  of  militarily  significant  tests  and  the  allowable  degree 
of  uncertainty,  are  ultimately  based  upon  political,  not  technical,  judgements.  It  is  generally 
accepted  among  scismologi.sts  that,  using  seksmic  methods,  compliance  with  a  low-yield  ITBT 
can  be  successfully  monitored,  with  the  maximum  allowable  yield  set  between  about  10  kt  and 
15  kt,  and  quite  possibly  lower  depending  on  the  di.stance  of  the  monitoring  network  to  the 
te.sting  area  (U.S.  Congress,  1988). 

Despite  the  absence  of  formal  negotiations,  beginning  in  August  of  1985  and  extending 
for  19  months,  the  Soviet  Union  ob.scrvcd  a  unilateral  moratorium  on  nuclear  testing.  The 
United  States,  however,  declined  the  offer  to  join  the  moratorium.  In  May  of  1986  another 
significant  event  occurred  when  the  Natural  Re.sources  Dcfen.se  Council,  an  environmental 
organization,  and  the  Soviet  Academy  of  Sciences  signed  an  agreement  to  establish  jointly- 
manned  seismic  monitoring  stations  near  the  the  principal  testing  areas  in  each  country  (sec 


Berger  et  al. ,  1987).  The  success  of  this  venture  in  the  very  sensitive  area  of  arms  control 
came  as  a  welcome  surprise  to  many,  and  may  well  set  a  precedent  for  future  private  actions 
on  publicly,  though  not  necessarily  politically,  supported  international  issues.  Finally  in  Sep¬ 
tember  of  1987,  with  much  improvement  in  U.S.-Soviet  relations  and  under  great  political  and 
public  pressure,  a  joint  U.S. -Soviet  statement  was  issued  on  the  resumption  of  negotiations  on 
the  limitations  of  nuclear  testing.  The  first  goal  was  to  ratify  the  TTBT  and  PNET,  which 
required  agreement  on  effective  verification  measures.  Towards  this  end,  unprecedented  exper¬ 
iments  have  taken  place  in  which  U.S.  and  Soviet  scientists  have  visited  each  others’  testing 
sites  and  cooperated  in  a  collection  and  exchange  of  tc.st-moniioring  data  (e.g.,  Priestley  et  al. , 
1990).  It  now  appears  that  ratification  of  the  Threshold  Test  Ban  and  Peaceful  Nuclear  Explo¬ 
sion  Treaties  is  imminent. 

The  future  of  a  Comprehensive  Tc.st  Ban  is  much  less  certain.  During  ilic  1990  review 
of  the  NPT  a  large  amount  of  pressure  was  applied  to  the  nuclear  weapon  states  by  many  non- 
weapon  states  to  bring  a  halt  to  nuclear  testing  (Epstein,  1991).  Though  the  Soviet  Union 
claims  it  is  willing  to  do  so,  the  United  States  and  Britain  currently  remain  committed  to  test¬ 
ing  programs  and  regard  a  CTBT  as  a  long-term  objective.  Neither  Chiria  or  France  is  party  to 
the  NPT  and  their  role  in  a  future  CTBT  is  unknown.  An  additional  concern  are  the  six  other 
countries  with  advanced  nuclear  programs:  India,  Argentina,  Brazil,  South  Africa,  Pakistan, 
and  Israel.  Though  all  currently  deny  any  nuclear-weapon  intentions,  none  are  part  of  the  PNT. 
These  countries  arc  however  bound  to  the  1963  PTBT,  and  if  current  efforts  by  some  non- 
aligned  nations  towards  amending  the  1963  PTBT  into  a  CBTB  arc  eventually  successful,  lhe.se 
nuclear  nations  may  find  themselves  pany  to  a  CTBT  by  default. 

In  the  mean  time,  seismologists  arc  continuing  research  to  improve,  and  so  further  estab- 
li.sh,  the  ability  of  seismic  networks  to  monitor  compliance  with  a  low'-yield  testing  treaty.  This 
is  still  a  fairly  young  area  of  rc.scarch,  and  much  remains  to  be  done.  Some  general  background 
information  on  low-yield  seismic  monitoring  is  given  below.  A  more  complete  overview  can 
be  found  in  U.S.  Congress  (1988). 
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1.2  Regional  Seismic  Monitoring 


When  one  considers  that  a  buried  nuclear  explosion  will  generate  waves  in  the  earth 
much  in  the  same  way  as  an  earthquake,  it  is  not  surprising  that  the  science  of  seismology  has 
provided  the  principal  means  for  monitoring  underground  nuclear  explosions.  A  seismic  moni¬ 
toring  system  has  three  primary  tasks:  (1)  to  detect  that  a  seismic  event  has  taken  place,  (2)  to 
discriminate,  i.e.,  to  identify  that  seismic  event  as  either  a  nuclear  explosion  or  something  else, 
such  as  an  earthquake  or  mining  explosion,  and  (3)  if  it  is  a  nuclear  test,  to  determine  whether 
the  size  of  the  explosion  exceeds  that  allowed  by  the  prevailing  nuclear  test  ban  treaty.  A 
great  advantage  of  seismic  measurements  is  that  they  can  be  made  well  outside  of  the  immedi¬ 
ate  te.sting  area.  Other  so-callcd  “on-site”  monitoring  techniques  exist,  such  as  radiochemical 
and  hydrodynamical  methods,  but  iicsc  methods  are  much  more  intrusive  and  so  arc 
inherently  more  difficult  to  negotiate  into  a  treaty.  To  date,  U.S.  and  Soviet  compliance  with 
the  1974  TTBT  has  been  .successfully  monitored  through  the  u.sc  of  tclescismic  measurements 
made  around  the  world.  To  reliably  monitor  a  treaty  in  which  the  maximum  allowable  explo¬ 
sive  size  is  in  the  range  of  say  1-10  kt,  much  less  than  the  150  kt  currently  allowed  by  the 
TTBT,  .seismic  recordings  must  be  made  much  clo.scr  to  the  testing  sites,  at  so  called  regional 
distances,  i.e.,  at  distances  less  than  about  2(XX)  kilometers,  and  new  seismic  d  i. sc  ri  mi  nation 
methods  must  be  employed. 

In  addition  to  placing  seismic  stations  at  regional  distances  to  monitor  low-yield  explo¬ 
sions,  it  is  advantageous  to  measure  the  high  frequencies  of  ground  motion.  This  is  becau.se  (1) 
for  .small-magnitude  events  the  signal  to  noi.se  ratio  typically  increa.ses  with  increasing  fre¬ 
quency,  before  eventually  dropping  back  to  very  low  values,  (2)  the  most  likely  way  to  hide  a 
weapons  test  is  by  placing  the  explosion  in  a  large  cavity,  which  mufUes  the  explosion,  reduc¬ 
ing  the  amount  of  ground  motion  produced  at  seismic  monitoring  stations:  however,  theory 
predicts  (e.g.,  Evemden  et  ai,  1986)  and  observations  have  conlirmed  (Garbin,  1986;  Glen 
et  ai,  1987)  that  the  amount  of  muflling  is  reduced  greatly  for  high  frequencies  of  ground 
motion,  and  (3)  potentially  u.seful  methods  may  exist  to  discriminate  weapons  tests  from  non- 
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weapons  rclaied  seismic  sources,  such  as  earthquakes  and  mining  blasts,  based  on  difierences 
in  relative  Irequency  content  at  high  Irequencies  of  ground  motion  (see  rel'crenccs  below). 

The  most  imponant  regional  phases  are  referred  to  as  the  Pn,  Pg,  Sn,  and  Lg  waves. 
Pigure  1  IS  an  example  of  high-frequency  wavefomi  data  Irom  NTS  being  used  in  this  stutly. 


Pn  Pg  Lg 


0.0  20  0  40.0  60.0  80  0  100  0 

Seconds 


FIGURE  1:  Example  of  a  regional  high-frequency  seismic  recording. 

Regional  wavelields  are  influenced  very  much  by  the  regional  characteristics  of  the  crustal 
structure  through  which  they  propagate.  This  is  why,  lor  example,  Sn  is  not  a  prominent  phase 
in  the  tectonically  active  western  U.S.,  but  it  is  in  the  more  ecologically  stable  eastern  U  S.  and 
other  shicld-likc  environments.  Notice  that  regional  phases  do  not  correspond  to  discrete 
arrivals.  They  consist  instead  of  a  train  of  arrivals  lasting  from  a  few  to  tens  of  seconds. 
Beyond  a  cross-over  distance  of  about  l()()-2()()  km.  depending  on  the  cmstal  thickness,  Pn  is 
the  first  arriving  waveirain.  Its  pha.se  velocity  of  between  atxiul  7.h  and  8.4  km/scc,  depend¬ 
ing  on  the  regional  are:i.  indicates  that  Pn  derives  from  critically  refracted  enercy  traveling 
along  or  just  below  the  Mohorovicic  disconlmuily.  Multi;de  relleciions  due  to  cmstal  layering 
will  continuously  feed  energy  into  Pn  refractions  and  so  prolong  the  duration  of  the  wavetrain. 
Additionally  there  arc  indications  that  much  of  the  Pn  coda,  and  m  fact  the  coda  of  the  other 


regional  phases  as  well,  is  derived  from  scaliering  of  the  wavefields  (e.u.,  Baumgardl,  1990, 
Dainty  and  Toksoz,  1990,  this  study).  Pg  is  a  crustal  wave,  arriving  after  Pn  beyond  the  cross¬ 
over  distat.ce.  Its  apparent  velocity  is  typical  of  crustal  P  velocities.  P-SV  rellcctions  due  to 
crustal  layering  contribute  to  the  prolonged  and  complex  structure  of  Pg.  Lg  is  regarded  as 
being  made  up  of  multiply-rellcctcd  post-critical  S  waves,  and  has  phase  velocities  typical  of 
crustal  S  velocities.  Lg  is  usually  the  largest  amplitude  arrival  of  the  regional  phases  and  it 
shows  great  promise  as  a  reliable  c.stimator  of  explosive  yield  (Nuttli  1986,  Hansen  et  al. , 
1990).  All  of  the  regional  pha.scs  show  promi.se  as  effective  discriminants,  though  again 
because  of  their  .sensitivity  to  variations  in  crustal  .structure,  the  effectiveness  of  many  discrim¬ 
inants  is  region,  or  even  site,  spccilic.  The  subject  of  regional  discriminants  will  not  be  pur¬ 
sued  in  this  study,  but  further  information  can  be  found  in  Nuttli  (1981),  Pomeroy  ct  al., 
(1982),  Evemden  et  al. ,  (1986),  Bonnet  and  Murphy  (1986),  Pulli  and  Dysan  (1987),  Taylor 
et  al. ,  (1988),  Taylor  et  al. ,  (1989),  Bonnet  et  al. ,  (1989),  and  Baumgardl  and  Young  (1990). 

1.3  Regional  Arrays 

Below  is  some  background  on  sei.smic  arrays  and  a  brief  discussion  of  the  motivation 
behind  this  study.  More  detailed  discussions  will  be  found  in  the  ensuing  chapters.  The  fun¬ 
damental  point  here  is  that  low-yicid  monitoring  capability  can  be  significantly  enhanced 
through  the  use  of  arrays  of  instruments  at  regional  monitoring  sites.  Seismic  arrays  have  a 
number  of  advantages  over  single  instruments  for  the  purposes  of  monitoring  low-magnitude 
seismic  events.  By  applying  various  signal  processing  methods  to  array  data,  improved  detec¬ 
tion  capability  of  a  monitoring  site  can  be  achieved.  Specifically,  array  processing  methods  can 
be  u.sed  to  increase  the  amplitude  of  seismic  signals  relative  to  the  background  noise.  There¬ 
fore,  an  array  can  potentially  detect  much  smaller  events  than  can  a  single-instrument  site.  In 
addition,  there  a  number  of  ways  in  which  an  array  can  improve  discrimination  and  yield- 
estimation  capability.  For  example,  as  we  can  seen  in  Figure  1.1,  a  regional  wavefonn  con¬ 
tains  different  types  of  seismic  waves,  traveling  along  different  paths  through  the  earth  to 
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arrive  at  the  recording  site.  Seismic  discrimination  and  yield  estimates  require  that  the  regional 
wave  type  be  known;  this  becomes  more  difficult  as  the  signal  amplitude  decreases  relative  to 
the  background  noise  level,  as  it  would  for  very  low-magnitude  cvenus.  However,  the  propaga¬ 
tion  velocity  of  the  different  wave  types  across  the  earth’s  surface  differ,  and  with  an  array  this 
velocity  can  be  measured.  Therefore,  the  wave  type  can  be  identified  by  its  velocity  and  the 
appropriate  spectral  discriminate  and  yield  measure  can  be  applied.  Secondly,  array  processing 
methods  can  also  be  used  to  estimate  the  epicentral  location  of  seismic  sources,  much  more 
reliably  than  a  single  three-component  monitoring  site.  A  reliable  estimate  of  location  alone  can 
be  a  very  useful  discriminate  if  the  geographic  areas  of  weapons  testing  and  earthquake  and 
mining  activity  are  known.  This  location  ability  is  especially  important  if  the  event  is  so  small 
that  it  is  recorded  at  only  one  or  two  monitoring  sites,  thus  making  network  location  impossi¬ 
ble.  Additionally,  estimation  of  seismic  source  type  through  the  use  of  discriminants,  and  esti¬ 
mation  of  source  size  will  be  statistically  more  reliable  given  the  greater  number  of  recordings 
and  lowered  noise  levels  that  can  be  achieved  through  array  processing.  Arrays  also  provide  a 
more  detailed  sampling  of  the  wavcficld  and  so  allows  us  to  better  understand  the  nature  and 
extent  of  propagation  effects  which  tend  to  obscure  the  source-related  information  by  distorting 
and  randomizing  the  waveficld. 

These  kinds  of  advantages  in  monitoring  became  apparent  over  twenty  years  ago  with 
large  aperture  arrays  such  as  LASA  in  Montana  (  525  short-period  seismometers  over  an  aper¬ 
ture  of  200  km)  and  NORSAR  in  Norway  (132  short-period  seismometers  over  an  aperture  of 
100  km).  The  large-aperture  for  these  arrays  is  appropriate  for  telcseismic  monitoring,  how¬ 
ever  much  smaller  apertures  on  the  order  of  a  few  kilometers  are  required  to  sustain  the  signal 
correlation  required  for  successful  regional  high-frequency  array  processing.  Serious  work  into 
the  development  of  regional  arrays  began  in  1979  with  the  NORESS  array,  located  in  the  Bal¬ 
tic  Shield  in  .southern  Norway.  This  array  supports  a  25-scn.sor  geometry  over  an  aperture  of  3 
km.  The  Noress  array  was  later  complemented  by  the  near-identical  ARCESS  array  in  northern 
Norway  in  1987,  the  15-element  2-km  apenure  FINES  A  array  in  .southern  Finland  in  1985,  and 
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the  25-elcment  4-km  aperture  GRESS  array  in  the  Bavarian  Forest  area  of  Germany  in  1990. 

To  date,  the  bulk  of  the  work  done  in  high-frequency  seismic  array  processing  has  been 
done  on  data  recorded  at  the  Scandanavian  arrays.  These  arrays  arc  located  on  a  geologic 
shield,  which  is  thought  to  represent  the  best  propagation  environment  for  high-frequency 
waves.  In  the  tectonically  active  western  United  States  the  crustal  structure  is  profoundly 
different,  and  one  should  therefore  expect  significantly  different  high-frequency  signal  charac¬ 
teristics  for  a  regional  array  monitoring  underground  weapons  tests  from  the  Nevada  Test  Site. 
However,  little  work  has  been  done  to  study  the  performance  characteristics  of  high-frequency 
regional  arrays  within  areas  like  the  Basin  and  Range.  Given  this  lack  of  information  this 
study  was  undertaken,  wherein  we  will  examine  the  array  performance  characteristics  of  the  Pn 
and  Pg  wavcficlds  at  three  distinct  array  sites  located  at  regional  distances  from  the  Nevada 
Test  Site.  Research  of  this  kind  is  required  for  optimally  designing  future  high-frequency 
arrays  in  similar  geologic  settings,  and  for  studies  which  attempt  to  estimate  the  overall  moni¬ 
toring  capability  of  a  hypothetical  network  of  seismic  arrays  located  in  the  western  United 
States. 

In  this  study  we  will  find  that  working  against  effective  array  processing  are  compli¬ 
cations  and  irregularities  in  the  geologic  structure  which  can  scatter  and  distort  the  wave  field, 
causing  less  than  optimal  signal  enhancement,  atypical  propagation  velocities,  mislocations  of 
the  source  epicenter,  and  therefore  degraded  detection  and  discrimination  capability.  Thc.se  geo¬ 
logic  irrcgulaniics  can  exist  on  a  variety  ol  scales,  from  the  very  local  structure  surrounding 
the  array  site,  such  as  hills  and  valleys,  to  more  deep-seated  structural  variations  extending  tens 
of  kilometers  into  the  earth’s  cru.st.  Much  of  this  dissertation  will  be  concerned  with  the  extent 
and  manner  in  which  the  wavcficlds  arc  distorted  and  randomized  at  the  array  sites.  We  will 
begin  in  Chapter  2  with  a  description  and  brief  discussion  of  the  regional  data  .sets  used  in  this 
study.  In  Chapter  3  will  be  an  analysis  of  how  spectral  amplitude  estimates  vary  over  tfn_ 
aperture  of  the  arrays  and  between  the  different  array  locations.  An  appreciation  and  under¬ 
standing  of  the  sensitivity  of  spectral  amplitude  estimates  to  changes  in  receiver  location  is 
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important  given  their  use  as  regional  source  discriminants  and  yield  measures.  Also  in  Chapter 
3  we  will  examine  briefly  the  division  of  energy  between  Pn  and  Pg  and  indirectly  compare 
two  explosions  to  a  simple  explosive  source  model.  In  Chapter  4  we  will  focus  in  on  one 
array  having  a  design  similar  to  that  of  the  NORESS-type  arrays,  but  again,  located  in  a  very 
different  environment.  We  will  examine  the  Pn  and  Pg  signal  correlation  and  propagation 
characteristics  at  this  array  and  suggest  relationships  to  local  site  effects  and  scattering.  A 
more  detailed  signal  correlation  analysis  follows  in  chapter  5,  where  we  wilt  study  the  spatial 
and  frequency  dependence  of  intersensor  wavefield  coherence  at  each  of  the  three  array  sites. 
The  coherence  is  parameterized  and  used  to  simulate  array  processing  characteristics  for  arbi¬ 
trary  array  configurations.  Chapter  6  offers  recommendations  based  on  the  findings  of  this 
study.  Four  appendices  (A  through  D)  are  also  included.  Appendices  A  and  C  describe  and 
explore  the  advantages  of  the  spectral  estimation  method  used  for  amplitude  spectra  and  coher¬ 
ence.  Appendix  B  displays  an  example  of  the  relative  frequency  content  of  regional  phases 
through  bandpass  filtering.  And  finally.  Appendix  D  displays  the  coherence  data  analyzed  in 
Chapter  5. 
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Chapter  2 


Regional  Data  Sets 


2.1  Introduction 

In  this  chapter  we  will  first  describe  the  recording  sites  and  recording  system  parameters 
used  in  this  study.  Next  the  Pn  and  Pg  regional  wavefield  data  will  be  displayed  with  brief 
comments  on  waveform  differences  due  to  differences  in  recording  location,  source  parameters, 
and  wavetype.  More  detailed  analyses  of  spectral  characteristics  and  coherence  structure  will 
follow  in  subsequent  chapters.  All  the  wavefields  displayed  were  recorded  at  nearly  the  same 
epicentral  distance,  between  340  km  and  387  km. 

2.2  Recording  Sites 

Three  temporary  arrays  sited  in  the  Basin  and  Range  province  of  the  western  U.S.  are 
included  in  this  study.  Each  array  recorded  explosions  from  the  Nevada  Test  Site.  The 
configurations  of  the  array  sites  and  the  location  of  the  arrays  relative  to  the  recorded  explo¬ 
sions  are  shown  in  Figure  2.1.  We  will  also  look  briefly  at  data  from  one  explosion  recorded 
at  a  permanently-installed  regional  instrument  site  operated  by  Lawrence  Livermore  National 
Laboratory  (LLNL).  Receiver  and  source  specifications  are  given  in  Tables  1  and  2,  respec¬ 
tively. 

The  Savaliia  Mountain  array  and  the  Rice  Valley  array  were  deployed  in  the  Mojave 
Desert  of  southern  California  as  part  of  a  CALCRUST  (California  Consortium  for  Crustal  Stu¬ 
dies)  crustal  exploration  survey.  The  Savahia  Mountain  array  consisted  of  145  vertical- 
component  stations.  The  intersensor  spacing  at  this  array  site  was  25  meters.  The  maximum 
difference  in  station  elevation  was  69  meters.  The  Rice  Valley  array  consisted  of  96  vertical- 
component  stations.  The  intersensor  .spacing  at  this  array  site  was  100  meters.  The  maximum 
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difference  in  station  elevation  here  was  59  meters. 


Descriptions  of  the  geology  and  tectonic  evolution  of  the  area  surrounding  the  Mojave 
Desert  arrays  can  be  found  in  Davis  et  al.  (1982,  1988),  Howard  and  John  (1987),  a  number 
of  studies  contained  in  Coward  et  al.  (1987),  and  Wang  et  al.  (1989).  A  geologic  map  of 
this  area  is  shown  in  Figure  2.2  (after  Howard  and  John,  1987).  The  area  shown  is  dominated 
by  extensional  fault  systems  which  resulted  in  the  fonnation  of  many  of  the  mountain  ranges. 
The  Savahia  Mountain  array  was  sited  near  Savahia  Mountain,  between  the  Turtle  Mountains 
to  the  west  and  the  Whipple  Mountains  to  the  east.  The  Rice  Valley  array  was  sited  within 
Rice  Valley,  which  lies  at  the  southern  base  of  the  Turtle  Mountains.  The  difference  in  eleva¬ 
tion  between  each  array  site  and  the  surrounding  mountains  is  nominally  between  about  400 
and  1000  meters. 

Advance  notice  of  detonation  times  allowed  these  two  arrays  to  record  the  two  NTS 
explosions  SALUT  (m^  =  5.5)  and  VILLE  (m^  =  4.4).  SALUT  was  detonated  at  a  depth  of 
608  meters  in  the  Pahute  Mesa  area  of  NTS,  371  km  and  385  km  to  the  north  of  the  Savahia 
Mountain  and  Rice  Valley  arrays,  respectively.  VILLE  was  detonated  at  a  depth  of  293  meters 
in  the  Yucca  Valley  area  of  NTS,  340  km  and  355  km  to  the  north  of  the  Savahia  Mountain 
and  Rice  Valley  arrays,  respectively.  These  two  explosions  were  separated  by  a  distance  of  40 
km.  The  system  response  at  these  two  arrays  was  flat  to  velocity  from  approximately  9  Hz  up 
to  the  12-pole  anti-alias  filter  at  62.5  Hz.  The  normalized  amplitude  response  to  velocity  at 
these,  and  the  other,  instrument  sites  is  shown  in  Figure  2.3.  Not  shown  in  Figure  2.3  for 
these  two  arrays  is  a  notch  filter  applied  at  60  Hz.  The  data  were  sampled  at  250  samples  per 
second.  The  two  arrays  were  separated  by  a  distance  of  about  35  km.  The  large  number  of 
instruments  and  peculiar  configurations  make  the  Savahia  Mountain  and  Rice  Valley  arrays 
unlikely  models  for  future  low-yield  TTBT  recording  system  configurations.  However,  the 
dense  spatial  sampling  of  these  two  arrays  will  allow  us  to  study  regional  wavefield  charac¬ 
teristics  in  detail. 
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The  Ruby  Valley  array  was  a  temporary  12-station,  three-component  array  with  an  aper¬ 
ture  of  1.5  km.  The  array  was  designed  in  the  same  manner  as  the  NORESS  array,  that  is,  sta¬ 
tions  placed  in  concentric  rings  at  log-periodic  intervals  in  radius  (Kvaema,  1989).  However 
the  aperture  of  this  array  is  half  that  of  NORESS.  This  array  was  located  in  the  Basin  and 
Range  province  in  northern  Nevada  within  Ruby  Valley,  approximately  54  km  southeast  of  the 
city  of  Elko,  Nevada.  Eiescriptions  and  analyses  of  the  geology  surrounding  this  area  are  given 
by  Snoke  (198U)  and  Howard  (1980).  A  regional  geologic  map  of  the  area  is  shown  in  Figure 
2.4  (after  Snoke,  1980).  Immediately  to  the  west  of  the  array  site  lie  the  northern  Ruby  Moun¬ 
tains,  which  represent  one  of  a  number  of  metamorphic  complexes  mapped  in  this  region. 
These  mountains  rise  approximately  14(X)  meters  in  elevation  above  the  array  site.  There  was 
no  significant  change  in  station  elevation  across  the  array.  The  array  recorded  the  explosion 
HARDIN  (mb  =  5.5),  detonated  at  a  depth  of  625  m  in  the  Pahute  Mesa  area  of  the  Nevada 
Test  Site,  387  km  to  the  south  of  the  array.  The  system  response  at  this  array  site  was  flat  to 
velocity  from  approximately  5  Hz  up  to  the  5-pole  anti-alias  filter  at  50  Hz  (see  Figure  2.3). 
The  data  were  sampled  at  2{X)  samples  per  second.  The  Ruby  Valley  array,  with  its  small 
number  of  three-component  high-frequency  instruments  spread  over  1.5  km,  may  be  representa¬ 
tive  of  future  recording  system  configurations  emplaced  to  monitor  nuclear  test  sites  at  regional 
distances. 

The  explosion  HARDIN  was  also  recorded  by  the  LLNL  station  ELKO  0at-=  40.745, 
long.=  -115.239)  located  approximately  15  km  to  the  northwest  of  the  Ruby  Valley  array  at  a 
hard  rock  site  at  the  base  of  the  Ruby  Mountains.  The  system  response  at  this  site  is  flat  to 
velocity  between  approximately  0.05  Hz  and  10  Hz  and  was  sampled  at  a  rate  of  40  samples 
per  second.  Referring  to  Figure  2.3,  the  LLNL  response  is  essentially  a  low  pass  filter  of  velo¬ 
city  ground  motion  relative  to  the  array  site,  which  is  effectively  a  high  pass  filter.  The  siting 
of  the  Ruby  Valley  array  near  this  LLNL  station  was  intentional.  The  seismic  path  from  NTS 
to  Ruby  Valley  is  a  well  recorded  one  and  the  ELKO  station  is  the  quietest  of  the  LLNL 
regional  network  (Rodgers  and  Rohrer,  1987). 
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2.3  Explosions  SALUT  and  VILLE 

We’ll  begin  by  displaying  the  SALUT  and  VBLLE  Pn  and  Pg  wavefields  recorded  at  the 
Savahia  Mountain  and  Rice  Valley  arrays.  TTiough  we  will  be  analyzing  only  the  first  few 
seconds  of  these  wavefields  in  the  later  chapters,  we  will  nevertheless  display  them  here  in 
their  entirety.  Keep  in  mind  that  because  the  Savahia  and  Rice  arrays  recorded  the  same 
sources,  any  differences  in  wavefield  between  the  arrays  for  either  SALUT  or  VILLE  is  due 
only  to  the  slight  differences  in  propagation  path  arising  from  the  35  km  separation  of  the 
arrays  and,  peihaps  more  significantly,  to  differences  in  the  geologic  structure  surrounding  the 
arrays.  Conversely,  any  differences  between  SALUT  and  VILLE  for  either  array  site  are  due 
to  differences  in  the  seismic  source  and  the  difference  in  propagation  path  arising  from  the  40 
km  separation  between  the  two  explosions.  These  comparisons  will  be  further  quantified  in  the 
spectral  analysis  which  follows  in  Chapter  3. 

Before  focusing  in  on  the  Pn  array  recordings,  refer  to  Figure  2.5,  which  shows  typical 
recordings  of  both  the  Pn  and  Pg  wavetrains  at  the  Savahia  Moimtain  and  Rice  Valley  arrays. 
Contrasting  explosions,  there  are  two  noticeable  differences.  First  is  the  somewhat  lower  qual¬ 
ity  of  the  VILLE  data.  Due  to  recording  difficulties,  the  first  couple  seconds  of  the  VBLLE  Pn 
wavefield  were  not  recorded  at  Savahia.  Also,  the  signal  to  noise  ratio  is  lower  for  this  smaller 
yield  explosion;  the  Pn  arrival  is  only  barely  visible  above  the  noise  at  Rice  Valley.  Secondly, 
the  time  to  maximum  Pg  amplitude  from  the  time  of  the  Pg  arrival  differs  between  events.  For 
e’f  ample,  at  Rice  Valley  the  SALUT  Pg  wavetrain  increases  to  its  peak  amplitude  in  about  two 
seconds,  while  for  VILLE  the  time  is  doubled,  to  about  four  seconds.  This  behavior  is  similar 
for  the  Savahia  Mountain  recordings  and  may  be  a  result  of  differences  in  source  properties 
and/or  slight  differences  in  propagation  path.  Contrasting  recording  sites,  the  principal 
difference  appears  to  be  that  the  duration  of  relatively  large  amplitudes  is  greater  at  Savahia  for 
both  the  Pn  and  Pg  waves.  For  example,  at  Rice  the  SALUT  Pg  wave  decays  from  its  peak 
amplitude  to  a  fairly  constant  background  level  in  about  8  seconds,  while  at  Savahia  Pg  is  still 
above  the  background  level  after  15  seconds.  Also,  the  peak  amplitudes  are  somewhat  greater 
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at  Savahia,  particularly  for  VILLE.  These  figures  suggest  significant  amplification  and  possibly 
greater  scattering  effects  at  Savahia  Mountain  relative  to  Rice  Valley. 

2.3.1  Pn  Array  Recordings 

The  complete  SALUT  Pn  wavefield  recoided  at  Savahia  is  shown  in  Figure  2.6.  The 
figure  also  includes  the  onset  of  Pg.  Three  more  Pn  waveform  panels  like  this  will  follow,  fol¬ 
lowed  by  four  similar  panels  for  Pg  in  section  2.3.2.  In  all  cases  the  waveform  data  are 
displayed  as  recorded;  no  processing  has  been  done  other  than  trace  alignment  The  traces 
have  been  aligned  using  the  known  source  azimuth  and  reduction  velocities  of  8  km/sec  for  Pn 
and  6  km/sec  for  Pg.  Also,  for  visual  clarity,  each  of  the  panels  has  been  scaled  differently  in 
amplitude.  True  relative  amplitudes  can  be  inferred  from  Figure  2.5.  The  trace  separation  in 
Figure  2.6  is  25  meters,  the  greatest  station  separation  being  about  4  km.  All  145  recordings 
are  shown.  Note  that  there  is  a  ten-station  break  in  recordings  beyond  station  #96.  The  onset 
of  Pn  is  at  about  0.8  seconds  and  Pg  near  10  seconds  in  the  figure.  The  Pn  onset  is  soon  fol¬ 
lowed  by  a  coherent,  larger  amplitude  arrival  approximately  0.5  seconds  later.  This  arrival  is 
more  prominently  displayed  in  Figure  2.7,  which  shows  typical  recordings  of  the  first  4 
seconds  of  the  SALUT  Pn  wave  at  the  Savahia  and  Rice  arrays.  Based  on  the  near-surface 
velocity  at  Pahute  Mesa  (Leonard  and  Johnson,  1987),  this  secondary  arrival  occurs  near  the 
predicted  time  for  pPn.  Alternatively,  Doombos  and  Kvaema  (1987),  examining  a  group  of 
mining  explosions  in  Norway,  have  interpreted  delayed  and  prolonged  regional  Pn  energy  of 
this  kind  as  due  to  scattering  by  topographic  relief  of  the  Moho.  Therefore,  care  must  be  taken 
not  to  mistake  possibly  Moho  scattered  energy  as  due  to  pPn,  as  this  would  result  in  erroneous 
estimates  of  source  depth,  which  is  an  important  parameter  in  seismic  discrimination.  The 
advantage  of  an  array  here  is  that  the  phase  velocity  of  this  delayed  arrival  can  be  measured, 
and  if  it  differs  significantly  from  that  for  Pn  then  Moho  scattering  may  be  indicated. 

The  SALUT  Pn  and  early  Pg  wavefields  at  Rice  Valley  are  shown  in  Figure  2.8.  Here 
the  trace  separation  is  100  meters,  four  times  the  trace  separation  at  Savahia.  Stations  #17  and 
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#49  through  #96  did  not  properly  record  the  SALUT  wavefleld  and  so  are  omitted  from  the 
figure.  Here  the  maximum  separation  of  displayed  waveforms  is  again  about  4  km.  The  onset 
of  Pn  is  at  about  0.5  seconds  in  the  figure,  again  quickly  followed  by  coherent,  larger  ampli¬ 
tude  ground  motion.  The  onset  of  Pg  is  near  1 1  seconds;  notice  that  relative  to  the  SALUT 
recordings  at  Savahia,  the  Pg  wavefield  at  Rice  is  much  more  emergent  (see  also  Fig.  2.5). 
Also  notice  in  Figure  2.7  the  relative  difference  in  frequency  content  between  the  airay  sites, 
indicating  a  relative  site  effect  between  the  two. 

The  VILLE  Pn  and  eariy  Pg  wavefields  at  Savahia  Mountain  are  shown  in  Figure  2.9. 
The  onset  of  Pn  is  absent  due  to  late  recording,  but  would  have  occurred  near  1.0  seconds  in 
the  figure.  The  VILLE  Pn  and  early  Pg  wavefields  at  Rice  Valley  are  shown  in  Figure  2.10. 
Here  94  of  the  96  stations  recorded  useful  data  and  the  maximum  waveform  separation  is 
about  7  km.  The  onset  of  Pn  is  just  barely  discernible  at  about  0.5  seconds.  The  amplitudes 
of  Pn  are  quite  variable  across  the  array.  A  pocket  of  relatively  large  Pn  amplitude  is  apparent 
in  the  lower  left  of  Figure  2.10. 

2.3J2  Pg  Array  Recordings 

The  first  17  seconds  of  the  SALUT  Pg  wavefield  at  Savahia  is  shown  in  Figure  2.11. 
Only  the  first  96  recordings  are  shown;  the  remaining  49  recordings  display  similar  characteris¬ 
tics  and  have  been  omitted  from  the  figure  for  clarity.  As  is  true  for  all  of  the  following  Pg 
panels,  the  first  two  seconds  of  the  Pg  waveforms  shown  overlap  with  approximately  the  last 
two  seconds  of  data  displayed  in  the  corresponding  previous  Pn  data  panel,  with  changes  in 
amplitude  scale,  and  are  aligned  using  a  reduction  velocity  of  6  km/sec.  Beyond  about  three 
seconds  in  Figure  2.11,  the  Pg  wavefield  becomes  quite  complex  and  maintains  significant 
amplitudes  well  into  the  wavetrain.  Relative  to  Pn,  Pg  appears  to  be  richer  in  low-frequency 
energy.  The  Salut  Pg  wavefield  at  Rice  Valley,  shown  in  Figure  2.12,  is  also  quite  complex, 
though  as  noted  in  Figure  2.5,  it  decays  in  amplitude  faster  than  at  Savahia.  Like  Pn,  Pg  at 
Rice  Valley  appears  to  lack  low-frequency  energy  relative  to  Savahia  Mountain. 
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The  first  96  VILLE  Pg  waveforms  recorded  at  Savahia  are  shown  in  Figure  2.13.  Again, 
the  remaining  recordings  display  similar  characteristics  and  have  been  omitted  from  the  figure 
for  clarity.  This  wavefield  appears  similar  in  its  complexity  to  the  SALUT  Pg  wavefield  in 
Figure  2.11  (note  that  Figure  2.13  is  normalized  to  a  lower  amplitude  scale  than  Figure  2.1  Ij. 
The  VILLE  Pg  wavefield  at  Rice  Valley  is  shown  in  Figure  2.14.  Notice  that  the  high- 
amplitude  pocket  we  saw  for  the  Pn  wavefield  is  also  visible  for  the  Pg  wavefield  between 
about  two  and  six  seconds  in  Figure  2.14,  suggesting  a  similar  site  effect. 

2.4  Explosion  HARDIN 

As  described  earlier,  the  explosion  HARDIN  (  mj,  =  5.5)  was  recorded  by  the  Ruby  Val¬ 
ley  array.  Recall  that  the  LLNL  ELKO  station,  located  15  km  away,  also  recorded  this  event 
The  three-component  recordings  of  HARDIN  at  the  LLNL  station  are  shown  in  Figure  2.15. 
The  regional  phases  Pn,  Pg,  Lg  and  longer  period  surface  waves  are  evident  on  these  record¬ 
ings.  A  typical  3-component  recording  of  HARDIN  at  the  Ruby  Valley  array  is  shown  in  Fig¬ 
ure  2.16.  Because  of  the  different  instrument  response,  the  longer  period  surface  waves  are  not 
visible  at  the  array  site.  As  at  the  LLNL  Elko  station,  Pg  and  Lg  are  of  comparable  amplitude. 

Both  the  LLNL  ELKO  site  and  the  array  site  recorded  significant  transverse  motion, 
presumably  generated  by  scattering  and  multipathing.  The  array  recording  shows  a  strong 
amplification  of  the  horizontal  ground  motion  relative  to  the  LLNL  site.  The  amplification  at 
the  array  site  is  most  likely  due  to  the  fact  that  the  array  is  a  soft  rock  site  and  its  position  in 
the  center  of  the  valley  makes  it  more  susceptible  to  complicating  resonance  and  interference 
effects  resulting  from  the  boundaries  of  the  valley.  These  differences  attest  to  the  strong  sensi¬ 
tivity  of  regional  wavefields  on  local  site  effects.  Similar  basin  amplification  effects  on 
regional  phases  have  been  observed  and  modeled  by  Barker  et  al.  (1981). 
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2.4.1  Pn  Array  Recordings 

For  reference,  the  configuration  of  the  Ruby  Valley  array  is  shown  expanded  in  Figure 
2.17  with  station  numbers  indicated.  The  first  four  seconds  of  the  vertical,  radial,  and 
transverse  components  of  the  Pn  wavefield  at  all  twelve  stations  are  shown  in  Figures  2.18a, 
2.18b,  and  2.18c  respectively.  Station  numbers  are  given  to  the  left  of  the  figure.  The 
waveforms  are  manually  aligned  on  the  early  arrivals  on  the  the  vertical  component  and  the 
resulting  array  beam  is  shown  at  the  bottom.  The  LLNL  Elko  recording  is  also  included  for 
comparison;  it  has  been  convolved  with  the  instrument  response  at  Ruby  Vr’tey  to  facilitate 
comparison.  The  first  motion  of  Pn  at  the  array  takes  place  at  about  0.2  seconds  in  the  figure 
and  has  very  low  amplitude.  Like  the  recordings  of  SALUT  and  VILLE  discussed  earlier,  we 
see  a  larger  arrival  about  0.5  seconds  after  the  first  motion,  near  the  expected  pPn  arrival  time. 
The  waveforms  have  been  roughly  grouped  by  similarity  to  more  clearly  show  an  interesting 
variation  in  amplitude,  namely  that  there  is  a  gradual  attenuation  of  amplitude  moving  towards 
the  eastern  side  of  the  array.  This  effect  will  be  examined  further  in  subsequent  chapters.  The 
first  1.5  seconds  or  so  of  Pn  ground  motion  at  the  array  site  is  predominantly  vertical,  unlike 
the  LLNL  recording,  which  has  significant  radial  ground  motion.  This  difference  is  likely  a 
consequence  of  a  higher  velocity  gradient  within  the  low-velocity  valley  sediments  beneath  the 
array  site.  Beyond  about  1.5  seconds  at  the  array  is  an  onset  of  increased  radial  and  transverse 
ground  motion.  The  low-amplitude  array  beams  indicate  that  either  this  motion  is  incoherent, 
and  is  therefore  due  to  random  scattering,  or  is  actually  coherent  but  has  a  much  different 
phase  velocity,  as  would  be  the  case  for  multipathing.  This  also  will  be  explored  in  the  fol¬ 
lowing  chapters. 

1A2  Pg  Array  Recordings 

The  early  portions  of  the  Pg  wavefield  recorded  at  all  twelve  stations  at  the  array  site  are 
shown  in  Figures  2.19a,  2.19b,  2.19c  for  the  vertical,  radial,  and  transverse  components, 
respectively.  Here  there  is  no  obvious  grouping  so  the  waveforms  are  simply  ordered 
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sequentially  by  station  number.  This  time  window  corresponds  to  the  10  to  20  second  window 
in  Figure  2.16.  The  waveforms  have  been  aligned  along  the  known  source -to-receiver  azimuth 
using  a  reducing  velocity  of  6  km/sec.  The  onset  of  Pg  is  at  about  1 .6  seconds  in  the  figure. 
These  waveforms  are  plotted  to  a  different  scale  than  the  Pn  waveforms  shown  earlier.  The 
relative  difference  in  scale  can  be  inferred  from  Figure  2.16  which  shows  Pg  having  a  peak 
amplitude  about  two  times  greater  than  Pn  on  the  horizontal  components.  Of  the  three  com¬ 
ponents,  the  superposition  of  the  vertical  displays  the  most  constructive  interference,  as  was  the 
case  for  Pn.  However,  for  all  components  the  degree  of  overlap  among  the  waveforms  is  not 
great,  as  reflected  in  the  low  amplitude  of  the  waveform  beams. 
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_ Table  1:  NTS  Explosions _ 

shot  GMT  latitude  longitude  depth  (m)  iii^ 

SALUT  1985  163:15:15:0.1  37.248  -116.489  608  5.5 

VILLE  1985  163:17:30:0.1  37.088  -116.084  293  4.4 

HARDIN  1987  120:13:30:0.1  37.233  -116.423  625  5.5 


_ Table  2;  Recording  Sites _ 

site  name  latitude  longitude  SALUT  A  VILLE  A  HARDIN  A 

Ruby  Valley  40.603  -115.191 *  *  *  387 

Rice  Valley  34.000  -114.756  385  355  * 

vahia  Mountain  34.250  -114.588  371  340  * 


*  array  did  not  record  explosion 
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FIG,  2.3  Recording  system  amplitude  response  to  velocity  ground  motir>n  at  each  of  the  array 
sites  and  at  the  Lawrence  Livermore  National  Laboratory  (LLNL)  ELKO  station.  The  responses 
have  been  normalized  for  this  figure. 


ihc  known  source  azimuth  and  a  phase  velocity  of  8  km/sec.  The  onset  of  Pn  is  near  0.8  seconds  and  Pg  near  10  seconds  in  the 
figure.  The  bottom  recording  corresponds  to  station  #1.  the  southern-most  station  in  Figure  2.1.  Subscr^uent  recordings  arc 
spiiccd  in  25  meter  segments.  Following  station  #96  a  ten-.station  break  occurs  where  no  data  were  recorded.  The  maximum  sta¬ 
tion  separation  is  approximately  4  km  between  stations  #1  and  #155. 
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SALUT  Pn  -  Savahia  Mountain  35 
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SALUT  Pn  “  Rice  Valley  48 
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FIG.  2.7  Early  SALUT  Pn  at  station  #35  of  the  Savahia  Mountain  array  and  station  #48  of  the  Rice 
Valley  array.  Note  the  onset  of  the  relatively  large  amplitude  arrival  near  0.7  seconds.  This  arrival  coin¬ 
cides  with  the  expected  arrival  time  of  pPn. 
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SALUT  Pn  at  Rice  Valley 
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nc.  2.8  Pn  and  early  Pg  wavefields  of  SALUT  recorded  at  the  Rice  Valley  array.  The  wavefield  has 
been  aligned  using  the  known  source  azimuth  and  a  phase  velocity  of  8  km/sec.  The  onset  of  Pn  is  near 
0.6  seconds  and  Pg  near  11  seconds  in  the  figure.  The  bottom  recording  corresponds  to  station  #1,  the 
southwestern-most  station  in  Figure  2.1.  Subsequent  recordings  are  spaced  in  100  meter  intervals.  For 
SALUT  only  the  first  48  stations  properly  recorded  the  wavefield,  and  recordings  at  stations  #49 
through  #96  are  therefore  omitted.  The  maximum  station  separation  is  approximately  4  km  between 
stations  #1  and  #48. 
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VILLE  Pn  at  Savahia  Mountain 
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no.  z.y  nti  and  early  yg  wavciicids  ot  viLLh  recorded  at  the  bavahia  Mountain  array.  The  wavericld  has  been  aligned  using 
the  known  source  azimuth  and  a  phase  velocity  of  8  km/sec.  The  early  portion  of  Pn,  approximately  3  seconds  of  it,  was  not 
recorded  at  the  array  due  to  recording  problems.  The  onset  of  Pg  is  near  9  seconds  in  the  figure. 
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VILLE  Pn  at  Rice  Valley 
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FIG.  2.10  Pn  and  early  Pg  wavefields  of  VILLE  recorded  at  the  Rice  Valley  array.  Recordings 
from  all  96  stations  are  shown  here.  The  wavefield  has  been  aligned  using  the  known  source 
azimuth  and  a  phase  velocity  of  8  km/sec.  The  onset  of  Pn  is  near  0.5  seconds  and  Pg  near  9 
seconds  in  the  figure. 
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PIG.  2.11  Portion  of  the  SAI.UT  Pg  wavcfield  recorded  at  Savahia  Mountain.  Only  stations  #1 
through  #96  are  displayed.  The  wavcfield  has  been  aligned  using  the  known  source  azimuth  and 
a  phase  velocity  of  6  km/sec.  The  first  two  seconds  displayed  correspond  to  the  last  two  seconds 
in  Figure  2.6,  however  the  amplitude  scale  in  these  two  figures  differ. 
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FIG.  2.13  Portion  of  the  VILLE  Pg  waveficld  recorded  ai  Savahia  Mountain.  Only  stations  #1  through 
#96  are  displayed.  The  wavefield  has  been  aligned  using  the  known  source  azimuth  and  a  phase  velo¬ 
city  of  6  km/sec.  The  first  two  seconds  displayed  correspond  to  the  last  two  seconds  in  Figure  2.9, 
however  the  amplitude  scale  in  these  two  figures  differ. 
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H(>.  2.14  Portion  of  the  VILLE  Pg  wavcfitid  recorded  at  stations  #1  through  96  at  the  Rice  Valley 
array.  The  wavefield  has  been  aligned  usinj  the  known  source  a/imuth  and  a  phase  velocity  of  6 
km/sec.  The  first  two  seconds  displayed  correspond  to  the  last  two  .seconds  in  Figure  2.10,  however  the 
amplitudes  in  these  two  figures  differ. 


HARDIN  at  LLNL  ELKO  Station 
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2,15  Three-component  broadband  recording  of  HARDIN  at  the  LLNL  seismic  station.  The 
regional  phases  Pn,  Pg,  and  Lg  are  indicated.  Z  =  vertical,  R  =  radial,  and  T  =  transverse. 
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Seconds 

FIG.  2.16  Example  of  a  three-component  recording  of  HARDIN  at  the  Ruby  Valley  array.  The 
regional  phases  Pn,  Pg,  and  Lg  are  indicated,  as  is  the  amplitude  scale  in  digital  counts. 
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Ruby  Valley  Array 


Distance  (km) 

FIG.  2.17  Configuraiion  of  the  Ruby  Valley  airay.  Station  numbers  are  indicated. 
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HARDIN  Pn  Vertical  at  Ruby  Valley 
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FIG.  2.18a  Vertical  component  Pn  recordings  of  HARDIN  to  four  seconds.  Waveforms  recorded 
at  each  station  of  the  Ruby  Valley  array  are  shown,  with  the  station  number  to  the  right  of  each 
trace.  The  traces  are  aligned  on  the  early  Pn  arrival  and  are  plotted  with  correct  relative  ampli¬ 
tudes.  The  superposition  of  the  array  recordings  on  this  Pn  alignment  is  also  shown,  as  is  the 
resul'ing  Pn  beam.  At  the  top  of  the  figure  is  the  corresponding  portion  of  the  LLNL  recording  in 
FIG.  2.14  scaled  to  the  early  Pn  amplitude  at  the  array. 
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FIG.  2.18b  Same  as  FIG.  18a,  but  for  the  radial  component. 
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HARDIN  Pn  Transverse  at  Ruby  Valley 
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FIG.  2.18c  Same  as  FIG.  2.18a,  bul  for  the  transverse  component. 
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HARDIN  Pg  Vertical  at  Ruby  Valley 
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FIG.  2.19a  Vertical  component  Pg  recordings  of  HARDIN.  The  time  window  displayed  coiiesponds  to 
the  10  to  20  second  time  window  in  FIG.  2.16.  Waveforms  recorded  at  each  siauon  of  the  Ruby  Valley 
array  are  shown,  with  the  station  number  to  the  left  of  each  trace.  The  traces  are  aligned  based  on  the 
known  source  azimuth  and  a  phase  velocity  of  6  km/sec.  The  superposition  of  the  array  recordings  is 
also  shown,  as  is  the  resulting  Pg  beam.  At  the  top  of  the  6gure  is  the  corresponding  portion  of  the 
LLNL  recording  in  FIG.  2.15  scaled  to  the  Pg  amplitude  at  the  array. 
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Chapter  3 


Spectral  Amplitude  Characteristics 


3.1  Introduction 

In  this  chapter  we  will  study  the  spectral  amplitude  characteristics  of  the  regional  data 
sets  described  in  the  previous  chapter.  There  arc  a  number  of  reasons  for  doing  this.  First  of 
all,  because  we  arc  concerned  with  high-frequency  monitoring,  we  need  to  know  the  frequency 
limit  of  the  signal.  Secondly,  to  the  extent  that  Pn  and  Pg  can  be  useful  for  discrimination  and 
yield  estimation,  it  is  important  to  be  aware  of  the  variation  in  spectral  amplitude  levels  over 
small  distances  in  receiver  location  which  result  from  the  effects  of  small-scale  variations  in 
geologic  structure.  Knowledge  of  these  stochastic  variations  is  particularly  important  for  a 
prospective  single-station  site  which  may  not  otherwise  be  aware  of  the  spectral  uncertainty 
caused  by  small-scale  spatial  variations.  We  will  therefore  examine  variations  in  spectral 
amplitude  over  the  aperture  of  each  of  our  three  arrays.  Also,  using  spectral  ratios  for  a  fixed 
source,  we  will  look  briefly  at  how  seismic  energy  is  partitioned  between  Pn  and  Pg  and  how 
sensitive  this  partitioning  is  to  slight  changes  in  propagation  path.  Finally,  because  we  have 
recorded  at  both  the  Savahia  Mountain  and  Rice  Valley  arrays  two  explosions  with  nearby 
source  locations,  we  arc  in  a  position  to  remove  the  major  path  effects  using  spectral  ratios  and 
examine  how  SALUT  and  VILLE  source  characteristics  compare  with  a  simple  explosion 
source  model. 

All  spectral  plots  which  follow  were  computed  from  the  waveforms  as  displayed  in  the 
figures  in  Chapter  2.  The  spectra  were  computed  using  the  first  five  37t  tapers  in  the  multiple- 
taper  spectral  estimation  method  with  adaptive  weighting  (Thomson,  1982).  The  principal 
advantage  of  this  method  is  its  ability  to  form  relatively  low-bias,  low-variance  estimates.  A 
description  and  examples  of  this  method  are  given  in  Appendix  A.  For  consistency  with  the 
HARDIN  data  set,  the  SALUT  and  VILLE  waveforms  were  resampled  at  200  samples  per 
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second.  The  Pn  and  Pg  spectral  time  windows  begin  near  the  onset  of  each  phase  and  extend 
for  2,6  seconds  (512  data  points).  The  sampling  in  frequency  is  therefore  0.39  Hz.  The  Pn 
and  Pg  spectra  which  follow  represent  instrument-corrected  displacement  amplitude  spectra 
with  no  correction  for  path  effects. 

3.2  Spectra  at  Savahia  Mountain  and  Rice  Valley 
3.2.1  Pn 

The  spectra  of  the  SALUT  Pn  wavefield  at  Savahia  Mountain  and  Rice  Valley  are  shown 
in  Figure  3.1.  Spectra  from  each  of  the  145  waveforms  in  Figure  2.6  at  Savahia  and  the  47 
waveforms  in  Figure  2.8  at  Rice  are  plotted.  Though  the  number  of  spectra  shown  for  each 
array  site  differs,  recall  that  the  maximum  sensor  separation  is  approximately  the  same,  about  4 
km.  Over  this  distance,  spectral  levels  vary  up  to  an  order  of  magnitude  or  greater  at  both 
array  sites.  As  a  measure  of  the  variation  in  spectral  amplitude  the  coefficient  of  variation  at 
each  frequency  is  shown  beneath  the  spectra.  This  is  given  by  o/p,  the  ratio  of  the  standard 
deviation  of  the  amplitude  spectra,  o,  to  the  mean  value  of  the  spectra,  p.  We  will  be 
interested  to  see  how  o/p  depends  on  frequency,  wavetype,  array  site,  and  seismic  source.  The 
coefficient  of  variation  in  Figure  3,1  is  somewhat  greater  at  Savahia,  though  at  the  lower  fre¬ 
quencies  it  behaves  similarly  at  both  array  sites,  with  local  minimum  in  variation  near  1  Hz 
and  6  Hz.  By  about  30  Hz,  o/p  has  reached  a  value  of  one  at  both  array  sites.  At  higher  fre¬ 
quencies,  o/p  reaches  a  peak  value  of  4.5  near  60  Hz  at  Savahia,  while  at  Rice,  o/p  remains 
below  1.3  up  to  100  Hz.  The  means  of  these  spectra,  computed  assuming  a  log-normal  distri¬ 
bution,  are  shown  in  Figure  3.2.  Also  shown  is  the  the  mean  pre-event  noise  and  same  plot  of 
spectral  variation.  At  Savahia,  we  see,  not  suiprisingly,  that  the  sharp  increase  in  variation 
occurs  right  where  the  signal  runs  into  the  noise,  at  about  25  Hz.  This  is  also  the  case  at  Rice, 
though  the  increase  is  not  as  dramatic.  Beyond  25  Hz,  the  signal  and  noise  spectra,  though 
each  widely  variable,  have  regressed  to  nearly  identical  mean  values.  Refer  ahead  to  Figure 
3.5a  which  superimposes  the  mean  SALUT  Pn  spectra  from  Savahia  and  Rice.  Beyond  6  Hz 
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the  two  mean  signal  spectra  are  very  similar  both  in  amplitude  level  and  spectral  decay,  and 
both  .ave  similar  comer  frequencies  near  6  and  12  Hz.  From  1  to  6  Hz,  however,  the  Pn 
wavefield  at  Rice  is  greatly  depleted  in  energy  relative  to  Savahia,  indicating  significant 
differences  in  site  response  between  these  locations  within  this  frequency  band.  Recall  we  saw 
this  difference  in  frequency  content  in  the  waveform  plots  of  Figure  2.7 

The  spectra  of  the  VILLE  Pn  wavefield  at  Savahia  and  Rice  arc  shown  in  Figure  3.3.  In 
Figure  3.3a,  spectra  from  each  of  the  145  waveforms  in  Figure  2.9  at  Savahia  are  shown.  In 
Figure  3.3b,  only  the  spectra  from  47  of  the  94  waveforms  in  Figure  2.10  at  Rice  are  plotted. 
The  recording  stations  omitted  are  the  same  ones  that  failed  to  properly  record  the  SALUT 
explosion  -  this  is  to  provide  consistency  with  all  of  the  other  spectral  plots,  which  show  spec¬ 
tra  from  waveforms  that  are  separated  by  no  more  than  approximately  4  km.  At  Savahia  the 
spectral  time  window  begins  about  2.5  seconds  into  the  Pn  wavefield  since,  as  shown  in  Figure 
2.9,  the  instruments  did  not  start  recording  until  this  time.  The  mean  spectral  plots  in  Figure 
3.4  show  the  signal  to  noise  ratio  to  be  much  less  for  VILLE,  the  signal  and  noise  merging 
near  13  Hz  at  both  array  sites.  (Recall  that  the  pre-event  noise  at  Savahia  is  not  available; 
instead,  as  an  approximation,  the  mean  noise  spectmm  at  Savahia  for  SALUT  is  plotted.) 
Given  the  similarity  of  propagation  path,  we  would  expect,  at  a  given  array  site,  that  the  spec¬ 
tral  variation  measured  for  VILLE  to  be  similar  to  that  measured  for  SALUT.  This  is  in  fact 
the  case  at  both  arrays  except  below  about  2  Hz  where,  because  of  a  lower  signal-to-noise 
ratio,  o/|i  for  VILLE  has  increased.  Also,  as  we  observed  for  SALUT,  the  spectral  variation  is 
somewhat  larger  at  the  Savahia  array  site.  In  Figure  3.5b  the  Savahia  and  Rice  Ihi  spectra  are 
overlaid  for  VILLE.  As  was  the  case  for  SALUT,  the  Savahia  and  Rice  spectra  are  nearly 
identical  beyond  about  6  Hz  up  until  where  the  signal  and  noise  merge,  though  the  structure  of 
high-frequency  is  different  from  that  for  SALUT.  Again,  below  about  6  Hz,  the  Rice  Valley 
wavefield  is  depleted  in  low-frequency  energy. 
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3.22  Pg 


SALUT  Pg  spectral  plots  from  the  waveforms  recorded  at  Savahia  (see  Fig.  2.11)  and  at 
Rice  (see  Fig.  2.12)  are  shown  in  Figure  3.6.  Again  the  spectral  variation  is  somewhat  greater 
at  Savahia.  Also  at  Savahia,  the  SALUT  Pg  coefficient  of  variation  is  nearly  identical  to  that 
for  the  corresponding  Pn  wavefield,  while  at  Rice  Valley  o/p.  for  Pg  is  consistently  less  than 
that  for  Pn  up  to  about  20  Hz.  The  mean  spectra  are  shown  in  Figure  3.7  along  with  the  pre¬ 
event  noise,  which  for  both  array  sites  merges  with  the  signal  near  20  Hz.  Referring  ahead  to 
Figure  3.10a,  we  see  that,  like  for  Pn,  these  two  mean  spectra  merge  near  6  Hz  and  remain 
nearly  identical  with  increasing  frequency.  We  again  see  comer  frequencies  near  6  and  12  Hz, 
though  here  the  rate  of  decay  between  6  and  12  Hz  is  greater  than  it  was  for  Pn.  Also  similar 
to  Pn,  the  mean  Pg  spectrum  at  Rice  is  relatively  deficient  at  low  frequencies. 

VILLE  Pg  spectral  plots  from  the  148  recordings  at  Savahia  (see  Fig.  2.13)  and  the  first 
47  of  the  96  recordings  at  Rice  (see  Fig.  2.14)  are  shown  in  Figure  3.8.  The  corresponding 
spectral  means  are  displayed  in  Figure  3.9.  The  signal  and  noise  merge  near  13  Hz  at  both 
array  sites.  There  is  only  a  small  difference  in  the  Pg  coefficient  of  variation  between  the 
array  sites.  The  mean  VILLE  Pg  spectra  are  superimposed  in  Figure  3.10b.  We  again  see  the 
lack  of  energy  below  6  Hz  at  Rice  and  an  overlapping  of  the  spectra  beyond.  The  high- 
frequency  decay  shows  the  same  scalloping  near  6  and  8  Hz  that  we  saw  for  Pn.  The  fact  that 
we  do  not  see  this  scalloping  for  SALUT  suggests  it  is  due  to  a  VILLE  source  effect. 

Contrasting  the  spectral  variation  between  Pn  and  Pg  for  VILLE,  at  Savahia  the  Pg  varia¬ 
tion  remains  noticeably  less  than  that  for  Pn  up  to  about  5  Hz,  beyond  which  the  two  arc 
almost  identical.  At  Rice  there  is  no  significant  difference  in  variation  between  Pn  and  Pg 
beyond  2  Hz,  which  is  about  where  the  VILLE  Pn  signal  climbs  out  of  the  noise.  Contrasting 
the  two  sources  SALUT  and  VILLE,  the  Pg  spectral  variations  at  Savahia  are  essentially  ident¬ 
ical  except  between  about  1  and  4  Hz,  where  the  variation  is  up  to  twice  as  great  for  SALUT. 
At  Rice  the  Pg  variation  for  SALUT  is  slightly  less  than  that  for  VILLE  up  to  about  10  Hz, 
beyond  which  the  difference  is  small. 
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The  characteristics  of  the  coefficient  of  variation  can  be  summarized  as  follows: 

(1)  The  spectral  variations  computed  arc  frequency  dependent,  varying  up  to  a  factor  of  two 
over  the  signal  bandwidth.  The  variation  increases  where  the  signal-to-noise  ratio  is  small. 
The  least  amount  of  variation  occurs  typically  between  1  and  2  Hz,  though  comparably  low 
values  often  exist  at  higher  frequencies. 

(2)  Contrasting  the  two  explosions  SALUT  and  VILLE  for  a  fixed  recording  site  (Savahia  or 
Rice)  and  a  fixed  arrival  type  (Pn  or  Pg),  o/p.  is  generally  only  slightly  larger  for  VILLE  up  to 
from  5  to  10  Hz,  beyond  which  o/p  is  about  the  same  for  the  two  explosions.  The  difference 
may  be  due  primarily  to  the  much  lower  signal-to-noise  ratio  for  VILLE. 

(3)  Contrasting  Pn  and  Pg  for  a  fixed  explosion  and  fixed  recording  site,  the  spectral  variation 
of  Pg  is  typically  either  equal  to  or  somewhat  less  than  the  spectral  variation  of  Pn.  Part  of  the 
difference  here  may  be  due  to  the  lower  signal-to-noise  ratio  for  Pn. 

(4)  Contrasting  recording  sites  for  a  fixed  explosion  (SALUT  or  VILLE)  and  a  fixed  arrival 
type  (Pn  r  Pg),  a/p  is  consistently  larger  at  the  Savahia  Mountain  site  over  almost  all  of  the 
signal  bandwidth.  This  suggests  that  decorrelating  effects,  such  as  scattering,  are  more  severe 
at  the  Savahia  site. 

3.3  Spectra  at  Ruby  Valley 

The  Ruby  Valley  spectra  shown  below  were  computed  in  the  same  manner  as  the  Savahia 
and  Rice  spectra  shown  above.  That  is,  the  first  five  37t  tapers  from  the  multi-taper  estimation 
procedure  with  adaptive  weighting  were  used  to  estimate  the  spectra  from  a  2.6  second  win¬ 
dow  of  data  near  the  beginning  of  the  Pn  and  Pg  waveforms,  as  displayed  in  the  figures  in 
Chapter  2.  As  before,  instrument-corrected  displacement  amplitude  spectra  are  plotted  and  the 
sampling  in  frequency  is  0.39  Hz. 
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3.3.1  Pn 


The  three -component  Pn  spectra  at  Ruby  Valley  are  shown  in  Figure  3.11.  The  vertical- 
component  spectral  estimates  come  from  the  2.6  second  window  beginning  near  0.2  seconds  in 
Figure  2.18a.  The  radial  and  transverse  spectral  estimates  come  from  2.6  second  windows 
beginning  near  1.4  seconds  in  Figures  2.18b  and  2.18c.  This  later  window  is  used  because  of 
the  lack  of  energy  earlier  in  the  horizontal  ground  motion.  For  each  component,  the  maximum 
difference  in  spectral  amplitude  among  the  twelve  spectra  is  about  a  factor  of  four,  or  about  12 
db.  The  coefficient  of  variation  is  typically  about  a  factor  of  two  lower  on  the  vertical  com¬ 
ponent  as  compared  to  the  horizontal  components  up  to  about  5  or  6  Hz.  At  higher  frequen¬ 
cies  the  horizontal  variation  is  generally  equal  to  or  slightly  less  than  the  vertical  variation. 
The  difference  in  variation  between  the  two  horizontal  components  is  not  great,  though  it  is  the 
radial  component  which  typically  has  the  lower  variation.  Contrasting  the  vertical-component 
Pn  coefficient  of  variation  here  with  those  shown  earlier  for  the  larger  Savahia  Mountain  and 
Rice  Valley  array  sites,  we  would  expect  the  variation  here  to  be  smaller  because  of  the 
smaller  array  aperture,  1.5  km  as  opposed  to  4  km.  This  is  in  fact  the  case,  but  only  up  to 
about  7  Hz,  beyond  which  the  variation  is  about  the  same.  It  therefore  appears  that  at  the 
higher  frequencies  the  variation  in  spectral  amplitude  is  less  sensitive  to  increases  in  sensor 
separation  than  at  the  lower  frequencies.  The  horizontal  spectral  variation  at  Ruby  Valley  is 
comparable  to  the  vertical-component  variation  of  the  Savahia  and  Rice  recordings. 

The  mean  Pn  spectrum  for  each  component  are  superimposed  in  Figure  3.1  Id.  The  noise 
spectrum  shown  is  that  due  only  to  the  array  site  recording  system.  Pre-event  noise  was  not 
available  due  to  the  insufficient  dynamic  range  of  the  recording  system.  The  differences 
between  the  components  arc  not  great.  The  horizontal  spectra  merge  with  the  system  noise  just 
beyond  20  Hz,  while  the  vertical  component  is  richer  in  high-frequency  energy,  merging  with 
the  noise  near  30  Hz.  The  spectral  notch  near  4  Hz  may  be  a  consequence  of  resonance 
behavior  within  the  valley  sediments.  Refer  ahead  to  Figure  3.13a  where  the  three-component 
Pn  spectra  at  the  nearby  LLNL  ELKO  .station  for  HARDIN  arc  plotted  up  to  the  Nyquist 
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frequency  of  10  Hz.  Here,  unlike  at  the  array  site,  the  transverse  spectrum  is  generally  a  good 
deal  lower  than  tiie  vertical  and  radial  spectra.  As  a  comparison,  in  Figure  3.14a  the  LLNL 
radial  Pn  spectra  is  plotted  with  the  mean  radial  Pn  spectra  at  the  array  site.  With  the  excep¬ 
tion  of  the  3-5  Hz  frequency  band,  the  array  spectrum  is  richer  in  energy  up  to  about  10  Hz, 
no  doubt  due  to  amplification  by  lower  density  sediments  underlying  the  array.  The 
differences  between  the  spectra  between  1  and  3  Hz  can  be  reduced  somewhat  if  a  higher  reso¬ 
lution  taper  is  used;  this  is  shown  in  Figure  A9  of  Appendix  A. 

3.3,2  Pg  Spectra 

Tlie  three-component  Pg  spectra  at  Ruby  Valley  are  shown  in  Figure  3.12.  These  spec¬ 
tral  estimates  all  come  from  the  2.6  second  windows  beginning  near  2.9  seconds  in  Figure 
2.19.  The  range  of  spectral  levels  is  again  about  12  db,  and  the  coefficient  of  variation  here  is 
almost  identical  to  that  for  Pn  on  each  of  the  components.  The  only  notable  exception  is 
below  2  Hz  on  the  radial  component  where  the  Pg  variation  is  about  a  factor  of  two  less  than 
that  for  Pn.  The  vertical-component  variation  here  is  again  somewhat  less  than  that  obser.ed 
at  the  Savahia  and  Rice  arrays,  but  as  for  Pn,  less  so  with  increasing  frequency. 

The  mean  Pg  spectra  for  the  three  components  are  shown  in  Figure  3.12d,  along  with  the 
recording  system  noise  level.  Notice  that  compared  to  Pn  there  is  much  less  structure  to  these 
spectra.  Below  about  4  Hz  the  radial-component  spectrum  has  slightly  larger  values,  though  in 
general  the  three  spectra  are  very  similar  and  can  be  characterized  simply  by  a  low-frequency 
slope  of  about  -2.5  up  to  5  Hz,  then  dropping  sharply  to  about  a  -8  slope  at  higher  frequencies. 
The  signal  merges  with  the  noise  near  20  Hz.  The  corresponding  Pg  spectra  at  the  LLNL  sta¬ 
tion  are  shown  in  Figure  3.13b.  Again  the  transverse  component  is  deficient  in  energy  only  at 
the  lower  frequencies.  In  Figure  3.14b  the  LLNL  radial  Pg  spxictra  is  plotted  with  the  mean 
radial  Pg  spectra  at  the  array  site.  Here,  the  principal  difference  is  the  more  gradual  low- 
frequency  spectral  .slope  of  -1  at  the  LLNL  site.  Any  array  site  amplification  effects  arc 
confined  to  frequencies  below  1  Hz,  and  between  1  and  5  Hz  the  ELKO  site  actually  recorded 
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more  Pg  energy.  The  two  Pg  spectra  merge  between  5  and  10  Hz.  Contrasting  the  spectral 
differences  between  the  LLNL  station  and  the  array  for  Pn  and  Pg,  it  is  clear  that  the  site 
responses  are  very  dependent  on  the  wavetype. 

In  Figure  3.15  is  a  summary  plot  of  all  the  mean  Pn  and  Pg  spectra  from  the  two 
mb  =  5.5  explosions.  In  general,  we  find  the  smallest  differences  occur  at  the  higher  frequen¬ 
cies.  Pn  at  the  two  Mojave  arrays  produced  almost  identical  mean  spectra  beyond  about  6  Hz. 
The  mean  HARDIN  Pn  spectrum,  derived  from  a  very  different  propagation  path,  merges  with 
these  two  Pn  spectra  near  12  Hz.  The  mean  Pg  spectra  from  all  three  array  sites  are  nearly 
equivalent  beyond  about  6  Hz.  If  we  contrast  Pn  and  Pg  at  each  array  site  we  find  much 
greater  differences  in  spectral  amplitude  at  low  frequencies.  The  high-frequency  decay  rate  of 
these  two  wavetypes  is  the  same  beyond  about  10  Hz  at  the  Mojave  arrays  and  beyond  about  7 
Hz  at  the  Ruby  Valley  array.  The  similarity  of  Pn  and  Pg  high-frequency  decay  rates  is  also 
seen  at  the  NORESS  array  (e.g.,  Baumgardt  and  Young,  1990).  By  about  12  Hz  the  difference 
among  all  six  Pn  and  Pg  spectra  become  minimal  and  they  come  to  share  a  common  high- 
frequency  slope  of  about  -8. 

3.4  Pg/Pn  Spectral  Ratios 

Here  we  will  use  the  mean  Pn  and  Pg  amplitude  spectra  shown  above  to  brielly  examine 
the  relative  energy  content  between  the  Pn  and  Pg  wavefields.  We  will  do  this  by  computing 
the  Pg/Pn  power  (or  energy)  spectral  ratios  for  a  given  source-receiver  pair.  Combinations  of 
these  pairs  wilt  tell  us  the  extent  to  which  the  energy  partitioning  between  Pn  and  Pg  is  main¬ 
tained  over  the  separation  distance  of  the  two  array  sites  on  the  receiver  side,  and  over  the 
comparable  separation  distance  of  the  two  explosions  on  the  source  side. 

These  spectral  ratios  can  be  equivalently  thought  of  as  ratios  of  propagation  effects.  In 
general,  once  a  wavcfield  reaches  a  distant  recording  site,  the  source  displacement  spectrum, 
IS(0l.  has  been  modified  by  propagation,  which  includes  the  effects  of  geometric  spreading, 
attenuation,  scattering,  and  station  site  response.  This  can  be  expressed  in  terms  of  an  earth 
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transfer  function  G(f,d.8)  for  the  source,  where  d  is  the  epicentral  distance  and  0  is  the 
source-to-receiver  azimuth.  If  we  also  include  the  signal  distortion  due  to  the  recording  instru¬ 
ment,  1(0,  and  the  ambient  noise,  N(0,  which  we’ll  assume  to  be  uncorrelated  with  the  signal, 
then  the  measured  energy  spectrum  is  given  by 

lA(0P  =  IS(0  G(f,d.e)  1(01^  +  IN(0  l(0l^  (3.1) 

So  if  we  have  two  spectra  with  the  same  source  spectrum  and  instrument  response,  then, 
after  subtracting  out  the  noise  energy,  the  ratio  of  the  measured  spectra  will  be  equivalent  to  a 
ratio  of  earth  transfer  functions.  Therefore,  Pg/Pn  spectral  ratios  derived  from  a  single  source 
and  single  recording  site  display  the  relative  differences  in  propagation  effects  between  the  two 
phases. 

3.4.1  Savahia  Mountain  and  Rice  Valley 

The  Pg/Pn  energy  spectral  ratios  from  the  SALUT  and  VILLE  wavefields  are  shown  in 
Figure  3.16  over  the  frequency  band  for  which  there  was  appreciable  signal  for  both  phases. 
The  ratios  are  based  on  the  mean  spectra  shown  earlier.  On  the  average,  at  both  array  sites,  as 
frequency  increases  the  power  of  Pg  relative  to  Fh  decreases  until  finally,  between  12  to  13 
Hz,  Pn  actually  becomes  the  more  powerful  signal,  though  its  strength  relative  to  Pg  levels  off 
with  increasing  frequency.  For  a  given  explosion,  the  sensitivity  of  the  division  of  energy  to 
the  35  km  difference  in  receiver  location  is  reflected  in  the  differences  between  the  two  ratios 
for  that  explosion.  In  Figure  3.16a  the  Rice  energy  ratio  is  consistently  greater  than  that  for 
Savahia  below  6  Hz,  by  as  much  as  a  factor  of  40.  Referring  back  to  Figure  3.15,  this  is  due 
to  the  fact  that  the  overall  drop  in  signal  amplitude  at  Rice  relative  to  Savahia  is  much  more 
significant  for  Pn  than  for  Pg.  The  energy  balance  between  Pn  and  Pg  is  therefore  not  main¬ 
tained  for  SALUT  at  the  lower  frequencies.  This  is  presumably  due  to  the  fact  that  the  rela¬ 
tionship  between  the  Pn  and  Pg  site  responses  differs  between  array  sites.  Beyond  6  Hz,  the 
differences  between  the  two  array  sites  arc  much  less.  For  VILLE  in  Figure  16b,  the  Rice 
ratio  is  again  greater  at  low  frequencies  but  here  the  greatest  differences  are  confined  between 
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3  and  4  Hz  and  are  no  greater  than  about  a  factor  of  ten.  The  fact  that  the  ratio  discrepancies 
between  array  sites  differ  so  much  between  slightly  different  source  locations  reflects  an 
azimuthal  dependence  of  low-frequency  site  effects. 

For  a  given  recording  site,  the  sensitivity  of  the  energy  ratio  to  the  40  km  difference  in 
source  location  is  also  reflected  in  the  differences  between  the  two  ratios  for  that  recording  site. 
At  Rice  the  SALUT  Pg/Pn  ratio  is  noticeably  greater  than  for  VILLE  up  to  about  6  Hz.  This 
is  also  true  at  Savahia  between  about  2  and  6  Hz.  This  again  may  reflect  an  azimuthal  depen¬ 
dence.  Another  possibility  are  source-related  differences,  such  as  SALUT  being  a  more 
effective  generator  of  Pg  energy  (independent  of  its  larger  yield)  due  to  a  difference  in  proper¬ 
ties  surrounding  the  two  sources.  However  it  is  not  clear  what  these  differences  in  source  pro¬ 
perties  would  be. 

In  general,  over  the  frequency  range  examined,  the  division  of  energy  between  Pn  and  Pg 
is  most  sensitive  to  slight  changes  in  path  at  frequencies  between  approximately  1  and  5  Hz. 
Similar  to  the  findings  on  absolute  spectral  amplitudes,  the  spectral  energy  ratios  become  less 
variable  at  higher  frequencies.  Appendix  B  provides  an  example  of  how  Pn  and  Pg,  as  well  as 
Lg,  compare  in  the  time  domain  within  various  bandpasses. 

3.4.2  Ruby  Valley 

The  Pg/Pn  energy  spectral  ratios  from  the  explosion  HARDIN  are  plotted  in  Figure  3.17. 
Ratios  from  all  three  components  of  ground  motion  at  both  the  array  site  and  the  LLNL  ELKO 
station  are  shown.  The  LLNL  ratio  values  fluctuate  more  widely  and  differ  more  between 
components.  This  is  perhaps  due  to  stochastic  variation  at  the  higher  frequencies,  which  has 
been  reduced  at  the  array  site  due  to  spectral  averaging.  The  overall  behavior  of  the  ratios  at 
the  two  sites  is  similar.  The  exception  is  near  2  Hz,  however  as  mentioned  earlier,  the 
discrepancy  here  can  be  reduced  with  a  higher-resolution  data  taper.  Contrasting  the  HARDIN 
spectral  ratios  with  those  shown  earlier  for  SALUT  and  VILLE,  notice  that  while  the  general 
structure  is  similar,  that  is  decreasing  ratio  values  with  increasing  frequency,  the  magnitude  of 


the  ratio  values  is  typically  much  less.  Whereas  Pn  dominates  Pg  in  spectral  amplitude  begin¬ 
ning  near  6  Hz  for  the  NTS-to-Ruby  Valley  path,  Pn  does  not  exceed  Pg  until  near  12  Hz  for 
the  NTS-to-Mojave  Desert  path.  These  differences  reflect  the  differences  in  P-wave  attenuation 
properties  between  these  northern  and  southern  portions  of  the  Basin  and  Range. 


3.5  VILLE/SALUT  Spectral  Ratios 

In  this  section,  we  will  examine  source  properties  from  Pn  and  Pg  energy  spectral  ratios. 
Spectral  models  of  the  explosive  source  have  been  described  in  a  number  of  studies.  These 
include,  for  example,  work  by  Sharpe  (1942),  Haskell  (1967),  Mueller  &  Murphy  (1971),  von 
Seggem  &  Blandford  (1972),  and  Helberger  &  Hadley  (1981).  Evemden  et  al.,  (1986)  find 
that  the  Sharpe  model  approximates  observed  explosive  source  spectra  from  NTS  reasonably 
well.  In  this  model,  if  we  approximate  the  source  time  function  by  a  step-function  in  pressure, 
then  the  displacement  amplitude  spectrum  at  the  source,  IS(f)l,  of  an  explosion  of  yield  Y  can 
be  parameterized  in  terms  of  a  low-frequency  asymptote  aY,  a  comer  frequency  f^  and  a 
high-frequency  decay  having  a  constant  log-log  slope  of  -2  .  The  relationship  can  be 
expressed  as 


ls(0f  = 


V.  J 


(3.2) 


where  the  constant  a  is  a  function  of  the  material  properties  surrounding  the  explosion. 

If  we  have  two  energy  spectra  lAi(0l^  and  lA2(f)l^  with  the  same  earth  and  instrument 
transfer  functions,  then,  after  subtracting  out  the  noise  energy,  the  ratio  of  the  measured  spectra 
will  be  equivalent  to  a  ratio  of  source  spectra.  For  the  spectral  model  above  we  can  predict 
what  this  is.  If  we  assume  that  the  comer  frequency  can  be  related  to  yield  through 

fc  =  6Y''^  (3.3) 


(e.g.,  Evemden  et  al.,  1986),  where  is  a  constant  which  depends  on  material  properties, 
then  the  ratio  of  energy  spectra  of  two  explosions  having  yields  Yj  and  Y2  is  given  by 
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where  fd  is  the  comer  frequency  corresponding  to  Ai(0. 

Though  not  identical,  the  propagation  paths  for  the  SALUT  and  VILLE  explosions  in  this 
study  are  similar  enough  to  warrant  a  brief  investigation  of  source  spectral  ratios.  Specifically, 
we  can  compute  VILLE/SALUT  energy  spectral  ratios  at  a  fixed  array  site  for  both  the  Pn  and 
Pg  wavefields,  contrast  these  with  the  Sharpe  model  ratios,  and  speculate  on  the  causes  of  any 
differences.  The  mean  amplitude  spectra  used  to  compute  the  observed  energy  spectral  ratios 
are  shown  again  in  Figures  3.18  and  3.19  for  the  Pn  and  Pg  wavefields,  respectively.  The 
noise  corrected  ratios  are  shown  in  Figures  3.20a  and  3.20b.  In  each  figure  the  VILLE  to 
SALUT  energy  spectral  ratio  is  shown  for  both  the  Rice  VaUey  and  Savahia  Mountain  array 
sites  over  the  frequency  band  for  which  there  was  appreciable  signal.  Also  shown  are  the 
theoretical  curves  of  spectral  ratio  generated  from  Equation  (3.4)  above;  adjacent  curves  are 
separated  by  a  factor  of  ten  difference  in  yield  ratio.  Spectral  ratios  for  other  explosive  source 
models  such  as  those  mentioned  earlier,  would  display  similar  behavior,  namely  a  constant 
value  at  low  frequencies,  a  transition  to  larger  values  around  the  comer  frequency,  and  finally 
leveling  off  at  higher  frequencies.  The  comer  frequency  f^  in  (3.4)  is  taken  to  be  1  Hz  for 
SALUT,  which  is  about  that  expected  for  mj,  =  5.5  explosions  at  Pahute  Mesa  (Murphy  et  al. , 
1989). 

The  differences  between  the  observed  spectral  ratios  and  those  predicted  by  the  Sharpe 
model  are  considerable.  Based  on  magnitude- yield  relations  for  NTS  (Evemden,1987),  the 
actual  VILLE/SALUT  yield  ratio  should  be  about  0.03.  Except  for  the  Pg  ratio  at  Savahia 
near  1  Hz,  the  observed  ratios  fall  well  below  the  0.03  contour.  In  fact,  rather  than  displaying 
ratio  values  which  monotonically  increa,se  with  increasing  frequency,  the  observed  ratios 
increase  only  up  to  at  most  a  bit  past  2  Hz.  Beyond  this,  the  ratios  in  general  decrea.se  with 
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increasing  frequency,  reaching  minimal  values  between  about  2  and  4  orders  of  magnitude  less 
than  those  predicted,  depending  on  the  array  site  and  wave  type.  The  difference  between  the 
two  Pn  ratios,  as  well  as  the  two  Pg  ratios,  is  most  likely  due  to  the  small  difference  in  propa¬ 
gation  path  arising  from  the  35  km  separation  of  the  two  array  sites.  The  differences  between 
the  Pg  and  Pn  ratios,  for  example  the  somewhat  greater  Pn  ratio  values  between  about  2  and  6 
Hz,  are  most  likely  due  to  propagation  differences  between  Pn  and  Pg.  For  a  given  array  site, 
these  are  differences  in  earth  transfer  function  ratio  G2 !  Gj,  which,  recall,  are  not  entirely  can¬ 
celed  out  due  to  the  40  km  separation  of  the  two  sources. 

Now  let’s  examine  possible  explanations  for  the  basic  observation  that,  for  each  of  the 
four  computed  ratios,  relative  to  the  Sharpe  model  predictions,  the  VILLE  wavefield  appears 
deficient  in  spectral  energy,  and  increasingly  so  above  1  to  3  Hz,  depending  on  the  array  site 
and  wave  type.  For  our  data  set,  discrepancies  between  the  predicted  ratios  and  those  observed 
should  arise  principally  for  two  reasons:  (1)  source-to-receiver  propagation  path  differences  due 
to  the  40  km  separation  of  the  two  sources  and  (2)  source  depth  effects.  The  fact  that  the 
differences  between  the  spectral  ratios  arising  from  the  35  km  separation  of  recording  sites  are 
relatively  small  argues  against  major  discrepancies  with  the  model  due  to  propagation  path 
differences  arising  from  the  40  km  separation  of  the  two  sources.  With  regard  to  source  depth 
effects,  the  model  does  not  account  for  the  effects  of  spall,  amplitude-dependent  attenuation, 
and  differences  in  rr'''*i''rn  p'^perties  between  two  shot  points.  The  effect  of  spall  is  to  act 
as  a  delayed  secondary  source  which  will  add  energy  to  the  signal  spectrum.  Consequently, 
low  spectral  ratios  values  could  result  over  the  bandwidth  where  the  SALUT  spall  effect  is 
more  pronounced  than  the  VILLE  spall  effect.  However,  for  an  explosion  the  size  and  depth 
of  SALUT,  the  addition  of  significant  spall  energy  is  restricted  primarily  to  frequencies  below 
about  1  Hz  (Taylor  and  Randall,  1989).  Therefore,  at  best,  it  appears  than  spall  can  account 
for  the  low  spectral  ratio  values  only  near  1  Hz  in  Figures  3.20a  and  3.20b.  At  higher  fre¬ 
quencies,  attenuation  effects  arising  from  the  311  m  separation  in  source  depth  may  be  respon¬ 
sible,  as  described  below. 
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In  the  region  surrounding  an  explosion  where  large  strains  result  in  a  non-linear  behavior 
of  the  wavefield,  the  attenuation  may  be  modeled  as  a  linear  function  of  strain,  e, 

Q-*(f,e)  =  Qo''(f)  +  Ye  (3.5) 

(e.g.,  Mavko,  1979,  Minster  and  Day,  1986).  Tire  sensitivity  to  strain  is  determined  by  the 
constant  y.  which  is  a  function  of  rock  type,  microstructure  (crack  density,  porosity,  grain  size), 
confining  pressure,  and  volatile  content,  and  can  range  in  value  from  near  1  to  1000  (Bonner 
et  al. ,  1989).  The  question  here  is  whether  Q“’(f,e)  surrounding  the  VILLE  source  is 
sufficiently  larger  than  that  surrounding  the  SALUT  source  such  that  by  the  time  the 
wavefields  leave  the  source  area,  the  VILLE  wavefield  is  depleted  in  1  to  20  Hz  energy  more 
so  than  SALUT. 

While  the  strain  amplitude,  e,  is  certainly  greater  for  the  the  larger  explosion  SALUT, 
and  would  therefore  contribute  to  greater  attenuation,  this  may  be  over  compensated  by  other 
attenuation  effects.  For  instance,  due  to  the  lower  confining  pressure  at  shallower  depth,  the 
strain-independent  attenuation  Qo(f)"'  experienced  by  the  downgoing  VILLE  wavefield  above 
the  depth  of  SALUT  should  be  greater  than  that  experienced  the  downgoing  SALUT  wavefield. 
In  addition,  the  strain  sensitivity  y  is  larger  for  shallower,  relatively  higher  attenuating,  lower 
modulus  materials,  but  can  be  decreased  greatly  with  increased  confining  pressure.  Also,  being 
above  the  water  table,  the  source  coupling  is  less  for  VILLE  than  for  SALUT,  which  was 
below  the  water  table.  Therefore,  in  addition  to  a  larger  Qo(f)  for  the  VILLE  explosion  we 
should  also  expect  that  y  surrounding  the  VILLE  explosion  to  be  significantly  greater  than  that 
surrounding  the  deeper  SALUT  explosion.  While  these  parameters  are  difficult  to  quantify,  it 
seems  likely  that  the  values  of  Qo(f).  y.  and  e  should  have  resulted  in  a  greater  net  attenuation 
for  VILLE,  which  could  explain  much  of  the  anomalous  spectral  ratio  behavior.  The  indica¬ 
tion  here  is  that  source  depth  appears  to  have  a  significant  effect  on  the  emitted  wavefield  and 
that  discriminates  brsed  on  the  spectral  characteristics  of  explosions  must  bear  in  mind  attenua¬ 
tion  effects  which  may  be  strongly  depth,  as  well  as  yield,  dependent 
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3.6  Summary  Discussion 

In  this  chapter  we  have  examined  a  number  of  spectral  characteristics  of  regional  Pn  and 
Pg  wavefields  recorded  outside  the  Nevada  Test  Site.  The  principal  findings  are  listed  below. 

Signal  Bandwidth 

For  the  mj,  =  5.5  explosions,  SALUT  and  HARDIN,  we  find  the  vertical  component  Pn 
signal  recorded  above  the  noise  to  approximately  30  Hz.  The  Pn  horizontal  components  at  the 
Ruby  Valley  array  merge  with  the  noise  closer  to  20  Hz.  The  Pg  high-frequency  spectral 
amplitude  levels  are  some'  hat  less  than  for  Pn.  The  SALUT  Pg  wavefield  and  all  three  com¬ 
ponents  of  the  HARDIN  Pg  wavefield  reach  the  noise  near  20  Hz.  For  the  smaller  m^  =  4.4 
VILLE  explosion  both  the  Pn  and  Pg  signals  reach  the  noise  near  only  13  Hz.  In  the  context  of 
a  low-yield  test  ban  treaty,  effective  high-frequency  monitoring  of  small-magnitude  events, 
smaller  still  than  VILLE,  will  require  good  signal  to  noise  ratios  -  the  13  Hz  signal  limit  for 
VILLE  suggests  that  we  may  have  to  depend  on  low-noise  borehole  instruments  to  increase  the 
signal  bandwidths  for  small  events. 

Spatial  Variation  in  Spectral  Amplitude 

Over  the  4  km  aperture  of  the  Rice  Valley  and  Savahia  Mountain  arrays,  spectral  ampli¬ 
tude  levels  vary  by  as  much  as  30  db.  However,  it  is  significant  that  at  high  frequencies  the 
mean  of  these  widely  varying  spectra  regress  to  very  similar  means  between  array  sites  (Figure 
3.15).  This  underscores  the  usefulness  of  spatial  averaging  for  reliable  estimates  of  regional 
spectra.  Over  the  1.5  km  aperture  of  the  Ruby  Valley  array  the  maximum  difference  in  spec¬ 
tral  levels  drops  to  about  12  db  for  each  of  the  three  components.  The  means  of  the  three- 
component  spectra  are  almost  identical. 

The  coefficient  of  variation  was  found  to  be  most  consistently  dependent  on  recording 
site  characteristics  rather  than  on  wavetype  or  source  properties.  The  spectral  variation  is 
greater  over  the  larger  aperture  arrays,  but  only  up  to  about  7  Hz,  beyond  which  there  is  little 
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difference  in  variation  between  the  1 ,5  km  and  4  km  arrays.  The  greatest  variation  exists  at 
the  Savahia  Mountain  array  site.  At  the  Ruby  Valley  site,  the  spectral  variation  of  the  horizon¬ 
tal  components  becomes  as  low  or  lower  than  the  vertical  variation  as  the  frequency  increases 
beyond  about  5  Hz. 

If  we  consider  the  spectral  variation  over  the  much  larger  distance  between  array  sites, 
we  find  much  greater  variation  below  6  Hz.  Beyond  6  Hz  the  differences  in  spectra  between 
the  three  array  sites  becomes  very  small. 

Spectral  Energy  Ratios 

The  Pg/Pn  ratios  we  examined  display  the  common  feature  that  the  amount  of  Pg  energy 
relative  to  Pn  decreases,  on  the  average,  with  increasing  frequency.  Pn  eventually  becomes  the 
dominant  signal.  The  magnitude  of  the  ratio  values  is  strongly  dependent  on  the  source-array 
site  combination  with  the  greatest  differences  taking  place  at  the  lower  frequencies.  Finally, 
the  anomalous  behavior  of  the  VILLE/SALUT  spectral  ratios  indicates  that  depth-dependent 
attenuation  effects  on  these  regional  waveflelds  is  significant. 

An  important  point  to  make  from  this  chapter  is  that  variations  in  mean  spectral  ampli¬ 
tude  and  energy  ratios  due  to  array-site  separation  appear  more  stable  at  the  higher  frequencies. 
This  observation  is  somewhat  reassuring  given  that  it  is  the  high-frequency  spectral  characteris¬ 
tics  which  are  particularly  important  in  low-yield  seismic  monitoring. 
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FIG.  2.8  arc  ploilcd.  Each  .specirum  is  computed  from  the  first  2.6  seconds  following  the  onset  of  Pn.  The  sjxxlral  csiimaies  arc 
made  using  the  multiple  3rt  tapering  described  in  Appendix  A.  The  maximum  sensor  separation  for  both  groups  of  spectra  is 
approximately  4  km.  Also  shown  is  the  coefficient  of  variation  for  each  group  of  spectra. 
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FIG.  1.3  Instrument-corrected  Pn  displacement  amplitude  spectral  estimates  for  VILLE  at  the  Savahia  Mountain  and  Rice  Valley 
arrays,  (a)  For  Savahia,  spectra  from  tlic  148  waveforms  in  FIG.  2.9  arc  plotted,  (b)  For  Rice,  spectra  from  the  47  wavefonns  in 
FIG.  2.10  arc  plotted.  Each  spectrum  is  computed  from  the  first  2.6  seconds  following  the  onset  of  Pn.  The  spectral  csiimaics  are 
made  using  the  multiple  3n  tapering  described  in  Appendix  A.  The  maximum  sen.sor  .separation  for  both  groups  of  spectra  is 
approximately  4  km.  Also  shown  in  the  figure  is  the  coefficient  of  variation  of  the  spectra. 
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Also  shown  is  the  coefficient  of  variation  of  the  Pn  spectra,  as  in  Figure  3.3. 


SALUT  Pn 


VILLE  Pn 


FIG.  3.5  Superimposed  mean  instrumem-corrected  Pn  displacement  amplitude  spectra  (reploaed  from 
Figures  3.2  and  3.4). 
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each  station. 
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FIG.  3.10  Superimposed  mean  instrument-corrected  Pg  displacement  amplitude  spectra  (rqplotted  from 
Figures  3.7  and  3.9). 
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coefficient  of  variations  of  the  spectra. 
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HARDIN  Pfl  VERTICAL  AT  RUBY  VALLEY  HARDIN  Pfl  RADIAL  AT  RUBY  VALLEY 

DISPLACEMENT  AMPLITUDE  SPECTRA  DISPLACEMENT  AMPLITUDE  SPECTRA 


o 


72 


coefficient  of  variation  of  the  spectra. 
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HARDIN  Pn  AT  IINL  ELKO  STATION  HARDIN  Pg  AT  LLNL  EUCO  STATION 

DISPUCEMENT  AMPLITUDE  SPECTRA  DISPUCEMENT  AMPLITUDE  SPECTRA 
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RELATIVE  AMPLITUDE 


Mb  =  5.5  Pn  AND  Pg 
DISPLACEMENT  AMPLITUDE  SPECTRA 


HZ 


FIG.  3.15  Mean  m,,  =  5.5  Pn  and  Pg  spectra  recorded  at  the  Savahia  Mountain  and  Rice  Valley  arrays 
(explosion  SALUT)  and  at  the  Ruby  Valley  array  (explosion  HARDIN).  Each  Pg  spectrum  can  be  dis¬ 
tinguished  from  its  corresponding  Pn  spectrum  by  its  larger  low-frequency  amplitudes.  Note  the  simi¬ 
larity  of  high-frequency  decay  among  the  six  spectra. 
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SALUT  Pg  /  Pn  ENERGY  SPECTRAL  RATIO  VILLE  Pg  /  Pn  ENERGY  SPECTRAL  RATIO 
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HARDIN  Pg  /  Pn  ENERGY  SPECTRAL  RATIO  HARDIN  Pg  /  Pn  ENERGY  SPECTRAL  RATIO 

AT  RUBY  VALLEY  ARRAY  AT  LLNL  ELKO  STATION 


Pn  AT  SAVAHIA  MOUNTAIN 


Pn  AT  RICE  VALLEY 


HZ 


FIG.  3.18  Superimposed  mean  SALUT  and  VEXE  Pn  displacement  amplitude  spectra  (replotted  from 
Figures  3.2  and  3.4). 
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Pg  AT  SAVAHIA  MOUNTAIN 


HZ 

FIG.  3.19  Superimposed  mean  SALUT  and  VILLE  Pn  displacement  amplitude  spectra  (replotted  from 
Figures  3.7  and  3.9). 
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Chapter  4 


Pn  and  Pg  Propagation  Characteristics 
at  the  Ruby  Valley  Array 


4.1  Introduction 

The  purpose  of  this  chapter  is  to  study  the  propagation  characteristics  of  the  Hardin  Pn 
and  Pg  wavefields  recorded  at  the  Ruby  Valley  array.  Recall  that  a  major  aim  of  array  pro¬ 
cessing  is  to  use  the  redundancy  of  the  recordings  to  raise  the  signal  to  noise  ratio  of  the 
wavefield.  This  can  be  extremely  important  for  detecting  small-amplitude  signals.  Addition¬ 
ally,  the  relative  phase  delays  of  the  wavefield  as  it  propagates  across  the  array  can  be  used  to 
estimate  phase  velocity  and  source  azimuth.  The  wavetype  can  be  identified  by  its  phase  velo¬ 
city  and,  if  the  locations  of  weapons  testing  sites  are  known,  the  azimuth  can  aid  in  discrimina¬ 
tion.  However  it  is  well  known  that  a  wavefield  undergoes  some  amount  of  decorrelation  as  it 
propagates  through  the  earth.  This  is  particularly  true  of  regional  wavefields,  which  propagate 
largely  through  the  crust  and  are  subjected  to  all  its  structural  complexities.  Scattering,  mul- 
tipathing,  warping,  dispersion,  and  attenuation  of  the  wavefield  due  to  geologic  structure  can 
result  in  amplitude  and  phase  fluctuations  across  the  array.  Wavefield  distortions  of  this  kind 
will  degrade  the  performance  of  an  array.  In  this  chapter  we  will  we  use  frequency- 
wavenumber  array  processing  to  examine  the  array  processing  pei  '  rmance  of  Pn  and  Pg  in 
terms  of  array  gain,  coherence,  phase  velocity,  and  source  azimuth,  paying  particular  attention 
to  the  extent  and  manner  in  which  these  performance  measures  are  affected  by  wavefield  dis¬ 
tortions.  Knowledge  and  anticipation  of  the  array  processing  capabilities  of  the  Pn  and  Pg 
wavefields  bear  upon  the  design  and  siting  of  regional  arrays. 

Also  of  interest  is  the  frequency  limit  to  which  useful  array  processing  can  be  done. 
While  signal  amplitude  may  extend  up  to  30  Hz  and  beyond,  viable  array  processing  is  depen¬ 
dent  on  the  wavefield  being  spatially  coherent  over  the  aperture  of  the  array.  As  frequency 
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increases,  the  wavefield  will  lose  spatial  coherence  over  a  given  distance  as  heterogeneities  are 
sampled  over  more  wavelengths.  At  the  Noress  array  Claassen  (1985)  found  that  significant 
signal  enhancement  is  possible  up  to  20  Hz  for  Pn  and  up  to  10  Hz  for  Pg.  It  is  not  surprising 
that  Pg  loses  signal  correlation  at  lower  frequencies  than  Pn  because  Pg  spends  much  more  of 
its  time  in  the  crust.  The  frequency  limit  for  source-azimuth  estimation  at  NORESS  is  about 
16  Hz  for  Pn  (Bame  et  al. ,  1990)  and  about  8  Hz  for  Pg  (Qaassen,  1985).  However,  a  stable 
geologic  shield,  upon  which  NORESS  is  sited,  is  one  of  the  best  environments  for  signal  pro¬ 
pagation.  In  the  more  tectonically  active  Basin  and  Range  we  should  expect  lower  frequency 
limits  for  these  wavefields.  It  is  worth  noting  that  Bame  et  al.,  (1990)  found  that  the  sub¬ 
configuration  of  the  NORESS  array  that  provides  the  best  regional  azimuth  estimates  is  the 
13-station,  1.5  km  aperture  combination  of  its  center  element  Aq  and  rings  B  and  C.  This 
sub-configuration  is  very  similar  to  the  12-station,  1.5  km  aperture  configuration  of  the  Ruby 
Valley  array. 

We  are  also  concerned  with  the  utility  of  deploying  three-component  sensors.  While  a 
three-component  sensor  alone  may  not  perform  as  well  as  a  single-component  array,  for  exam¬ 
ple  in  terms  of  velocity  and  azimuthal  estimates  (e.g.,  Harris,  1990),  there  can  be  advantages  to 
forming  a  three-component  multi-sensor  array  if  the  horizontal  ground  motion  is  sufficiently 
coherent  (e.g.,  Jurkevic,  1988;  Jepsen  and  Kennet,  1990).  Therefore,  we  will  also  be  examin¬ 
ing  the  frequency-wavenumber  spectra  of  the  radial  and  transverse  components  of  ground 
motion  recorded  at  the  array. 

Before  presenting  results,  we  will  begin  with  a  description  of  the  computational  methods 

used. 

4.2  Array  Power  Estimation 

The  basic  array  processing  tool  we  will  use  here  is  the  so-called  "conventional"  estimate 
of  the  array  power  spectrum.  A  description  of  this  estimate  more  detailed  than  given  here  can 
be  found  in  Abrahamson  and  Bolt  (1987).  The  anay  power  is  based  upon  cross-spectral 
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estimates  for  all  station  pairs  at  an  N-station  array.  If  we  denote  the  frequency  domain 
representation  of  a  time  series  recording  at  the  j*  station  by  Uj(0,  then  the  cross-spectra  foraa 
an  NxN  matrix  S(0.  with  its  elements  given  by 

Sjk(f)  =  E  [Uj(0  Uk(f)*]  ,  (4.1) 

where  E  is  the  expectation  operator  and  *  denotes  the  complex  conjugate.  The  estimate  of  the 
cross-spectral  matrix  will  be  computed  from  a  weighted  average  of  cross-spectra  over  nearby 
frequencies.  If  we  assume  that  the  cross-spectral  values  at  each  frequency  within  the  averaging 
bandwidth  are  statistically  independent  and  identically  distributed,  then  the  error  of  the  estimate 
should  decrease  as  the  averaging  bandwidth  is  widened.  The  validity  of  this  assumption 
degrades  if  there  arc  significant  variations  in  cross-spectral  amplitude  and  phase,  as  would 
likely  be  the  case  with  large  averaging  bandwidths. 

The  conventional  estimate  of  the  frequency-wavenumber  array  power  spectrum  is  given 
by 

P(f,  k)  =  ^  b(k)  S(f)  b(k)  (4.2) 

where  the  overbar  denotes  the  conjugate  transpose.  The  vector  b(k)  is  called  the  beamsteering 
vector  and  represents  the  spatial  phase  delays  which  steer  the  cross-spectra  into  the 

iJc'X' 

wavenumber  vector  k,  i.e.,  bj  =  e  '.  In  computing  P(f,k)  over  a  range  of  k,  we  are  essen¬ 
tially  searching  for  that  value  of  k  which  represents  the  best  plane-wave  fit  to  the  wavefield,  as 
characterized  by  a  maximum  in  the  power  spectrum.  The  value  of  k  at  the  maximum  teUs  us 
both  the  phase  velocity  and  azimuth  of  propagation. 

Theoretically,  the  amplitude  of  the  maximum  array  power  should  be  greater  than  the 
power  level  of  any  of  the  individual  recordings.  A  measure  of  this  is  the  array  signal  gain 
which  is  proportional  to  the  maximum  of  the  total  power  in  the  array  beam  divided  by  the 
average  power  of  the  individual  recordings, 

G(n  =  ~  max  [b(k)  S(0  b(k)] 

trS(0  ^ 

N 
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where  tr  denotes  the  trace  of  the  matrix.  The  gain  is  normalized  so  that  its  maximum  value  is 
1.0.  Zero  correlation  between  all  stations  produces  a  gain  of  1/N;  lower  values  are  possible  if 
stations  pairs  are  negatively  correlated.  Values  of  G(f)  <  1  result  from  anomalous  variations  in 
both  amplitude  and  phase.  We  can  isolate  the  effects  of  phase  alone  by  normalizing  the  cross- 
spectral  matrix.  It  then  becomes  a  matrix  of  intersensor  phase  differences  and  tr  S(f)  =  N. 
When  normalized  in  this  way  equation  (4.3)  effectively  becomes  a  coherence  estimate.  We 
will  refer  to  this  as  the  array  coherence  C(0-  Its  range  of  values  for  positive  correlation  is  also 
1/N  <  C(f)  <1.  Other  methods  of  measuring  signal  enhancement  based  on  signal -to-noise 
ratios  exist  which  we  cannot  examine  here  because  our  instmments  lacked  the  dynamic  range 
to  record  the  pre-event  noise  field.  These  methods  take  advantage  of  the  correlation  structure 
of  the  noise  and  use  station-weighting  schemes  to  improve  the  detection  performance  of  an 
array  (e.g.,  Mykkeltveit  et  ai,  1983;  Ingat  et  ai,  1985;  Der  et  al,  1988;  and  Kvrema, 
1989).  Because  we  are  without  the  advantages  of  a  recorded  noise  field,  the  analysis  here 
repres-^nts  a  lower  bound  on  array  performance. 

An  important  aspect  of  the  frequency-wavenumber  spectrum  is  that  it  is  shaped  by  the 
array  response,  or  beam  pattern,  which  is  detennined  by  the  array  geometry.  This  is  simply 
the  power  spectral  response  to  perfectly  correlated  input  and  describes  the  resolution  and  spa¬ 
tial  aliasing  characteristics  of  the  array.  The  Ruby  Valley  array  response  is  shown  in  Figure 
4.1.  The  resolution  is  characterized  by  the  central  peak  width,  which  is  a  function  of  the  array 
aperture.  The  resolution  improves  (a  narrower  peak)  as  the  aperture  increases.  Spatial  aliasing 
is  a  function  of  the  intersensor  spacing  within  the  array.  The  aliasing  effect  decreases  (a  greater 
bandwith  up  to  the  Nyquist  wavenumber)  as  the  average  intersensor  spacing  decreases. 
Because  the  Ruby  Valley  array  is  two-dimensional  the  Nyquist  wavenumber  varies  with 
azimuth,  but  is  typically  near  about  3  cycles/km.  Our  ability  to  resolve  low  apparent  velocity 
wavefronts,  resulting  perhaps  from  multipathing  or  coherent  local  scattering,  is  limited  to  the 
lower  frequencies.  At  higher  frequencies,  low  velocity  energy  becomes  alia.sed.  For  example, 
for  this  array  a  3  km/sec  wavefront  becomes  alia.sed  above  about  9  Hz.  Higher  velocities  will 
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have  higher  frequency  limits.  A  Pn  wavefront  traveling  at  7  km/sec  can  be  array  processed 
without  aliasing  up  to  about  20  Hz.  In  general,  the  tradeoff  to  constructing  a  large  aperture, 
many-sensor  array  is  that  the  signal  becomes  increasingly  decorrelated  as  the  aperture  increases 
and  the  noise  becomes  increasingly  correlated  as  the  intersensor  spacing  is  decreased  with  the 
addition  of  more  sensors. 

In  the  analysis  of  the  Pn  and  Pg  wavehelds  below,  we  will  be  obtaining  estimates  of 
array  power,  coherence,  phase  velocity,  and  source  azimuth,  all  as  a  function  of  frequency. 
The  power  spectra,  upon  which  all  depend,  are  calculated  as  follows.  First,  to  reduce  the  bias 
due  to  spectral  leakage,  the  time  series  are  pre-whitened.  The  cross-spectral  estimates  are  com¬ 
puted  using  multi-taper  cross-spectral  averaging  (Thomson,  1982).  This  method  of  cross- 
spectral  estimation  and  its  advantages  are  discussed  in  Appendix  C.  Here  we  use  the  five 
lowest-order  4n  tapers  applied  over  the  first  adjacent  frequency  to  either  side  of  each  discrete 
frequency  fn,  for  a  total  of  15  cross-spectra  in  the  average.  For  the  T=1.3  second  time  win¬ 
dows  we  will  be  considering  below  this  corresponds  to  an  averaging  bandwidth  of  about  1.6 
Hz  (fn  ±  1/T).  For  simplicity,  the  spectra  are  computed  in  the  frequency- slowness  (f,s)  domain 
rather  than  the  frequency-wavenumber  (f,k)  domain.  The  two  are  simply  related  by  k=fs. 
Power  spectra  are  first  computed  on  the  raw  data  to  produce  a  preliminary  estimate  of  the 
slowness  propagation  vector.  Then,  in  attempt  to  remove  the  propagation  effects,  the  slowness 
propagation  vector  is  used  to  realign  the  waveforms  to  near  infinite  velocity  and  the  power 
spectra  are  then  recomputed.  Aligning  the  waveforms  removes  extraneous  noise  and  uncorre- 
lated  signal  and  it  generally  results  in  greater  beam  power.  Typically,  after  this  temporal  align¬ 
ment,  slight  propagation  effects  remain  for  some  frequencies,  i.c.,  the  initial  alignment  is  not 
optimal  becau.se  the  revised  spectral  peak  does  not  lie  precisely  at  zero  slowness.  Appropriate 
phase  correctioas  arc  made  in  the  frequency  domain  to  each  cross-spectrum  and  the  power 
spectrum  is  recalculated.  The  process  is  repeated  until  the  spectral  peak  converges  to  zero 
slowne.ss,  usually  requiring  only  one  or  two  iterations.  The  net  slownc.ss  vector  is  simply  the 
sum  of  the  initial  slowness  vector  with  which  the  waveforms  were  aligned  plus  the  sub.scqucnt 
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adjustments.  An  iterative  approach  very  similar  to  this  is  described  by  Wang  and  Kaveh 
(1985).  A  boxcar  weighting  function  for  cross-spectral  averaging  is  used  on  the  aligned 
waveforms.  In  all  calculations,  the  spectra  are  computed  over  a  slowness  grid  sampled  in 
increments  of  0.0025  sec/km.  Finally,  because  phase  velocity  and  azimuth  are  essentially 
based  upon  cross-spectral  phase  information,  the  cross-spectra  from  which  they  are  determined 
are  normalized  prior  to  computing  the  array  power  spectrum  to  remove  the  effects  of  amplimde 
variation.  This  normalization  takes  place  both  before  and  after  the  frequency  averaging  pro¬ 
cess. 

4.3  Pn  Wavefield 

Since  our  concern  is  achieving  high  signal  gain  and  reliable  estimates  of  slowness,  we 
will  concentrate  on  just  the  first  few  seconds  of  the  Pn  and  Pg  wavefields,  avoiding  the  later- 
arriving  coda.  As  we  will  see  below,  the  coherence  quickly  degrades  in  later  time  windows. 
The  3-component  Pn  waveforms  we’ll  examine  are  displayed  in  Figure  4.2  and  ordered  from 
bottom  to  top  in  order  of  increasing  epicentral  distance  from  the  NTS  explosion.  As  we  noted 
in  Chapter  2,  the  earliest  ground  motion  is  essentially  vertical,  indicating  a  high  near-surface 
velocity  gradient.  Significant  horizontal  ground  motion  is  delayed  and  does  not  begin  until 
about  1.3  seconds.  This  delay,  combined  with  the  large  transverse  amplitudes,  suggests  an 
onset  of  scattered  and/or  multi-pathed  energy  beyond  1.3  seconds.  To  examine  this  further  we 
will  divide  the  waveforms  into  two  sequential  T=1.3  second  (256  point)  time  windows, 
denoted  W1  between  0  and  1.3  seconds  and  W2  between  1.3  and  2.6  seconds.  Note  that  even 
in  the  earlier  time  window,  we  can  sec  anomalies  in  amplitude  and  phase.  For  example,  look¬ 
ing  at  the  vertical-component  in  window  1,  in  addition  to  the  eastern  amplitude  attenuation  we 
noted  in  Chapter  2  (stations  4,  5,  6,  11,  and  12),  those  station  pairs  which  lie  at  identical  epi¬ 
central  distances  and  should  therefore  have  no  relative  time  delays  (i.e.,  .stations  6&8,  3&5, 
9&12,  and  lO&ll)  appear  to  have  delays  up  to  a  tenth  of  a  second  or  so.  Tlicsc  anomalies 
will  further  reveal  themselves  in  the  analysis  of  power  spectra  below. 
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4.3.1  Signal  Correlation 

We  will  first  examine  the  signal  correlation  of  Pn  in  terms  of  the  maximum  array  power 
and  array  gain.  The  results  for  windows  W1  and  W2  are  shown  in  Figure  4.3.  No  instrument 
correction  has  been  applied  to  the  power.  Recall  that  the  array  power  is  influenced  by  signal 
strength  and  signal  correlation,  and  that  high  power  does  not  necessarily  translate  into  high 
array  gain.  Array  gain  is  a  measure  of  the  ability  of  the  array  to  pull  a  signal  further  above  the 
background  noise  level,  and  is  theoretically  independent  of  overall  signal  strength.  In  practice 
however,  array  gain  can  can  be  degraded  by  a  low  signal  to  noise  ratio.  Within  window  1 ,  the 
horizontal  power  levels  fall  typically  at  least  two  orders  of  magnitude  below  the  vertical.  In 
contrast,  in  W2  the  horizontal  power  levels  are  up  to  an  order  of  magnitude  greater  than  the 
vertical  at  low  frequencies.  However  beyond  8  Hz  there  is  very  little  difference  in  power 
among  the  three  components  in  W2.  Looking  at  the  array  gain,  in  window  W1  the  vertical- 
component  reaches  the  0.32  mean  gain  level  for  white  noise  near  10  Hz.  This  general  drop  in 
gain  with  increasing  frequency  is  as  expected  for  a  scattering  media  since  inhomogeneities 
become  sampled  over  more  wavelengths.  The  gain  of  the  horizontals  is  significantly  less  than 
the  vertical,  no  doubt  largely  due  to  the  much  lower  signal  to  noise  ratio.  By  window  2  the 
vertical  gain  has  dropped  dramatically,  falling  below  the  noise  level  by  5  Hz.  The  horizontal 
gains  fall  near  or  below  the  noise  level  over  the  entire  bandwidth,  indicating  that  the  increase 
in  horizontal  array  power  is  due  to  signal  strength  and  not  signal  correlation.  This  increase  in 
horizontal  energy  combined  with  a  drop  in  array  gain  relative  to  WI  is  consistent  with  an  onset 
of  increased  scattering  within  this  delayed  time  window. 

The  array  coherence  estimates  for  the  two  windows  are  shown  in  Figure  4.4.  In  general, 
we  should  expect  these  values  to  be  .somewhat  'arger  than  for  array  gain  since  we  arc  ignoring 
the  effects  of  amplitude  variation.  Here,  the  mean  coherence  value  for  white  noi.se  is  approxi¬ 
mately  0.43.  Within  Wl,  as  was  the  case  with  array  gain,  the  array  coherence  is  greatest  on 
the  vertical  and  least  on  the  transverse  at  the  lower  frequencies.  The  vertical  and  radial  coher¬ 
ence  values  are  much  higher  than  the  gain  values  below  5  Hz.  However  beyond  5  Hz  there  is 
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a  rapid  transition  to  low  coherence  values.  The  vertical  gain  in  W2  is  much  lower  and  decays 
more  rapidly  than  in  Wl.  The  coherence  on  both  W2  horizontal  components  has  dropped  to 
extremely  low  values.  Overall,  we  see  a  drop  in  coherence  moving  out  of  window  Wl  and 
into  W2,  again  consistent  with  increased  scattering  in  W2.  As  we  see  below,  this  will  have  an 
effect  on  Pn  phase  velocity  and  source  azimuth  estimates. 

4.3.2  Phase  Velocity  and  Source  Azimuth 

The  Pn  phase  velocity  and  source  azimuth  estimates  we  will  focus  on  here  are  based  only 
on  the  most  coherent  component  of  ground  motion,  i.c.,  the  vertical.  Thev  are  determined 
from  the  slowness  location  of  the  peak  in  the  power  spectra,  computed  with  cross-spectral  nor¬ 
malization.  The  amplitude  of  the  peak  is  given  by  the  magnitude  of  the  coherence  in  Figure 
4.4.  The  results  are  shown  in  Figure  4.5  for  both  time  windows.  In  Wl  the  phase  velocity 
increases  smoothly  with  frequency  up  to  12  Hz,  increasing  from  5  to  7  km/sec.  In  W2  the 
velocity  ranges  between  5  and  6  km/sec  up  to  6  Hz,  beyond  which  the  coherence  is  very  low 
and  the  velocities  arc  for  the  most  part  aliased.  The  Wl  azimuth  estimates  display  a  systematic 
variation,  ranging  from  about  140  degrees  at  the  low  frequencies  to  within  about  5  degrees  of 
the  true  azimu  ’■  of  180  degrees  between  8  and  12  Hz.  Azimuthal  estimates  from  W2  also  fall 
east  of  the  true  azimuth  but  due  to  the  extremely  low  coherence  beyond  4  Hz,  we  cannot  see 
.similar  variations  at  higher  frequencies.  An  analysis  of  the  latter  0.65  seconds  of  Wl  alone  did 
not  reveal  velocities  or  azimuths  significantly  different  from  the  first  0.65  seconds  of  Wl.  This 
suggests  that  the  onset  of  energy  in  this  latter  half  is  due  to  pPn  and  not  Moho  scattering  as 
di.scusscd  in  section  2.3.1. 

The  strong  variation  of  azimuth  with  frequency  in  Wl  is  obviously  not  consistent  with  a 
non-dispersive,  plane-propagating  wavefront.  Also,  the  c.stimalcs  of  Pn  phase  velocity  arc  well 
below  the  7. 6-7. 8  km/scc  Pn  velocities  typical  for  this  region  of  Basin  and  Range  (Hearn 
et  al. ,  1990).  To  examine  this  further.  Figure  4.6  .shows  linearly  interpolated  travel-time  con¬ 
tours  based  on  picks  of  vertical  Pn  arrival  times  at  each  recording  station.  These  contours 
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represent  an  approximation  of  how  the  onset  of  the  Pn  wavefront  propagated  with  time.  The 
contours  show  a  counter-clockwise  warping  of  the  wavefield.  The  uniformity  of  surface  pro¬ 
perties  at  the  array  site  and  the  magnitude  of  the  relative  time  delays  between  stations  of  equal 
epicentral  distance  argue  against  the  warping  being  due  to  variations  in  shallow  soil  response. 
The  source  of  the  warping,  if  it  is  a  local  one,  extends  at  least  a  few  kilometers  beneath  the 
array. 

The  back-azimuthal  direction  of  the  contours  gradually  increases  from  about  140  degrees 
on  the  eastern  side  of  the  array  to  near  180  degrees  on  the  westernmost  side  of  the  array.  This 
is  the  same  range  of  W1  azimuths  we  saw  in  Figure  4.5.  Note  that  the  greatest  amount  of 
counter-clockwise  rotation  takes  place  among  the  attenuated  stations  to  the  east.  The  strong 
azimuthal  variation  with  frequency  arises  from  the  fact  that,  in  addition  to  there  being  an  east- 
west  transition  in  signal  amplitude,  there  is  also  a  transition  in  coherence  across  the  array.  An 
example  of  this  can  be  seen  in  Figure  4.7  which  contrasts  two-station  coherence  for  station 
pairs  7  &  8,  both  located  on  the  western  side  of  the  array,  6  &  12,  both  located  on  the  eastern 
side,  and  6  &  7,  which  overlaps  the  two  sides  (refer  back  to  Figure  2.17  for  the  array 
configuration).  The  station  separation  in  all  cases  is  650  meters.  The  mean  coherence  level 
for  white  noise  is  approximately  0.33.  The  details  of  two-station  coherence  estimation  arc  dis¬ 
cussed  in  the  next  chapter  and  in  Appendix  C.  This  example  shows  that  when  both  stations 
arc  located  on  the  same  side  of  the  array  the  coherence  remains  fairly  high  up  to  about  10  Hz. 
However  when  the  two  sides  arc  taken  together  the  coherence  drops  significantly  with  increas¬ 
ing  frequency,  up  to  about  7  Hz  in  this  example,  with  some  recovery  at  higher  frequencies. 
The  net  result  over  the  entire  array  is  that  as  frequency  increases,  the  array  power  spectra 
become  increasingly  weighted  by  the  greater  number  of  higher-amplitude  stations  to  the  west. 
Because  the  wavefield  becomes  less  warped  towards  the  west,  the  accuracy  of  the  source 
azimuth  .steadily  improves.  This  variation  of  coherence  also  explains  the  increase  in  velocities 
with  frequency.  Looking  back  at  Figure  4.6  we  can  see  that  the  distance  between  contours  gen¬ 
erally  increases,  i.c.,  velocity  increases,  moving  west  across  the  array. 
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One  explanation  for  the  generally  low  values  on  Pn  phase  velocity  could  be  the  effect  of 
a  local  northward -sloping  dip  in  the  basement  of  the  valley  sediments.  For  example,  a  dip  of 
ten  degrees  with  the  velocity  increasing  from  2  km/sec  above  the  dip  to  3  km/sec  below  would 
drop  the  apparent  velocity  from  8  km/sec  to  about  6.5  km/sec.  However  we  do  not  know  if 
such  a  dip  exits.  The  warping  of  the  wavefield  can  perhaps  be  explained  by  a  focusing  effect 
caused  by  the  three-dimensional  structure  of  the  valley.  Similar  effects  of  phase  velocity 
reduction  and  wavefield  rotation  were  seen  by  Chiou  (1991)  who  modeled  the  effects  of  a 
three-dimensional  basin  located  beneath  a  small-aperture  array.  Focusing  might  also  explain 
the  relatively  larger  amplitudes  on  the  western  side  of  the  array.  A  more  complete  understand¬ 
ing  of  the  true  nature  of  the  wavefield  distortion  and  an  anticipation  of  its  effects  would  com¬ 
pensate  to  some  extent  for  the  degraded  source-azimuth  and  phase-velocity  estimation  perfor¬ 
mance  at  this  array  site. 

4.3J  Slowness  Stacking 

A  further  examination  of  the  wavefield  can  be  made  by  displaying  the  complete  power 
spectrum  in  the  slowness  plane.  Rather  than  show  the  spectrum  at  each  frequency,  we  will 
sum,  or  stack,  the  power  spectra  over  frequency  in  the  slowness  plane  to  produce  an  average 
spectrum.  This  is  referred  to  as  slowness  stacking.  Note  that  an  ideal  non-dispersive  pro¬ 
pagating  plane  wave  will  display  a  power  peak  at  the  same  location  in  the  slowness  plane  as  a 
function  of  frequency.  Therefore  if  spectra  from  such  a  wave  are  summed,  or  stacked,  in  the 
slowness  plane,  spectral  peaks  will  sum  constructively,  while  spectral  sidelobes,  which  vary  in 
their  location  with  frequency,  and  random  noise  will  sum  destructively.  A  more  peaked 
response  will  result.  A  slowness  stack  should  provide  a  more  reliable  estimate  of  source 
azimuth  and  phase  velocity,  and  remove  spurious  peaks  in  the  spectrum.  Any  decorrelation  or 
dispersive  qualities  of  the  wavefield  will  show  up  as  a  broadening  or  smearing  of  the  peak  in 
the  stack. 
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Two  frequency  bands  will  be  used  in  the  stacks:  a  low-frequency  band  of  1.6-5. 5  Hz  and 
a  broadband  of  1.6-10.2  Hz.  The  low-frequency  band  encompasses  those  frequencies  where 
we  find  the  greatest  coherence  values.  The  broadband  encompasses  higher  frequencies  with 
low  coherence,  which  nevertheless  might  sum  constructively.  For  comparison  with  the  data 
results.  Figures  4.8a  and  4.8b  show,  in  both  two-  and  three-dimensional  perspectives,  the 
response  of  a  perfectly  correlated  plane  wave  and  the  response  of  a  random  wavefield  stacked 
over  these  frequency  bands.  Note  that  provided  the  signals  are  well  correlated,  the  spectral 
resolution  improves  with  frequency.  Given  what  we’ve  learned  about  the  Pn  wavefield,  the 
best  estimates  of  source  azimuth  and  phase  velocity  would  come  from  a  high-frequency  stack, 
or  a  broadband  stack  which  downweights  the  eastern  stations.  Nevertheless  we  will  display 
simple  low-band  and  broad-band  stacks  to  display  the  effects  of  the  wavefield  distortion.  The 
Pn  low-frequency  stacks  from  W1  are  shown  in  Figure  4.9;  all  three  components  are  included. 
The  spectra  are  plotted  at  slowness  increments  of  0.01  sec/km.  The  stacks  are  based  on  nor¬ 
malized  cross-spectra  and  so  arc  somewhat  more  peaked  than  they  would  be  if  amplitude 
decorrelation  were  included.  The  maximum  power  possible  when  normalized  like  this  is  1.0. 

Within  window  Wl,  the  vertical  and  radial  low-frequency  stacks  result  in  similar  esti¬ 
mates  of  velocity  and  source  azimuth,  however  the  vertical  component  displays  a  somewhat 
more  peaked  response.  The  azimuthal  estimates  are  in  error  by  about  25  degrees,  and  the 
phase  velocity  is  more  typical  of  Pg  than  Pn.  A  power  peak  also  shows  up  on  the  transverse 
component  though  its  amplitude  i.s  very  low  and  its  azimuth  differs  from  the  other  components 
by  about  10  degrees.  The  peak  power  levels  within  W2  (Figure  4.10)  are  all  less  than  those 
for  Wl.  TTie  vertical-component  azimuth  is  comparable  to  that  in  Wl,  however  the  transverse 
power  is  very  low  and  displays  no  coherent  energy  in  the  direction  of  the  source.  There  is 
also  a  greater  difference  between  vertical  and  radial  phase  velocities  and  azimuths.  Isolated 
peaks  representing  coherent  scattering  or  multipathing  arc  difficult  to  identify  within  either  time 
window.  This  suggests  that  the  onset  of  horizontal  ground  motion  and  corresponding  drop  in 
correlation  is  primarily  due  to  ncar-rcccivcr  random  scattering. 


The  broadband  stacks  for  W1  and  W2  are  shown  in  Figures  4.11  and  4.12,  respectively. 
Within  W1  all  peak  power  levels  have  dropped  relative  to  the  low-band  stacks.  The  vertical- 
component  peak,  however,  is  much  closer  to  the  true  source  azimuth  of  180  degrees,  though  a 
good  deal  of  energy  is  still  spread  to  the  east.  The  velocity  of  the  peak  has  also  increased 
from  6.1  km/sec  to  6.7  km/sec.  The  horizontal  stacks,  lacking  high-frequency  coherence, 
display  no  improvement  in  azimuth.  The  vertical  stack  within  W2  displays  a  similar  eastward 
spread,  but  its  peak  is  in  error  by  about  45  degrees,  has  a  much  lower  amplitude,  and  a  velo¬ 
city  of  only  5.3  km/sec.  Both  horizontal  power  levels  within  W2  are  near  the  noise  level  and 
may  not  be  significant.  In  general,  energy  away  from  the  source  direction  within  both  time 
windows  is  spread  throughout  the  slowness  plane  and  of  low  amplitude.  It  therefore  appears 
that  local  scattering  is  not  only  strong,  i.t..  large  signal  amplitude,  but  is  also  quite  random. 
As  we  will  no\/  see  below,  the  situation  is  similar  for  Pg. 

4.4  Pg  Wavefielcl 

Here  again  we  will  restrict  the  analysis  to  the  early  portion  of  the  wavefield.  Figure  4.13 
displays  the  first  2.6  seconds  for  all  three  components.  Unlike  Pn,  with  its  lack  of  early  hor¬ 
izontal  ground  motion,  the  Pg  wavefield  arrives  essentially  isotropic.  There  is  comparable 
ground  motion  on  all  three  components,  with  an  increase  in  overall  amplitude  in  the  latter  half 
of  the  2.6  window.  Given  the  horizontal  nature  of  Pg  propagation,  we  would  expect  a  good 
deal  of  initial  radial  motion.  However  the  initial  transverse  energy  indicates  an  earlier  presence 
of  scattering  than  for  Pn,  some  of  which  derives  from  the  Pn  coda.  We  do  not  see  for  Pg  the 
eastward  attenuation  we  saw  for  Pn.  This  is  presumably  because  by  virtue  of  their  different 
mode  of  propagation,  the  local  sites  effects  between  the  two  wave  types  differ.  As  before,  we 
will  divide  the  waveforms  into  two  sequential  1.3  second  time  windows,  again  denoted  W1  and 
W2.  The  computational  parameters  for  power  spectra  will  be  identical  to  those  used  for  Pn. 
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4.4.1  Signal  Correlation 

The  array  power  and  array  gain  for  the  two  time  windows  are  shown  in  Figure  4.14.  The 
power  is  plotted  with  the  same  relative  scale  as  for  Pn  in  Figure  4.3.  Within  W1  the  three- 
component  power  levels  are  fairly  similar,  in  contrast  to  W1  power  levels  for  Pn.  The  greatest 
variation  in  the  power  is  below  5  Hz,  where  the  vertical  power  is  up  to  about  three  times 
greater  than  the  horizontals.  Within  W2  the  power  levels  are  somewhat  greater  than  in  W1 
and  the  greatest  power  is  still  in  the  vertical  component  at  the  low  frequencies.  The  Pg  array 
gain  in  Wl,  like  Pn,  is  greatest  on  the  vertical,  but  is  lower  in  magnitude  and  drops  to  the 
noise  level  by  6  Hz.  The  horizontal  gains  fall  below  the  noise  level  over  the  entire  bandwidth. 
The  vertical  array  gain  in  W2  also  drops  to  the  noise  level  near  6  Hz,  however  the  gain  values 
are  less  in  this  later  time  window,  though  the  decrease  firom  Wl  is  not  as  dramatic  as  that  for 
Pn.  The  horizontal  gain  levels  remain  below  the  noise. 

The  Pg  array  coherence  estimates  are  shown  in  Figure  4.15.  The  vertical  coherence 
within  window  1  is  as  high  as  it  was  for  Pn  at  the  lowest  frequencies,  but  falls  off  more 
rapidly  with  increasing  frequency.  The  vertical  coherence  drops  in  W2.  There  is  no 
significant  horizontal  coherence  in  either  of  the  two  time  windows.  In  later  time  windows  sig¬ 
nal  correlation  drops  even  further  -  within  about  two  seconds  beyond  W2,  the  vertical  coher¬ 
ence  falls  to  the  mean  noise  level  over  the  entire  bandwidth.  The  principal  difference  between 
the  Pn  and  Pg  array  gains  and  coherence  are  the  much  lower  Pg  values  in  Wl,  particularly  on 
the  horizontal  components.  The  gain  and  coherence  in  window  2  of  Pn  and  Pg  are  comparably 
low  on  all  three  components.  In  terms  of  overall  signal  correlation,  Pn  and  Pg  share  two  basic 
characteristics,  namely  (1)  in  the  presence  of  scattering  the  vertical  component  remains  much 
more  coherent  than  the  horizontal  components  and  (2)  the  decorrelating  effect  of  scattered 
waves  is  greater  in  later  time  windows.  WeTl  now  see  how  they  compare  in  terms  of  phase 
velocity  and  source  azimuth. 
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4.4.2  Phase  Velocity  and  Source  Azimuth 


The  estimates  of  Pg  phase  velocity  and  source  azimuth  are  shown  in  Figure  4.16.  As  for 
Pn,  these  estimates  are  based  on  vertical -component  normalized  cross-spectra.  The  velocities 
in  W1  range  between  5  and  6  km/sec  up  to  about  6  Hz  and  display  an  overall  increase  with 
frequency,  similar  to  the  W1  Pn  velocity  estimates.  In  contrast  to  Pn,  the  maximum  phase 
velocities  of  close  to  6  km/sec  are  quite  reasonable.  Beyond  6  Hz,  where  the  coherence  is 
low,  the  velocity  becomes  aliased.  Though  the  velocities  in  W2  appear  stable  over  a  wider  fre¬ 
quency  band,  this  may  not  be  significant  given  the  low  coherence  values  beginning  near  5  Hz. 
Below  5  Hz,  these  velocities  range  between  about  5.0  and  5.5  km/sec  and  also  increase  with 
frequency.  The  W1  azimuthal  estimates  display  a  low-frequency  variation  similar  to  Pn,  i.e., 
they  generally  increase  with  frequency,  increasing  from  150  to  160  degrees.  Beyond  6  Hz  the 
estimates  are  very  erratic  due  to  the  low  coherence  and  aliased  velocities.  The  open  tiiangles 
denote  azimuths  that  fall  off  the  scale  of  the  plot.  In  contrast  to  Wl,  the  W2  azimuths  fall 
within  only  a  few  degrees  of  the  true  azimuth  below  3  Hz.  Beyond  3  Hz  the  azimuthal  error 
increases  with  frequency,  but  again,  at  5  Hz  and  beyond  the  coherence  is  not  significant. 

Because  of  its  different  mode  of  propagation,  we  would  not  have  necessarily  expected  the 
same  behavior  in  Pg  estimates  of  velocity  and  azimuth  that  we  found  for  Pn.  Nonetheless  the 
behavior  is  similar  within  the  onset  of  Pg  at  low  frequencies  where  the  coherence  is  high. 
However  in  the  later  time  window,  though  the  velocity  remains  similar,  the  azimuthal  estimates 
are  very  different  and  in  fact  improved,  falling  to  within  a  few  degrees  of  the  true  azimuth  at 
low  frequencies.  It  therefore  appears  that  the  onset  of  the  Pn  and  Pg  wavefields  are  subject  to 
similar  site  effects.  However  in  the  later  high-amplitude  Pg  time  window,  by  virtue  of  a 
slightly  different  path  of  propagation,  the  site  effect  has  changed  such  that  relatively  little 
azimuthal  bias  takes  place  at  low  frequencies. 
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4.4J  Slowness  Stacking 


The  low-frequency  Pg  stacks  for  windows  1  and  2  are  shown  in  Figures  4.17  and  4.18. 
Within  Wl,  only  the  vertical  component  displays  a  significant  source-directed  spectral  peak. 
The  azimuthal  estimate  of  156  degrees  is  identical  to  that  for  Pn  in  Wl.  The  Pg  phase  velocity 
estimate  is  5.6  km/sec.  The  horizontal  stacks  have  very  low  amplitude  but  nevertheless  fall 
above  the  0.19  noise  level.  If  significant,  the  low  velocities  and  large  azimuthal  errors  may  be 
indicative  of  locally  generated  P  to  Lg  or  Rg  scattered  energy.  Scattering  of  this  kind  is  also 
suggested  for  some  regional  events  recorded  at  the  Scandanavian  arrays  (Dainty  and  Toskoz, 
1990).  Within  W2  the  vertical  peak  power  levels  have  dropped  somewhat  but  the  phase 
city  and  azimuth  have  changed  very  little.  Both  horizontal  components  have  very  low  power 
and  resemble  the  random  wavefield  spectra.  Contrasting  these  low-frequency  Pg  stacks  with 
those  for  Pn  we  find  for  Pg  lower  power  levels  and  a  lack  of  source-directed  horizontal  energy. 
There  is  a  similar  eastward  bias  of  the  vertical-component  azimuthal  estimates,  however  this 
can  be  removed  for  Pg  by  restricting  the  stacking  in  W2  to  frequencies  under  3  Hz.  The  Pg 
phase  velocities  are  less  than  Pn,  but  not  by  much.  For  example,  the  vertical-component  Pg 
velocity  in  Wl  of  5.6  km/sec  is  very  close  to  the  corresponding  Pn  velocity  in  Wl  of  6.1 
km/sec. 

The  Pg  broadband  stacks  are  shown  in  Figures  4.19  and  4.20.  Within  Wl,  the  vertical- 
component  spectrum  shows  no  improvement  in  estimates  of  velocity  and  azimuth  and  the 
power  is  now  quite  low.  The  horizontals  resemble  the  broadband  random  stack  in  Figure  4.8b. 
The  situation  is  W2  is  similar,  with  the  vertical-component  azimuthal  error  increasing  to  43 
degrees.  The  horizontal  spectra  again  appear  random.  Contrasting  the  broadband  and  low- 
frequency  stacks  for  Pn  and  Pg,  we  find  that  adding  the  higher  frequency  spectra  generally 
degrades  estimates  of  phase  velocity  and  azimuth  due  to  the  low  coherence.  The  exception  is 
the  earliest  time  window  of  Pn,  where  broadband  stacking  improved  the  source  azimuth  esti¬ 
mate  and  resulted  in  a  greater  velocity  contrast  between  Pn  and  Pg,  5.6  vs.  6.7  km/sec. 
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4.5  Summary  Discussion 

In  this  chapter  we  have  examined  and  contrasted  the  array  processing  characteristics  of 
the  three-component  Pn  and  Pg  wavefields  recorded  at  the  Ruby  Valiey  array.  The  principal 
findings  are  summarized  and  discussed  below. 

Scattering  Effects 

The  analysis  has  focussed  on  only  the  fii^t  few  seconds  of  these  wavefields  because  we 
find  that  the  array  processing  performance  is  degraded  at  later  times  due  to  more  severe 
decorrelation  effects.  In  just  a  little  over  a  second  into  the  Pn  wavefield  there  is  an  abrupt 
onset  of  horizontal  ground  motion  and  a  sharp  drop  in  signal  correlation.  Frequency- 
wavenumber  analysis  within  this  later  time  window  was  quite  random  and  did  not  reveal  any 
coherent  sources  of  scattering.  Pg  differs  from  Pn  in  that  it  arrives  at  the  array  with 
incoherent,  high-amplitude  horizontal  ground  motion.  TTiere  is  however  significant  vertical  sig¬ 
nal  correlation  which,  like  Pn,  drops  in  the  later  time  window.  The  fact  that  vertical  correla¬ 
tion  is  sustained  and  horizontal  correlation  is  lost  indicates  that  the  wavefields  arc  being  com 
plicated  more  by  horizontal  structural  irregularities  in  the  crust  than  by  vertical  ones.  Given 
the  horizontal  nature  of  Pg  propagation,  it  is  not  surprising  it  arrives  at  the  array  with  poor  sig¬ 
nal  correlation.  An  additional  decorrclating  effect  for  Pg  is  that  is  superimposed  on  random  Pn 
coda. 

The  onset  and  persistence  of  incoherent  ground  motion  are  likely  due  to  local  effects  such 
as  reverberations  and  wavetype  conversions  in  and  around  Ruby  Valley.  The  strong  effect  of 
the  valley  was  seen  clearly  in  Chapter  2  in  a  comparison  of  the  array  recordings,  located  near 
the  center  of  the  valley,  with  the  LLNL  recording,  located  14  km  to  the  north  and  away  from 
the  interior  of  the  valley  (see  Figures  2.15  and  2.16).  High-amplitude  ground  motion  is  sus¬ 
tained  for  a  much  longer  period  of  time  within  the  valley.  The.se  random  effects  are  in  contrast 
to  more  systematic  peculiarities  which  may  be  related  to  local  site  effects. 
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Site  Effects 


In  this  analysis  we  have  come  upon  irregularities  in  phase  velocity  and  source  azimuth 
which  may  be  related  to  local  site  effects,  though  we  have  made  no  attempt  to  model  these 
effects.  Within  the  onset  of  Pn,  over  the  bandwidth  where  the  array  gain  and  coherence  are 
sufficiently  high,  the  azimuth  varies  systematically  from  eastward-biased  azimuths  40  degrees 
in  error  at  the  low  frequencies,  to  estimates  within  a  few  degrees  of  the  true  source  azimuth  at 
higher  frequencies.  The  Pn  phase  velocity  also  displays  variation,  increasing  slightly  with  fre¬ 
quency.  Still,  at  its  maximum  it  reaches  only  7.0  km/sec,  weU  below  the  expected  7.6-7.8 
km/sec  range.  The  behavior  in  the  later  time  window  is  similar  at  the  lower,  coherent  frequen¬ 
cies.  Estimates  of  average  phase  velocity  and  source  azimuth  from  slowness  stacking  are  com¬ 
plicated  by  this  frequency  dependence  of  the  power  spectra.  This  obviously  poses  a  problem 
for  velocity-based  phase  identification  and  azimuth-based  source  discrimination.  These  irregu¬ 
larities  in  azimuth  and  velocity  may  be  due  to  a  wavefield  distortion  effect  caused  by  the 
three-dimensional  structure  of  Ruby  Valley. 

The  onset  of  Pg  appears  to  be  subject  to  a  site  effect  similar  to  that  for  Pn.  However  in 
the  later  time  window  below  3  Hz  the  azimuthal  bias  drops  greatly,  with  estimates  falling  just 
a  few  degrees  shy  of  the  true  source  azimuth.  This  improvement  in  source  azimuth  is  good  to 
see,  however  it  is  perplexing  and  reflects  the  almost  capricious  nature  of  these  regional 
wavefields. 

Array  Processing  Capability 

The  array  processing  capability  is  of  course  degraded  by  the  random  scattering  and  site 
effects  discussed  above.  In  terms  of  signal  correlation,  the  array  gain  and  coherence  are 
greatest  at  the  onset  of  the  Pn  and  Pg  wavefields,  and  approach  the  mean  noise  level  near  5  or 
6  Hz.  This  is  significantly  less  than  the  frequency  limits  for  signal  enhancement  reported  at 
Noress  where,  for  example,  Qaassen  (1985)  reports  a  Pn  array  coherence  of  nearly  0.9  at  8 
Hz.  The  performance  at  the  Ruby  Valley  site  can  be  improved  to  some  extent  with  more 
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sophisticated  signal  enhancement  methods  which  take  advantage  of  noise  correlation  properties 
and  use  station  weighting  schemes. 

In  terms  of  source  azimuth  and  phase  velocity,  azimuthal  estimates  in  error  by  20  to  30 
degrees  and  Pn  phase  velocities  only  slightly  higher  than  Pg  velocities  were  not  uncommon. 
Nevertheless  there  were  isolated  segments  of  both  the  Pn  and  Pg  wavefields  which  produced 
azimuthal  estimates  only  a  few  degrees  in  error,  though  this  was  within  isolated  bandwidths. 
In  particular  we  had  to  rely  on  the  very  low,  yet  still  significant,  coherence  of  the  Pn  onset 
beyond  5  Hz.  For  later  portions  of  the  Pn  wavefield  and  all  of  the  Pg  wavefield,  velocity  and 
azimuth  estimates  are  not  significant  above  5  or  6  Hz.  It  is  clear  that  effective  routine 
azimuthal-based  source  discrimination  and  velocity-based  phase  identification  at  this  site  would 
require  the  further  experimental  operation  of  an  array  so  as  to  come  to  a  better  understanding 
of  the  local  Pn  and  Pg  azimuthal  and  velocity  bias  effects. 

Three-component  Processing 

The  extent  of  random  wave  propagation  within  Ruby  Valley  severely  limits  the  utility  of 
three-component  processing.  Only  within  the  onset  of  the  Pn  wavefield  is  the  horizontal  array 
gain  and  coherence  significantly  above  the  noise  level.  These  results  indicate  that  three- 
component  signal  processing  schemes  may  be  viable  for  Pn  at  this  array  site.  However,  due  to 
low  coherence  levels,  this  would  have  to  be  limited  to  frequencies  below  about  5  Hz  and  prob¬ 
ably  to  only  the  very  early  portion  of  the  Pn  wavefield. 
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Array  Response 


FIG  4.1  Array  response  of  ihe  Ruby  Valley  array  in  iwo-  and  three-dimensional  perspectives.  The 
power  contours  are  in  linear  10%  intervals  of  the  maximum  power. 
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W1:  Array  Power 


W1:  Array  Gain 


W2:  Array  Power 


W2:  Array  Gain 


FIG.  4  3  3-component  Pn  array  power  and  array  gain  for  Window  1  and  Window  2.  The  vertical  com¬ 
ponent  is  denoted  by  solid  circles,  the  radial  by  solid  squares,  and  the  transverse  by  crosses.  The 
instrument  response  has  not  been  removed  from  the  power.  The  mean  array  gain  level  for  white  noise 
is  approximately  0.32 


Pn 


W1;  Array  Coherence 


W2:  Arrf’y  Coherence 


FIG.  4.4  3-component  Pn  array  coherence  for  Windows  1  and  2.  The  vertical  component  is  denoted  by 
solid  circles,  the  radial  by  solid  squares,  and  the  transverse  by  crosses.  The  mean  coherence  level  for 
white  noise  is  approximately  0.43. 
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Wl;  Phase  Velocity  (km/sec)  W1:  Azimuth  (degrees) 


W2;  Phase  Velocity  (km/sec)  W2:  Azimuth  (degrees) 


FIG.  4.5  Estimates  of  Pn  phase  velocity  and  wavefront  azimuth,  determined  from  the  locations  of 
vertical -component  power  spectral  peaks  in  the  slowness  plane.  Azimuth  refers  to  the  azimuth  to  the 
source.  The  true  source  azimuth  is  180  degrees. 


-0.8 


-0.4 


0.0 

kilometers 


0.4 


0.8 


FIG.  4.6  Pn  amval  time  contours  interpolated  from  manual  picks  of  first  motion  on  the  vertical  com¬ 
ponent.  The  absolute  contour  values  are  times  relative  to  the  earliest  time  of  arrival  at  the  arrav  .  The 
contour  interval  is  0.1  seconds.  Each  element  of  the  array  is  denoted  by  an  asterisk. 
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FIG.  4.7  Two-siaiion  Pn  coherence  from  window  1  for  station  pairs  7  &  8,  both  located  on  the  western 
side  of  the  array,  6  &  12,  both  located  on  the  eastern  side,  and  6  &  7,  one  from  each  side.  The  station 
separation  is  650  meters  in  all  cases. 
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FIG.  4.8a  Slowness  stacks  over  the  frequency  range  1.6  to  5.5  Hz  and  1.6  to  10.2  Hz  for  a  perfectly 
correlated  planer  wavefront.  The  spectrum  is  centered  at  zero  slowness.  The  power  contours  are  in 
linear  intervals  of  0.1. 
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FIG.  4.8b  Frequency-slowness  spectral  stacks  over  the  frequency  range  1.6  to  5.5  Hz  and  1.6  to  10.2  Hz 
for  a  random  wavefield.  The  power  contours  are  in  linear  intervals  of  0.1. 
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FIG.  4.9  Observed  Pn  3-componeni  slowness  stacks  over  the  frequency  range  1.6  to  5.5  Hz  for  Window 
1.  The  peak  power  levels  and  corresponding  estimates  of  phase  velocity  and  azimuth  are  indicated 
above  the  two-dimensional  plots. 


Vertical  Radial  Transverse 


113 


FIG.  4.13  3-componenl  Pg  waveforms  used  in  the  analysis.  The  waveforms  are  ordered  from  bottom  to  top  in  order  of  increasing 
epiccntral  distance  from  the  NTS  explosion.  The  station  numbers  are  indicated  to  the  left  of  the  waveforms.  As  for  Pn,  two  time 
windows  arc  examined;  Window  1:  0  to  1.3  seconds  and  Window  2:  1.3  to  2.6  seconds. 
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W2:  Array  Gain 


FIG.  4.14  3-componenl  Pg  array  power  and  array  gain  for  Window  1  and  Window  2.  TTie  vertical 
component  is  denoted  by  solid  circles,  the  radial  by  solid  squares,  and  the  transverse  by  crosses.  The 
instrument  response  has  not  been  removed  from  the  power.  The  mean  array  gain  level  for  white  noise 
is  approximately  0.32. 
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FIG.  4.15  3-component  Pg  array  coherence  for  Windows  1  and  2.  The  vertical  component  is  denoted 
by  solid  circles,  the  radial  by  solid  squares,  and  the  transverse  by  crosses.  The  mean  coherence  level  for 
white  noise  is  approximately  0.43. 
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FIG.  4.18  Same  as  Figure  17  but  for  Window  2. 


Pg  window  l 
1.6-  10.2  Hz 


VERTICAL 


Pow=0.49  Vel=5.4  Az=l50 


Sx  (sec/km) 


FIG.  4.19  Same  as  Figure  17  with  stacking  over  the  frequency  range  1.6  to  10.2  Hz. 
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Chapter  5 


Coherence  Analysis 


5.1  Introduction 

In  this  chapter  we  will  study  the  spatial  correlation  properties  of  the  Pn  and  Pg 
wavefields  at  the  Ruby  Valley,  Rice  Valley,  and  Savahia  Mountain  array  sites.  It  is  spatial 
correlation  which  primarily  determines  the  effectiveness  of  array  processing  schemes  and  hence 
the  suitability  of  a  location  as  a  potential  array  monitoring  site  for  NTS  explosions.  Correla¬ 
tion  structure  of  course  affects  array  design  directly.  Recording  stations  placed  too  far  apart 
will  have  too  little  correlation.  Stations  placed  very  close  together,  though  having  high  signal 
coherence,  may  also  have  high  noise  coherence.  The  final  design  of  the  NORESSS  array,  for 
example,  was  preceded  by  a  reconnaissance  array  which  studied  signal  and  noise  correlation 
properties  (Mykkeltveit  et  al. ,  1983).  This  chapter  will  assess  the  relative  effectiveness  of  the 
three  array  locations  as  potential  monitoring  sites  in  terms  of  signal  correlation  and  in  doing  so 
provide  a  means  of  constructing  signal  correlation  models  which  can  be  used  to  simulate  array 
processing  characteristics. 

The  approach  here  is  to  first  construct  reliable  estimates  of  coherence  as  a  measure  of 
signal  correlation.  We  then  parameterize  these  estimates  in  terms  of  spatially  continuous 
models  using  a  least  squares  inversion.  Using  these  models  we  can  specify  a  common 
geometry  for  the  three  array  sites  and  directly  contrast  their  array  processing  capability.  The 
data  come  from  the  explosion  HARDIN  recorded  at  the  Ruby  Valley  array  and  the  explosion 
SALUT  recorded  at  the  Rice  Valley  and  Savahia  Mountain  arrays.  The  smaller  explosion 
VILLE  is  not  included  because  of  its  low  signal  to  noise  ratio.  In  the  process  of  the  inversions 
we  will  be  assessing  the  relative  efficacy  of  six  different  models  of  coherence.  The  relative 
performance  of  each  array  will  be  quantified  in  terms  of  array  signal  gain,  though  additional 
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measures  are  possible. 


The  high  station  density  at  Rice  Valley  and  Savahia  Mountain  will  allow  us  to  obtain 
very  statistically  significant  estimates  of  signal  decorrelation  properties.  A  particular  aspect  of 
decorrelation  which  will  be  examined  at  all  array  sites  is  directionality.  Studies  by  Mrazek 
et  al.  (1980)  and  Der  et  al.  (1984)  on  regional  wavefield  correlation  at  large-aperture  arrays 
find  the  Lg  phase  to  decorrelate  more  along  the  wavefront  (transversely)  as  compared  to  along 
the  direction  of  propagation  (longitudinally).  Longitudinal  decorrelation  of  Lg  is  expected 
because  of  dispersion.  However  the  more  significant  transverse  decorrelation  indicates  that 
inhomogeneous  effects  are  more  significant,  such  as  scattering  and  multipathing,  which  causes 
a  spread  in  the  direction  of  arrivals  and  so  decorrelate  the  wavefield  transversely.  Claassen 
(1985),  studying  a  regional  event  at  NORESS,  finds  transverse  decorrelation  to  be  more  severe 
for  all  the  regional  phase  he  examined,  Pn,  Pg,  Sn,  and  Lg.  As  we  will  see,  our  ability  to 
resolve  directional  differences  of  Pn  and  Pg  will  be  largely  dependent  on  array  geometry.  We 
will  precede  the  discussion  of  the  coherence  inversion  results  with  a  section  describing  the 
methods  of  coherence  estimation  and  model  parameterization. 


5.2  Coherence  Estimation  and  Modeling  Procedure 

The  measure  of  signal  correlation  we  will  examine  is  the  magnitude  of  the  complex 
coherence  spectrum.  For  a  pair  of  time  series  the  complex  coherence  is  defined  by 


Yjk(0  = 


SjR(0 

[Sjj(0  Skk(0]''^ 


(-‘5.1) 


where  the  cross-spectrum  Sj|j(f)  is  defined  as  in  Equation  4.3.  In  effect,  given  that  the  complex 
spectral  estimates  have  zero  mean,  the  coherence  is  the  correlation  coefficient  between  Uj(f)  and 
U|j(0*.  Its  absolute  value  lyjijl  is  commonly  used  as  a  measure  of  similarity,  reflecting  the 
extent  to  which  two  time  series  are  linearly  related,  lyjijl  ranges  in  value  from  0  to  1.  Values 
of  lYjkl  <  1  suggest  any  combination  of  the  following: 
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(1)  decorrelation  of  the  wavefield  due  to  scattering  from  heterogeneities  within 
the  earth 

(2)  improper  alignment  of  the  two  time  series;  the  relative  time  delays  from  an 
impulsive  propagating  plane  wave  should  be  removed  to  isolate  the  pulse  and 
concentrate  its  energy  within  the  time  segment  being  analyzed 

(3)  the  presence  of  uncorrelated  noise,  including  recording  system  noise  and 
background  earth  noise 

(4)  dispersion 

As  evident  in  (5.1),  the  coherence  estimates  require  a  procedure  for  estimating  cross-  and 
auto-spectra  of  two  time  series.  Here  again  we  use  the  multiple-taper  spectral  estimation  pro¬ 
cedure  of  Thompson  (1982).  As  in  chapter  4  with  the  frequency-wavenumber  analysis,  we  use 
the  five  lowest-ordcr  47t  tapers  averaged  over  the  first  adjacent  frequency  to  either  side  of  the 
center  frequency.  The  details  and  advantages  of  using  multiple-taper  estimation  are  discussed 
in  Appendix  C.  The  statistical  properties  of  lyjijl  are  rather  complex.  In  particular,  lyj^l  is  not 
normally  distributed  and  it  has  a  variance  which  depends  on  the  magnitude  of  its  value.  How¬ 
ever,  a  transformation  to  tanh'MyjijI  produces  a  distribution  that  is  approximately  normal,  has  a 
constant  variance,  and  a  constant  bias  which  can  be  removed  (Enochson  and  Goodman,  1965). 
For  these  reasons,  a  least  squares  fitting  procedure  is  more  appropriate  for  tanh^'lyji^l  rather 
than  lyjijl.  The  tanh"'lyj,jl  distribution  becomes  a  poor  approximation  for  small  values  of  lyjkl. 
A  more  detailed  discussion  of  the  tanh^'lyji^l  transformation  and  the  corresponding  statistics  is 
given  in  Appendix  C.  The  standard  deviation  of  the  tanh~Myj,^l  estimates,  based  on  our  cross- 
Sjjectral  estimation  procedure,  is  0.26  (Equation  C6  in  Appendix  C). 

We  will  consider  three  types  of  correlation  models:  exponential,  gaussian,  and  .self- 
similar.  Random  scattering  phenomena  described  in  terms  of  exponential  and  gaussian  models 
have  been  u.scd,  for  example,  by  Qiemov  (1960),  Aki  (1973),  Frankcl  and  Clayton  (1986), 
and  Flatte'  and  Wu  (1988).  Frankel  and  Clayton  (1986)  also  examined  a  self-similar  model. 
These  studies  used  these  models  to  characterize  the  correlation  properties  of  the  earth’s  crust. 
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Here,  functions  of  these  types  will  be  used  to  characterize  the  random  nature  of  the  ground 
motion  itself.  For  each  of  the  three  models,  we  will  fit  one-dimensional  (isotropic)  and  two- 
dimensional  (directional)  models.  The  isotropic  models  depend  on  the  absolute  separation 
between  sensors  and  are  characterized  by  a  correlation  distance  a.  The  directional  models 
depend  on  sensor  separation  in  the  longitudinal  and  transverse  directions  and  arc  characterized 
by  orthogonal  correlation  distances  aL  and  ax,  respectively.  The  functional  forms  of  the  lyl 
models  are  given  in  Table  1.  The  self-similar  model  is  derived  from  the  zeroth  order  Bessel 
function  Ko(r).  Because  Kq  equals  infinity  at  r  =  0,  this  function  is  used  to  model  tanh"'  iyi 
directly. 

The  three  isotropic  coherence  models  are  shown  in  Figure  5.1a  for  a  =  1.  Below  r  =  a, 
the  exponential  model  exhibits  the  greatest  decay  rate,  while  beyond  r  =  a,  it  is  the  gaussian 
model  which  decays  most  rapidly.  The  self-similar  model  is  a  compromi.se  between  these  two 
models,  behaving  most  like  the  gaussian  model  below  r  =  a,  and  most  like  the  exponential 
model  beyond  r  =  a.  If  we  assume  unit  amplitude  of  the  wavcfield  the  complex  coherence  in 
Equation  (5.1)  becomes  equivalent  to  the  cross-spectrum  and  we  can  use  Equation  (4.3)  to 
compute  the  corresponding  wavenumber  power  .spectrum  for  each  model.  These  are  plotted  in 
Figures  5.1b  and  5.1c  to  display  the  differences  in  decay  rate  and  in  Figure  5.2a  to  display  the 
full  three-dimensional  perspective.  Were  there  no  decorrelation  of  the  waveficld  the  power 
spectrum  would  simply  be  a  spike  located  at  the  wavenumber  of  propagation.  However, 
effects  such  as  scattering  produce  a  continuum  of  propagation  velocities  and  directions  result¬ 
ing  in  a  bandpass  of  power  in  k  space  centered  at  the  dominant  wavenumber.  The  bandpass  of 
the  one-dimensional  models  is  isotropic  in  k  .space,  i.e.,  decorrelation  effects  arc  uniform  in  all 
directions.  The  two-dimensional  models  allow  for  non-uniform  decorrelation  about  the  central 
wavenumber.  In  terms  of  array  processing,  transverse  decorrelation  will  tend  to  degrade 
azimuthal  c.stimatcs  while  longitudinal  decorrelation  will  tend  to  degrade  estimates  of  phase 
velocity.  Examples  of  the  directional  coherence  models  arc  shown  in  Figures  5.2b  for  equal 
decay  constanus  and  in  Figure  5.2c  for  decay  constants  differing  by  a  factor  of  two.  Note  that 
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the  isotropic  exponential  and  self-similar  models  are  not  equivalent  to  the  corresponding  direc¬ 
tional  models  with  equal  decay  constants,  while  the  gaussian  model  is. 

The  inversion  for  the  model  parameters  from  tanh“'  lyl  is  a  non-linear  one  and  here  we  a 
Levenberg-Marquart  inversion  procedure  (e.g.  Lawson  and  Hanson,  1974).  The  inversion 
requires  an  initial  estimate  of  the  decay  constants,  which  for  the  exponential  and  gaussian 
models  we  can  get  by  performing  a  linear  inversion  with  values  of  In  (lyl).  By  talcing  the 
natural  logarithm  the  exponential  and  gaussian  models  become  linear  in  1/a  and  1/a^,  respec¬ 
tively.  The  initial  gaussian  decay  constants  are  also  used  as  starting  values  for  the  self-similar 
inversion,  which  proved  sufficient  in  synthetic  tests.  The  non-linear  inversions  are  iterated 
until  no  further  significant  improvement  in  fit  can  be  obtained. 

Synthetic  random  noise  values  of  tanh~Myj|jl  at  each  array  are  shown  in  Figure  5.3.  The 
values  are  plotted  as  a  function  of  absolute  station-pair  separation  and  the  best-fitting  isotropic 
model  is  superimposed.  The  number  of  station  pairs  at  each  array  site  is  66  at  Ruby  Valley, 
1081  at  Rice  Valley,  and  10440  at  Savahia  Mountain.  The  90%  significance  level  for 
tanh'MyjijI  noise,  as  derived  in  Appendix  C,  is  approximately  0.7,  which  appears  consistent 
with  the  random  simulations.  The  noise  will  have  its  greatest  effect  on  inversions  for  low 
values  of  decay  constant,  tending  to  return  values  with  a  positive  bias.  An  additional  biasing 
effect  is  that  at  low  coherence  the  tanh“Myj,jl  transformation  becomes  one-sided  as  it  nears  its 
lower  limit  of  zero.  One  course  of  action  is  to  omit  from  the  inversion,  or  downweight, 
tanh'Myj^l  values  at  distances  beyond  which  they  fall  almost  exclusively  below  the  noise  level. 
However,  the  choice  of  a  cutoff  distance  is  not  necessarily  an  obvious  one,  as  wv  will  see.  In 
the  inversions  here  we  will  simply  weight  the  tanh^'lyji^l  values  equally  over  the  full  aperture 
of  the  arrays,  and  keep  the  biasing  effects  in  mind  when  we  make  use  of  the  resulting  models 
later  on. 

The  ability  of  the  tanh“'lyl  data  to  reliably  estimate  the  model  parameters  is  largely 
dependent  on  the  distribution  of  the  independent  variables,  i.e.,  the  spatial  separations.  Note  in 
Figure  5.3  that  differences  in  the  best-fitting  i.sotropic  decay  rates  arise  out  of  differences  in 


125 


station  distribution.  In  general,  estimates  of  the  decay  constant  improve  as  the  number  and 
range  of  the  spatial  separations  increase.  In  Figures  5.4,  5.5,  and  5.6  we  examine  the  disuibu- 
tion  of  spatial  separations  at  each  array  site.  The  distribution  of  the  intersensor  separations  at 
the  Ruby  Valley  array  shown  in  Figure  5.4  is  based  on  a  160  degree  back-azimuth  to  the 
HARDIN  event.  For  each  of  the  66  station  pairs  Figure  5.4a  displays  the  longitudinal  and 
transverse  components  of  intersensor  separation.  In  Figure  5.4b  is  a  histogram  of  the  number 
of  station  pairs  vs.  the  absolute  separation  between  them.  The  histogram  interval  is  0.2  km. 
Histograms  are  also  shown  in  Figures  5.4c  and  5.4d  for  the  components  of  longitudinal  and 
transverse  separation;  the  similarity  of  the  distributions  indicates  that  the  two  directional  decay 
constants  will  be  estimated  with  about  equal  uncertainty.  Similar  distribution  plots  for  Rice 
VaUey  array  are  shown  in  Figure  5.5  based  on  the  known  back-azimuth  to  SALUT.  Here  the 
range  of  longitudinal  separation  is  almost  a  factor  of  four  less  than  the  range  of  tra^'^verse 
separations  and  so  we  would  expect  greater  uncertainty  in  the  estimate  of  the  longitudinal 
decay  constant.  At  Savahia  Mountain  in  Figure  5.6  the  span  of  longitudinal  separation  is  again 
greater  than  the  span  of  transverse  separation,  though  the  difference  is  not  as  great  as  at  Rice 
Valley. 

The  effect  of  station  distribution  on  the  ability  to  resolve  two-dimensional  model  parame¬ 
ters  is  shown  in  Figures  5.7,  5.8,  5.9,  which  display  predicted  values  of  tanh^'lyji^l  at  each 
array  site  based  on  a  two-dimensional  exponential  model  with  decay  constants  differing  from 
between  a  factor  of  2  and  15.  In  anticipation  of  the  results  below  we  have  set  aL  >  ax  at  Ruby 
Valley  and  Savahia  Mountain  and  a^  <  aj  at  Rice  Valley.  The  propagation  directions  used  are 
the  same  as  those  used  in  the  inversions.  Best-fitting  isotropic  exponential  models  arc  also 
shown.  We  sec  that  as  the  decay  constant  ratio  increases  so  does  the  spread  of  tanh^'lyji^l 
values.  The  spread  would  change  very  little  had  we  chosen  larger  values  of  and  ax  with  the 
same  ratio.  The  important  point  to  make  here  is  that  for  small  decay  constant  ratios  the  varia¬ 
tion  in  tanh"'lYj,jl  values  away  from  a  one-dimensional  model  may  be  within  the  standard  devi¬ 
ation  of  the  data,  in  which  case  it  may  not  be  possible  to  resolve  the  true  directionality  of  the 
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wavefield  correlation.  This  is  particularly  true  for  Savahia  Mountain  where  the  spread  of 
values  is  almost  iiupciceptible  even  with  a  factor  of  ii  difference  in  uie  two  decay  constants. 
This  is  a  result  of  the  narrowness  of  the  intersensor  distribution  in  Figure  5.6a.  These  points 
about  the  effects  of  station  distribution,  as  well  as  the  bias  effects  caused  by  noise,  must  be 
kept  in  mind  when  interpreting  the  inversion  results  below. 

53  Pn  Coherence 

As  we  saw  in  Chapter  4,  the  choice  of  time  window  can  have  a  strong  effect  on  signal 
correlation.  If  the  window  is  too  long  the  coherence  can  be  greatly  reduced.  But  a  very  short 
window  will  result  in  poor  frequency  resolution.  To  reduce  the  contaminating  effects  of 
delayed  scattering,  and  as  an  acceptable  degree  of  frequency  resolution,  we  will  (as  we  did  in 
Chapter  4)  restrict  the  Pn  window  length  to  the  first  1.3  seconds  (Af  =  0.78  Hz)  following  the 
onset  of  Pn  at  each  of  the  three  array  sites.  The  recordings  are  aligned  prior  to  windowing. 

The  tanh'^Yj^l  estimates  used  in  the  inversions  are  displayed  in  Figures  Dl,  D2,  and  D3 
in  Appendix  D  for  each  array  site  at  frequencies  ranging  from  1.6  to  14.9  Hz;  the  best-fitting 
isotropic  models  are  also  plotted.  In  general,  by  15  Hz  there  appears  to  no  significant  coher¬ 
ence  beyond  a  station  separation  of  about  5(X)  meters.  Vernon  and  Fletcher  (1990),  in  a  coher¬ 
ence  study  of  local  earthquakes  at  Pinon  Flats,  California  also  find  P  waves  to  be  incoherent 
beyond  about  15  Hz  for  station  separations  over  500  meters.  We  will  make  relative  comparis¬ 
ons  of  the  data  in  Appendix  D  through  the  model  parameters  discussed  below. 

We’ll  preface  the  examination  of  the  the  model  parameters  by  referring  to  Table  5.2 
which  describes  the  relative  goodness  of  fit  for  each  of  the  six  models.  This  is  listed  at  each 
frequency  in  terms  of  the  F  statistic,  which  provides  a  relative  comparison  of  the  ability  of  the 
models  to  fit  the  data.  The  F  statistic  is  simply  the  ratio  of  the  variance  of  fit  of  each  model  to 
the  variance  of  fit  to  the  best-fitting  model.  Larger  F  statistics  represent  poorer  fits  to  the 
model.  The  best  fitting  models  are  indicated  by  an  asterisk  *;  in  some  cases  two  models  fit  the 
data  equally  well.  At  all  three  array  sites  it  is  generally  the  exponential  models  which  provide 
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the  best,  or  near-best,  fit  to  the  data.  At  Ruby  Valley  the  variance  of  the  directional  exponen¬ 
tial  model  is  about  a  factor  of  two  smaller  than  the  isotropic  variance  at  low  frequencies,  but  as 
frequency  increases  the  difference  in  fit  between  the  two  models  becomes  insignificant.  We 
see  similar  behavior  in  the  self-similar  and  gaussian  models  at  Ruby  Valley.  At  Rice  Valley 
the  difference  in  fit  between  the  isotropic  and  directional  exponential  models  is  very  small  over 
the  entire  frequency  band.  This  is  also  the  case  at  Savahia  Mountain.  Some  of  the  poorest 
model  fits  occur  at  Savahia  from  the  self-similar  and  gaussian  models.  The  actual  standard 
deviations  of  fit  for  the  best-fitting  models  arc  shown  in  Figure  5.10,  and  typically  fall  between 
about  0.3  and  0.5  tanh“M7j,jl  units. 

With  this  background  in  mind,  we  can  now  turn  to  Table  5.3  which  lists  the  Pn  decay 
constants  obtained  at  the  three  array  sites  for  each  of  the  six  models.  Because  of  their 
superiority  of  fit,  we  will  focus  the  discussion  primarily  on  the  exponential  models.  It  is 
extremely  unlikely  that  the  larger  decay  constants,  up  to  59  km  for  the  isotropic  model  at  Ruby 
Valley,  can  be  extrapolated  to  such  great  distances.  We  can  only  assert  that  they  arc  appropri¬ 
ate  over  the  aperture  of  the  array  from  which  they  are  derived.  In  general,  we  should  expect 
the  decay  constants  to  decrease  with  increasing  frequency  since  a  higher  frequency  wave 
raverses  subsurface  heterogeneities  over  more  wavelengths.  Referring  to  Figure  5.11a,  which 
7lots  the  exponential  isotropic  decay  constants  for  each  array  site,  we  see  that  this  is  generally 
he  case.  However,  there  are  pronounced  local  maxima  in  decay  constant  values  near  4  Hz  at 
Rice  Valley  and  near  10  Hz  at  Ruby  Valley,  which  are  probably  due  to  complex  site  effects. 
The  increase  near  10  Hz  at  Rice  Valley  and  Savahia  Mountain  is  due  to  noise  bias  mentioned 
earlier  (see  Figures  D.2  and  D.3  in  Appendix  D). 

An  alternative  parameterization  of  isotropic  coherence  is  an  exponential  one  used  by 
Menke  (1990),  which  has  the  advantage  of  also  incorporating  frequency.  This  model  is  given 
by  where  d  is  a  constant.  Equating  this  to  our  lyl=e“^'“  model,  and  using  the  fact  that 

f  =  v/X  where  v  is  velocity  and  X  is  wavelength,  implies  that  our  decay  constant  a  is  a  linear 
function  of  wavelength  and  that  the  quantity  a*f  is  a  con.stant.  fo  examine  the  extent  to  which 
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such  a  model  would  be  applicable  to  the  data  set  here.  Figure  5.11b  displays  a*f  for  Pn  at  each 
array  site.  The  strong  variation  over  frequency  at  each  array  site  indicates  that  the  relationship 
between  the  decay  constants  and  wavelength  is  more  complicated  that  a  linear  one  and  is  site 
dependent.  It  is  interesting  to  note  that  beyond  5  Hz  the  behavior  of  a*f  is  almost  identical 
between  Rice  Valley  and  Savahia  Mountain,  though  they  differ  by  a  factor  of  two. 

The  two-dimensional  exponential  decay  constants  are  more  closely  examined  in  Figure 
5.12.  For  each  frequency  the  ratio  of  the  two  decay  constants  is  plotted  such  that  the  result  is 
greater  than  one.  So  if  aL  <  a-j-  the  ratio  a^/aL  is  plotted  (the  cross-hatched  bars  in  the  figure) 
and  if  aL  >  at.  's  plotted  (the  blank  bars).  The  ratios  therefore  indicate  the  factor  by 
which  the  two  decay  constants  differ.  For  random  noise  the  ratio  is  within  a  factor  of  two.  At 
Ruby  Valley  the  difference  in  decay  constants  is  large,  with  aL  between  about  9  to  18  times 
greater  than  a-]  (i.e..  much  greater  decorrelation  along  the  wavefront).  The  difference  in  the  two 
decay  constants  eventually  decreases  with  increasing  frequency,  as  we  would  expect  as  we  gra¬ 
dually  fall  into  the  noise.  Beyond  9  Hz  the  ratios  lie  just  above  a  value  of  one  and  the  i.sotro- 
pic  model  fiLs  the  data  just  as  well. 

At  Rice  Valley,  the  differences  between  decay  constants  is  probably  not  significant.  With 
the  exception  of  the  2-5  Hz  range,  the  ratio  values  corre.spond  to  two-dimensional  tanh"'lY.,jl 
spreads  that  are  within  the  standard  deviation  of  the  data  (see  Figure  5.8).  The  large  ratio 
values  between  2  and  5  Hz  correspond  to  inversions  where  a^  was  unable  to  converge  to  a 
positive  value.  This  is  obviously  non-phy.sical  since  it  implies  that  signal  correlation  improves 
with  increasing  scn.sor  .separation.  The  decay  constant  at  these  frequencies  is  given  the  value 
of  infinity  in  Table  5.2.  The  cause  of  this  could  be  the  disparity  in  .scn.sor  separation  range. 
Less  constraint  is  placed  on  ai  becau.se  the  maximum  .sensor  separation  in  this  direction  is  four 
times  smaller  than  in  the  aj  direction.  The  F  statistics  in  Table  5.3,  obtained  by  .setting  ai  = 
actually  represent  relatively  good  fits  to  the  data,  however  F  statistics  only  10%  greater  can  be 
obtained  by  setting  ai  =  ap.  indicating  little  .sensitivity  of  3l  to  the  data.  In  any  case,  over  the 
entire  bandwidth  the  i.sotropic  model  fits  the  data  about  equally  well. 
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There  is  also  little  indication  of  directionality  at  Savahia  Mountain.  Below  5  Hz  ilic  ratio 
values  are  very  small  and  inconsistent  in  their  inequality.  At  higher  frequencies,  a[  again 
failed  to  converge  to  positive  values.  As  at  Rice  Valley,  the  quality  of  fit  is  little  affected  by 
changes  in  aL.  And  here  again,  the  fit  of  the  directional  model  is  no  better  than  the  isotropic. 
However,  because  of  the  insufficient  array  geometry  we  cannoi  exclude  the  possibility  of  direc¬ 
tional  correlation  at  these  two  array  sites.  However  at  Ruby  Valley  the  large  differences  in 
decay  constants  and  the  consistency  of  their  sense  of  inequality  with  frequency  indicates  that 
directionality  is  significant. 

5.4  Pg  Coherence 

As  was  the  case  for  Pn,  the  Pg  waveforms  arc  aligned  and  windowed  beginning  at  the 
onset  of  the  phase.  At  Ruby  Valley  the  time  window  is  1.3  seconds  long,  like  Pn.  At  Rice 
Valley  and  Savahia  Mountain  the  Pg  window  length  is  2.6  seconds  long.  A  longer  window 
length  is  used  here  since  it  docs  not  significantly  degrade  the  signal  correlation  ;ind  provides 
greater  frequency  resolution  (Af  =  0.39  Hz).  The  Pg  tanh^’lYjkl  values  used  in  the  inversions 
are  displayed  in  Figures  D4,  D5,  and  D6  of  Appendix  D. 

The  F  statistics  for  Pg  are  given  in  Table  5.4.  At  Ruby  Valley  and  Savahia  Mountain  the 
exponential  models  again  provide  the  best  fits.  This  is  true  also  at  Rice  Valley  beyond  6  Hz, 
however  below  6  Hz  the  self-similar  models  arc  superior.  The  standard  deviations  of  the  best 
model  fits  is  shown  in  Figure  5.13.  The  notable  improvement  in  high-frequency  fit  of  Pg  rela¬ 
tive  to  Pn  at  Ruby  Valley  is  not  significant.  Figure  D4  shows  the  Pg  data  to  be  near  or  below 
the  noise  level  at  high  frequencies  over  the  entire  distance  range.  The  standard  deviations  at 
Rice  and  Savahia  arc  comparable  to  those  for  Pn. 

The  Pg  decay  constants  are  li.stcd  in  Table  5.5.  Looking  back  to  Figure  5.11c,  tlic 
exponential  i.sotropic  decay  constants  arc  plotted  for  each  array  site.  Rice  Valley  consistently 
has  the  greatest  decay  constants.  The  decay  constants  at  Ruby  Valley  and  Savahia  Mountain 
are  very  similar,  with  the  exception  of  the  6-9  Hz  bandwith  where  the  Savahia  values  are 
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notably  larger.  The  a*f  values  in  Figure  5.1  Id  again  show  a  good  deal  of  variation,  though  at 
Ruby  Valley  the  values  are  much  closer,  relative  to  Pn,  to  the  constant  values  predicted  by  the 
linear  a=pX  relationship  (p  a  constant).  In  general  however,  it  does  not  appear  that  incorporat¬ 
ing  this  simple  linear  relationship  into  the  exponential  isotropic  model  would  be  appropriate  for 
our  data  set. 

The  two-dimensional  exponential  decay  constant  ratios  are  plotted  in  Figure  5.14.  At 
Ruby  Valley  aL  is  again  consistently  greater  than  ay,  however  the  ratios  are  about  five  times 
smaller  than  for  Pn.  Though  the  ratios  are  relatively  low,  their  consistency  of  magnitude  and 
sense  of  inequality  suggests  the  directionality  is  significant  below  7  Hz,  where  the  fit  is  also 
slightly  superior  to  the  isotropic  model.  At  Rice  Valley,  the  ratios  are  again  probably  not 
significant,  as  suggested  by  the  very  erratic  ratio  values.  Below  about  8  Hz,  the  three  direc¬ 
tional  models  often  inconsistent  in  their  sense  of  inequality,  unlike  Pg  at  Ruby  Valley  where 
the  directionality  appears  robust  enough  to  produce  consistent  inequalities  across  the  three 
model  types.  The  larger  Rice  Valley  ratios  beyond  10  Hz  are  also  suspect  given  that  so  much 
of  the  data  used  in  the  inversion  falls  below  the  noise  level.  Also,  the  fit  of  isotropic  models, 
the  self-similar  above  6  Hz  and  the  exponential  below  6  Hz,  is  comparable  to  the  directional 
models.  At  Savahia  Mountain  we  have  the  same  problem  as  we  had  for  Pn.  The  directional 
exponential  model  is  again  a  poorly  constrained  model,  with  aL  converging  to  negative  values. 
Restricting  the  inversion  to  station  separations  of  1.5  km  or  less,  thus  omitting  low-amplitude, 
possibly  noise-contaminated  tanh'Myjkl  values,  does  not  improve  the  convergence.  And  again, 
the  isotropic  model  provides  as  good  a  fit  over  the  entire  width. 

Overall  we  cannot  conclude  significant  directional  properties  in  signal  correlation  at  the 
Rice  Valley  and  Savahia  Mountain  array  sites  for  either  Pn  or  Pg.  This  is  not  to  say  that 
directional  characteristics  do  not  exist.  They  might  still  be  revealed  with  more  spatially  uni¬ 
form  distributions  of  stations.  It  does  appear,  however,  that  directional  properties  exist  at  the 
Ruby  Valley  array  site.  Specifically,  at  frequencies  below  10  Hz  for  Pn  and  below  about  7  Hz 
for  Pg,  the  wavefields  dccorrelate  more  rapidly  with  distance  along  the  transverse  dircctioa 


131 


The  directional  sense  of  decorrelation  agrees  with  the  studies  referenced  earlier.  We  will  now 
put  our  correlation  results  to  the  task  of  constructing  models  for  array  gain. 


Modeling  Array  Gain 

The  objective  here  is  to  use  our  spatial  measurements  of  lyjkl  to  construct  continuous 
models  for  the  cross-spectrum,  which  can  then  be  used  in  array  processing  simulations.  We 
will  form  the  cross-spectrum  using 

S(f.r)  =  ly(f,r)l  (5.2) 


where  r  is  the  relative  displacement  vector  between  sensor-pair  locations.  This  relation  is  a 
simplification  of  (5.1)  in  that  it  assumes  unit  amplitude  of  the  cross-spectrum.  The  phase  term 
is  modeled  from  a  propagating  plane  wave  having  wavenumber  vector  k,  and  so  is  given  by 
6(f,r)  =  k  •  r.  The  model  cross-spectra  are  constructed  by  inserting  the  parameterized  coher¬ 
ence  models  into  (5.2).  We  will  use  these  models  to  construct  predictions  of  array  signal  gain. 
Recall  from  Chapter  4  that  array  gain  reflects  the  overall  signal  correlation  at  an  array  and  is  a 
measure  of  how  well  the  detectability  of  a  phase  can  be  improved  through  array  processing. 
From  Equations  (4.3)  and  (5.2)  we  have  the  gain  expressed  as 


=  J_  max  [b(k)  S(0  b(k)]  _  _L  ^  iv  ml 

Em.  ' 

N 


(5.3) 


5J.1  Pn  Array  Gain 

Before  performing  array  processing  simulations  with  the  coherence  models  weTl  first 
examine  how  modeled  array  gain  compares  with  observed  array  gain,  and  use  this  information 
to  make  any  necessary  adjustments  to  the  models.  Beginning  at  Ruby  Valley,  Figure  5.15a 
shows  the  Pn  model  gain  (5.2)  based  on  the  decay  constants  in  Table  5.2.  The  directional 
exponential  model  is  used  up  to  9  Hz.  Beyond  9  Hz,  where  there  is  much  less  suggestion  of 
directionality,  the  exponential  isotropic  model  is  used.  The  observed  array  gain  computed 
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directly  from  the  cross-spectra  of  the  data  is  also  shown;  this  gain  is  equivalent  to  the  array 
gain  of  Chapter  4  (see  Equation  4.3  and  Figure  4.3).  The  two  are  close  in  agreement  up  to 
about  5  Hz.  At  higher  frequencies  the  data  gain  drops  to  much  lower  values  than  the  model 
gain.  Figures  5.15b  and  5.15c  explain  this  behavior.  In  Figure  5.15b  the  model  is  compared 
to  the  magnitude  coherence  gain  of  the  data,  i.e.,  observed  values  of  lyjkCOl  are  used  in  (5.3). 
These  gain  values  are  in  good  agreement  with  the  model  gain  values  over  the  entire  bandwidth. 
In  Figure  5.13c  the  model  values  are  contrasted  with  the  complex  coherence  gain,  i.e.,  we  sim- 

ply  include  the  cross-spectral  phase  term  and  compute  the  gain  using  G  =  — We 

see  that  when  we  add  the  phase  term  the  gain  values  falls  to  much  lower  levels  beyond  5  Hz. 
The  reason  for  this  is  that  at  higher  frequencies  the  phase  becomes  increasingly  more  variable, 
which  has  the  effect  of  scaling  down  the  gain.  At  lower  frequencies  the  cross-spectral  phase  is 
consistently  near  zero  and  so  has  little  effect  on  the  gain.  The  array  signal  gain  in  Figure 
5.15a  is  almost  identical  to  the  complex  coherence  gain,  the  slight  differences  owing  to  the  fact 
that  for  the  coherence  gain  cross-spectral  normalization  takes  place  before  the  matrix  multipli¬ 
cation  with  the  beamsteering  vectors  and  for  the  signal  gain  it  takes  place  after.  The  model 
gain  and  the  observed  signal  gain  in  Figure  5.15a  can  be  brought  into  agreement  with  the 
model  correction  factors  in  Figure  5.15d.  These  factors,  when  multiplied  by  the  model  gain 
values,  produce  the  observed  array  signal  gain.  Below  5  Hz  the  correction  factors  arc  near 
unity.  They  become  more  significant  with  increasing  frequency.  These  model  correction  fac¬ 
tors  will  be  used  to  scale  down  the  model  gains  that  we  simulate  in  the  following  section. 

A  similar  analysis  is  done  for  Pn  at  Rice  Valley  in  Figure  5.16.  Here,  rather  that  exam¬ 
ine  the  entire  4  km  aperture  of  the  array,  the  array  is  subdivided  into  all  possible  1.5  aperture 
sub-arrays.  This  allows  us  to  compare  the  results  with  those  from  the  1.5  km  aperture  Ruby 
Valley  array,  and  it  will  give  the  model  correction  factors  less  statistical  uncertainty.  In  all, 
there  arc  32  sub-arrays  with  between  15  and  17  elements  each.  The  observed  signal  gain  from 
each  array  is  plotted  in  Figure  5.16a;  there  is  a  great  deal  of  variation  above  5  Hz.  The  model 
gains  for  each  sub-array  arc  plotted  in  Figure  5.16b.  Because  the  existence  of  directional 
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decorrelation  properties  at  Rice  Valley  is  not  conclusive,  the  exponential  isotropic  model 
parameters  are  used  for  the  model  gains.  Plotted  with  these  gains  at  each  frequency  arc  the 
sub-array  magnitude  coherence  gain  values  which,  like  for  Ruby  Valley  in  Figure  5.15b, 
should  be  close  to  the  model  values.  The  differences  however  are  greater  here.  This  is  due 
primarily  to  the  variation  in  decorrelation  properties  among  the  sub-arrays  and  the  fact  that  the 
model  parameters  represent  an  average  of  these  properties  over  the  entire  4  km  aperture.  The 
sub-array  data  values  do,  however,  fall  within  the  range  of  model  gains.  There  is  a  relative 
offset  of  the  model  gains  at  higher  frequencies.  This  is  due  to  a  noise  effect  mentioned  earlier, 
i.e.,  the  inversion  for  these  small  decay  constants  receives  a  positive  bias  by  the  noise  values  at 
distances  beyond  1.5  km  (see  Figure  D2c  in  Appendix  D).  These  high-frequency  model  gains 
can  be  brought  to  the  mean  of  the  sub-array  values  if  the  inversion  is  restricted  to  station 
separations  of  within  1.5  km.  The  model  correction  factors  are  plotted  in  Figure  5.16c.  There 
is  a  value  at  each  frequency  for  every  sub-array  which  brings  the  model  gains  in  Figure  5.16b 
into  agreement  with  the  observed  gains  in  Figure  5.16a.  The  mean  of  the  correction  factors, 
which  is  also  plotted,  becomes  more  significant  (i.e.,  lower  values)  with  increasing  frequency. 
It  is  this  mean  of  the  model  correction  factors  which  will  be  used  in  the  later  array  gain  simu¬ 
lations.  The  mean  of  the  model  gains  in  Figure  5.16b,  when  scaled  by  the  mean  correction 
factors,  produces  the  mean  sub-array  gain  shown  in  Figure  5.16d.  This  mean  gain  can  of 
course  also  be  computed  from  an  average  of  the  observed  gain  in  Figure  5.16a.  Note  that  our 
overestimation  of  high-frequency  decay  constants  due  to  noise  contamination  is  compensated 
by  lower  model  correction  values. 

The  same  sequence  of  plots  for  Pn  at  Savahia  Mountain  is  shown  in  Figure  5.17.  Here 
there  are  a  total  of  97  1.5  km  aperture  sub-arrays,  having  between  52  and  63  elements  each. 
The  sub-array  gains  are  in  general  less  than  at  Rice  Valley.  Again,  because  directionality  is 
not  conclusive,  the  model  gains  are  ba.scd  on  the  isotropic  exponential  model  parameters  listed 
in  Table  5.3.  The  model  correction  factors  arc  more  significant  here  than  at  Rice  Valley  and 
the  mean  Pn  gain  in  Figure  5.16d  is  consistently  about  0.2  gain  units  less  than  the  Rice  Valley 
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Pn  gain. 

We  should  not  expect  that  the  model  correction  factors  derived  here  for  a  1.5  km  aperture 
to  be  applicable  to  other  array  aperture  dimensions.  For  example,  a  larger  aperture  array  will 
have  station  pairs  with  greater  decorrelation  and  more  erratic  cross-spectral  phases,  resulting  in 
more  significant  model  correction  factors.  3  km  aperture  arrays  are  examined  in  Figures  5.18 
and  5.19.  At  Rice  Valley  there  are  15  sub-arrays  with  this  aperture,  having  33  or  34  elements 
each.  The  variation  in  sub-array  gain  (Figure  5.18a)  is  greatly  reduced  relative  to  the  1.5  km 
gains.  The  mean  gain  levels  are  also  reduced.  Notice  that  there  is  generally  better  agreement 
between  the  model  gains  and  the  sub-array  gains  in  Figure  5.18b.  This  is  because  the  3  km 
aperture  over  which  the  sub-array  gains  were  computed  is  much  closer  to  the  4  km  over  which 
the  model  decay  constants  were  derived.  As  expected,  the  model  correction  factors  have  also 
become  more  significant,  particularly  beyond  8  Hz.  Similar  results  for  the  45  3  km  arrays  at 
Savahia  Mountain  (11!  to  113  elements  each)  are  displayed  for  Pn  in  Figure  5.19.  The  dispar¬ 
ity  of  coherence  gain  values  in  Figure  5.19b  beyond  7  Hz  is  due  to  the  large  number  of 
tanh'Myj^l  values  falling  below  the  noise  level  but  above  the  fit  of  the  model.  Recall  that  at 
these  low  decay  constant  values,  tanh"'lYj,jl  is  bounded  very  nearby  by  zero;  the  distribution  is 
no  longer  normal  and  values  of  tanh  get  "piled  up".  These  low  decay  constant  models 
decay  so  rapidly  with  distance  that  they  necessarily  underestimate  the  noise  gain  over  an  aper¬ 
ture  of  this  size,  (see  Figure  D3  in  Appendix  D). 

5,5.2  Pg  Array  Gain 

An  examination  of  Pg  array  gain  at  the  three  array  sites  paralleling  that  for  Pn  is 
displayed  in  Figures  5.20  to  5.24.  The  Pg  model  gain  at  Ruby  Valley  in  Figure  5.20a  is  com¬ 
puted  using  the  directional  exponential  model  up  to  7  Hz.  At  higher  frequencies,  where  there 
is  little  evidence  for  directionality,  the  isotropic  exponential  model  is  used.  The  difference 
between  the  model  and  data  increa.scs  greatly  beyond  5  Hz,  as  it  did  for  Pn.  The  relatively 
low  values  of  the  model  relative  to  the  coherence  gain  beyond  5  Hz  in  Figure  5.20b  is  again 
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due  to  the  small  decay  constants  and  the  larger  number  of  tanh'Myjjjl  values  falling  below  the 
noise  level  at  these  higher  frequencies  (see  Figure  D4  in  Appendix  D).  The  Pg  correction  fac¬ 
tors  are  a  bit  more  significant  relative  to  Pn.  but  the  behavior  with  frequency  is  very  similar. 

The  Pg  gains  at  the  1.5  km  sub-arrays  at  Rice  Valley  are  shown  in  Figure  5.21.  Like  Pn, 
there  is  a  wide  range  of  signal  gain.  The  model  gains  are  computed  using  the  best-fitting  iso¬ 
tropic  models.  These  are  the  self-similar  model  below  6  Hz  and  the  exponential  model  above.  • 

The  disparity  between  the  model  and  coherence  gains  beyond  6  Hz  can  again  be  remedied  by 
restricting  the  inversions  to  sensor  separations  to  within  1.5  km.  The  model  correction  factors 
are  a  bit  more  significant  than  for  Pn  and  the  mean  sub-array  gain  is  about  0.1  units  less  than 
for  Pn.  At  Savahia  Mountain  in  Figure  5.22  the  sub-array  gains  also  span  a  wide  range  of 
values.  The  mean  Pg  signal  gain  falls  between  about  0.1  and  0.2  gain  units  below  Pn  and 
about  0.3  gain  units  below  the  Pg  gain  at  Rice  Valley.  The  correction  factors  arc  comparable 
to  Pn,  except  between  5  and  10  Hz  where  the  Pg  values  are  up  to  0.2  units  lower.  For  the  3 
km  aperture  sub-arrays  in  Figures  5.23  and  5.24  we  find,  as  we  do  for  Pn,  that  with  respect  to 
the  1.5  km  sub-arrays  there  is  (a)  less  variation  in  signal  gain,  (b)  decreased  mean  signal  gain, 
and  (c)  more  significant  model  correction  factors. 

The  Pn  and  Pg  gain  calculations  above,  performed  as  a  check  on  the  applicability  of  our 
model  parameters,  have  been  invaluable.  We  find  that  due  primarily  to  stochastic  lluctuations 
in  cross-spectral  phase  away  from  that  predicted  by  the  best-filling  planar  wavefront,  our  model 
parameters  consistently  overestimate  the  array  signal  gains.  Correction  factors  had  to  be 
derived  to  bring  the  models  into  agreement  with  the  data.  We  also  find  extreme  variation  in 
signal  gain  from  the  1.5  km  sub-arrays  even  though  they  arc  themselves  confined  within  a  4  * 

km  aperture.  This  again  underscores  the  stochastic  nature  of  these  wavcficlds  and  the  need  to 
approach  them  in  a  statistical  fashion,  as  we  have  done  here  by  matching  average  gains  for  a 
given  aperture.  Now,  with  our  revised  models  in  hand,  we  arc  ready  to  proceed  to  simulations 
of  arrays  gain. 
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5^.3  Array  Gain  Simulations 


As  discussed  at  the  outset,  a  principal  reason  for  parameterizing  the  coherence  was  to  use 
the  results  to  construct  spatially  continuous  models  of  the  cross-spectral  matrix  (see  Equation 
5.?).  This  allow*:  js  tc  sim’ ’ate  array  procc.'sirig  characteristics  for  an  arbitrary  array 
configuration  and  wavefield  propagation  direction.  Here  we  will  use  this  fact  to  reference  our 
array  gains,  which  are  based  on  three  very  different  array  configurations,  to  a  single  array 
geometry.  This  will  remove  differences  in  gain  due  to  differences  in  the  arrangement  of  sen¬ 
sors  and,  if  significant,  azimuthally-dependent  differences  arising  from  directional  correlation 
properties.  The  array  design  we  will  use  is  that  used  for  the  NORESS  array.  The 
configuration  is  shown  in  Figure  5.25  and  is  very  similar  to  ARCESS,  FINESA,  and  CRESS, 
the  other  recently  deployed  regional  arrays  mentioned  in  Chapter  1.  The  NORESS  stations  are 
laid  out  in  four  concentric  rings  A,  B,  C,  and  D,  with  one  station  Aq  located  at  the  center. 
We’ll  use  two  geometries:  (1)  rings  A,  B,  and  C  with  Aq  (1.5  km  aperture),  and  (2)  rings  A, 
B,  C,  and  D  with  Aq  (3.0  km  aperture).  The  basic  difference  between  the  two  arrays  is  that 
the  larger  has  greater  resolving  power  in  frequency-wavenumber  space.  This  advantage  can  be 
negated,  however,  if  the  signal  looses  significant  coherence  over  the  larger  intersensor  separa¬ 
tions.  The  Pn  and  Pg  array  gains  for  these  two  array  geometries  are  shown  in  Figure  5.26. 
The  model  parameters  used  are  the  same  as  those  used  above  in  Figures  5.15  to  5.24.  The 
gains  have  been  scaled  by  the  model  correction  factors  appropriate  for  the  indicated  aperture. 
The  mean  correction  factors  are  used  for  the  Rice  Valley  and  Savahia  Mountain  gains. 

For  the  1.5  km  array  the  gain  is  typically  the  greatest  at  Rice  Valley;  Ruby  Valley  is  only 
slightly  greater  below  4  Hz.  The  lowest  Pn  gain  is  consistently  at  Savahia  Mountain.  Rice 
Valley  also  maintains  the  greatest  Pg  gain,  by  an  even  greater  margin  than  for  Pn.  The  Ruby 
Valley  Pg  gain  is  notably  greater  than  Savahia  below  5  Hz,  but  the  two  become  indistinguish- 
ably  low  at  higher  frequencies.  At  all  array  sites  the  Pn  gain  is  greater  than  the  Pg  gain. 
Analogous  3  km  gains  are  shown  in  Figure  5.26c  and  5.26d.  The  model  correction  factors 
used  for  Ruby  Valley  correspond  to  its  1.5  km  aperture,  and  .so  probably  overestimate  the  3 
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km  gain  simulations  -  the  correction  factors  for  Rice  Valley  and  Savahia  Mountain  are  typi¬ 
cally  more  significant  for  the  3  km  sub-arrays.  A  better  representation  might  be  to  reduce  the 

1.5  km  correction  factors  by  the  average  of  the  amount  that  these  factors  were  reduced  at  the 
two  other  array  sites.  The  3  km  gain  is  in  all  cases  lower  the  the  1 .5  km  gain.  The  Pg  gain 
has  fallen  extremely  low,  with  the  most  dramatic  drop  taking  place  at  Rice  Valley. 

As  a  final  example,  we  will  use  the  fact  that  we  can  simulate  the  array  gain  from  an 
arrival  having  an  arbitrary  azimuth  to  examine  the  effects  of  directional  decorrelation  on  array 
gain.  Figure  5.27  shows  Pn  gain  based  on  Ruby  Valley  exponential  decay  constants  and  the 

1.5  km  NORESS  configuration.  Three  gain  curves  are  plotted.  One  is  based  on  the  isotropic 
decay  constants  and  so  is  independent  of  wavefront  propagation  direction.  The  others  are 
based  on  the  directional  decay  constants  up  to  9  Hz,  (and  again  isotropic  beyond)  with  two 
orthogonal  directions  of  propagation,  one  propagating  due  north  (0  degrees)  and  the  other  pro¬ 
pagating  due  west  (90  degrees).  We  see  that  even  though  the  directional  decay  constants  differ 
by  up  to  a  factor  of  18,  the  orthogonal  gains  arc  practically  identical  and  only  slightly  lower 
than  the  isotropic  gain  below  9  Hz.  This  similarity  is  a  positive  attribute  of  the  NORESS  array 
configuration,  which  samples  longitudinal  and  transverse  station  separations  about  equivalently, 
independent  of  propagation  direction.  Were  the  configuration  a  linear  one  similar  to  those  of 
Rice  Valley  and  Savahia  Mountain,  we  would  see  a  more  significant  azimuthal  dependence 
with  these  values  of  directional  decay  constants.  Array  simulations  of  this  kind  can  be  an  aid 
to  effective  array  design  once  one  has  established  signal  correlation  models  for  a  potential 
array  site  of  intere.st. 

5.6  Summary  Discussion 

The  aim  of  this  chapter  was  to  contrast  the  signal  correlation  properties  at  our  three 
regional  array  sites  and  in  doing  so  provide  a  means  of  constructing  signal  correlation  models 
that  could  be  used  to  simulate  array  processing  characteristics.  Below  we  summarize  and  dis¬ 
cuss  the  principal  results. 
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Modeling  Procedure 

Several  notable  aspects  of  the  modeling  procedure  are  worth  summarizing.  One  is  the 
use  of  multiple-taper  spectral  estimation  techniques,  which  improved  the  frequency  resolution 
of  the  estimates  bv  decreasing  the  amount  of  frequency  avenging  needed  to  obtain  acceptab’e 
statistics  (see  Appendix  C).  Another  is  the  tanh'^lyl  transformation,  which  produced  a  more 
nearly  normal  distribution  of  data,  making  a  least  squares  inversion  more  appropriate.  In  addi¬ 
tion,  because  we  have  no  way  of  knowing  a  priori  how  the  spatial  coherence  behaves,  we 
were  careful  to  examine  a  range  of  correlation  models.  No  doubt,  even  somewhat  better 
parameterizations  than  those  used  here  can  be  found,  however  the  ability  of  the  models  to 
match  the  magnitude  coherence  gain  indicates  that  the  model  fits  are  adequate.  One  obvious 
refinement  to  the  modeling  procedure  is  to  restrict  the  inversion  to  sensor  separations  to  within 
the  array  aperture  being  simulated.  This  would  further  reduce  model  misfit  and  avoid  much  of 
the  noise  bias  that  comes  at  high  frequencies  and  large  sensor  separations.  Although  the 
configurations  of  the  Rice  Valley  and  Savahia  Mountain  array  geometries  precluded  an  under¬ 
standing  of  directional  correlation  properties,  we  were  nevertheless  very  fortunate  to  have  such 
a  high  density  of  stations.  The  large  number  of  recordings  at  these  sites  allowed  us  to  obtain 
very  statistically  significant  estimates  of  isotropic  decay  constants  and  mean  array  gain. 

Modeling  Results 

The  exponential  models  typically  provided  the  best  fits  to  the  tanh'Myji^l  data.  The  self¬ 
similar  model  gave  thf  best  fit  in  a  few  cases  and  the  gaussian  model  usually  gave  the  worst  fit 
of  all.  We  can  argue  for  the  directionality  of  decay  constants  only  at  Ruby  Valley,  where  we 
find  for  Pn,  and  to  a  lesser  extent  Pg,  that  the  wavefield  decorrelates  more  severely  transverse 
to  the  direction  of  propagation.  This  larger  transverse  decorrelation  reflects  the  significant 
spread  in  propagation  directions  that  we  found  in  Chapter  4  for  Pn  and  Pg  at  Ruby  Valley  (.see 
Figures  4.5  and  4.16).  Again,  we  cannot  exclude  the  possibly  that  directionality  aKso  exits  at 
the  Rice  Valley  and  Savahia  Mountain  sites.  The  greatest  correlation  of  the  Pn  and  Pg 
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wavefields,  measured  in  terms  of  the  magnitude  of  the  isotropic  decay  constants,  exits  at  the 
Rice  Valley  site.  This  could  be  due  to  the  combined  effects  of  less  significant  subsurface 
heterogeneity  and  topographic  scattering  effects. 

Simulations 

The  array  gain  simulations  were  performed  so  that  geometry-  and  azimuth-independent 
comparisons  could  be  made  betwc»:n  the  array  sites.  They  also  serve  as  an  example  of  the 
kinds  of  simulations  which  can  be  performed  as  an  aid  to  the  design  of  an  array  at  a  potential 
recording  site,  or  to  enhance  understanding  of  wavefield  properties  at  an  existing  array  site. 
However,  the  mechanics  of  array  gain  simulations  were  not  as  straight  forward  as  initially- 
thought.  The  .stochastic  behavior  of  the  cross-spectral  phase  at  Ruby  Valley  required  that 
adjustments  be  made  to  the  parameterized  models,  our  model  correction  factors.  Then,  even 
more  randomness  was  revealed  in  the  wide  range  of  sub-array  gains  at  Rice  Valley  and 
Savahia  Mountain.  Mean  correction  factors  were  computed  which  were  themselves  found  to 
depend  on  array  aperture.  These  stochastic  manifestations  of  Pn  and  Pg  underscore  the  fact 
that  study  of  these  wavefields  requires  good  statistical  sampling. 

The  analysis  of  simulated  array  gains  showed  that  effective  array  processing  of  Pn  much 
beyond  a  few  Hz,  and  of  Pg  perhaps  at  any  frequency  much  above  1  Hz,  may  not  be  possible 
for  an  array  dimension  as  large  as  3  km.  We  also  find  that  though  differences  in  directional 
decay  constants  will  result  in  differences  in  the  ability  to  resolve  phase  velocity  and  propaga¬ 
tion  azimuth,  they  do  not  result  in  a  strong  azimuth-dependency  of  array  gain  as  long  as  there 
is  a  fairly  uniform  two-dimensional  distribution  of  recording  stations. 

We  have  obviously  not  exhausted  all  possible  array  processing  simulations.  Any  compu¬ 
tations  involving  the  cross-spectral  matrix  can  be  explored.  This  includes  frequency- 
wavenumber  power  spectra,  as  in  Figures  5.1  and  5.2,  and  other  measures  of  signal  enhance¬ 
ment  which  incorporate  the  correlation  structure  of  noise.  Addi'ionally,  the  ability  to  simulate 
correlation  characteristics  can  facilitate  the  determination  of  complex  station  weighting 
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parameters  aimed  at  optimizing  array  processing  capability. 
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Coherence  Models 

lyi  model 

isotropic 

directional 

Exponential 

e-r/a 

g-x/aT  g-y/a^ 

Gaussian 

e-r^/a" 

Self-similar 

tanh(KQ(r/a)) 

tanh(Ko(x/ax))  tanhCKoCy/aJ) 

Table  5.1  Coherence  models  considered  in  this  study.  For  the  directional  models  ai.  is  the  decay  con¬ 
stant  in  the  longitudinal  direction  (i.e.,  in  the  direction  of  propagation)  and  Sj  is  the  decay  constant  in 
the  transverse  direction  (i.e.,  perpendicular  to  the  direction  of  propagation).  For  statistical  reasons  the 
inversion  for  decay  constants  takes  place  after  performing  inverse  hyperbolic  tangent  transformations  on 
the  coherence  data  and  the  models. 
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Isotropic  Cobereoce  Models 


Isotropic  Models:  Spectrsl  Power 


Isotropic  Models:  Spectral  Power 


FIG.  5.1  (a)  Examples  of  ihe  three  isotropic  coherence  models  used,  (b)  and  (c)  The  corresponding 
spectral  power  in  linear  and  log  plots.  The  decay  constant  here  is  1.0  for  all  three  models.  The  power 
for  each  model  is  normalized  to  a  value  of  one.  (El  =  Exponential.  G1  =  Gaussian,  and  SI  =  Self 
Similar) 
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—  2  3.'j'  —  1 


FIG.  5.2c  Same  as  Figure  5.2b  with  decay  constants  differing  by  a  factor  of  two,  aL  =  2aj. 
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Savahia  Mountain 


EX:  a  =  0.14 
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FIG.  5.3  Random  noise  tanh'Myl  values  for  all  station  pairs  at  each  array  site.  The  lanh^'lyl  90% 
significance  level  for  noise  is  approximately  0.7.  The  best-fitting  isotropic  model  is  also  plotted;  the 
decay  constant  and  model  type  are  indicated  (EX  =  Exponential). 
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FIG.  5.4  (a)  Distribution  of  the  longitudinal  and  transverse  components  of  inter-sensor  separations  for 
all  66  station  pairs  at  the  Ruby  Valley  array,  (b)  Histogram  of  absolute  intersensor  separations, 
summed  in  0.2  km  intervals,  (c)  Histogram  of  the  longitudinal-component  separations,  (d)  Histogram 
of  the  transverse-component  separations. 
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fig.  5.5  Same  as  FIG.  5.7  but  for  the  47  stations 
SALUT. 
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FIG.  5.6  Same  as  for  FIG.  5.7  but  for  ihe  155  stations  at  the  Savahia  Mountain  array  that  recorded  the 
explosion  SALUT. 
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tanh  Myl  tanh  'lyl 


2-D  Exponential  Model  at  Ruby  Valley 


QL/oT  =  2  al/al  =  5 


oL/oT  =  10 


qL/qT  =  15 


FIG.  5.7  Predicted  tanh  '  lyl  values  at  Ruby  Valley  based  on  the  directional  exponential  model  with  the 
indicated  values  of  decay  constant  ratios.  The  best-fitting  isotropic  exponential  model  is  also  plotted. 
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tanh  Myl  tarih  'I7I 


2-D  Exponential  Model  at  Rice  Valley 

oT/aL  =  2  qT/qL  =  5 


oT/aL  —  10  aJ/aL  =  15 


FIG.  5.8  Same  as  FIG.  S.Tbul  for  Rice  Valley. 
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tanh'MYl  tanh  Myl 


2-D  Exponential  Model  at  Savahia  Mountain 


aL/aT  =  2  oL/aT  =  5 


aL/aT  =  10 


oL/aT  =  15 


FIG.  5.9  Same  as  FIG.  5.7  but  for  Savahia  Mountain. 
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statistics  represent  poorer  model  fits.  (El  =  1-D  Exponential,  E2  =  2-D  Exponential,  SI  =  1-D  Self 
Similar,  S2  =  2-D  Self  Similar,  G1  =  1-D  Gaussian,  G2  =  2-D  Gaussian) 


Standard  Deviations  of  Best  Fitting  Models 


Pn  at  Ruby  Valley 


FIG  5.10  Standard  deviations  of  the  best  Pn  tanh  'lyl  model  6ts  at  each  of  the  array  sites. 
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Decay  Constants  for  Pn  at  Ruby  Valley 


Hz 

SI 

S2 

ax 

G1 

^1-  G2  SLj 

1.6 

59.4 

303.7 

38.8 

2.3 

13.6 

6.1 

6.8 

10.8 

4.9 

2.3 

35.1 

152.0 

23.8 

6.4 

9.8 

4.8 

5.2 

7.8 

3.9 

3.1 

13.5 

139.9 

7.9 

3.9 

8.4 

2.5 

3.2 

6.5 

2.1 

3.9 

68,8 

4.1 

2.8 

6.2 

1.8 

2.3 

4.8 

1.5 

4.7 

3.9 

32.7 

2.4 

4.4 

1.3 

1.7 

3,4 

1.2 

5.5 

2.7 

23.4 

1.6 

1.6 

3.7 

1.0 

1.4 

2.8 

0.9 

6.3 

17.2 

1.2 

1.3 

3.1 

0.8 

1.2 

2.3 

0.8 

7.0 

2.0 

14.1 

1.2 

1.3 

2.7 

0.9 

1.2 

2.1 

0.8 

7.8 

2.0 

11.5 

1.3 

1.3 

2.5 

0.9 

1.2 

2.0 

0.9 

8.6 

2.2 

8.2 

1.6 

1.4 

2.3 

1.0 

1.3 

1.8 

1.0 

9.4 

2.8 

5.3 

2.7 

1.6 

2.0 

1.4 

1.5 

1.7 

1.3 

3.4 

5.3 

2.7 

1.8 

1.7 

2.0 

1.6 

1.7 

1.3 

11.0 

2.5 

2.7 

4.0 

1.5 

1.3 

1.8 

1.4 

1.2 

1,5 

11.7 

1.8 

1.8 

2.8 

1.2 

1.1 

1.4 

1.1 

1.0 

1.3 

1.4 

1.4 

2.4 

0.9 

1.3 

1.0 

0.8 

1.1 

13.3 

1.2 

1.1 

2.1 

0.8 

1.2 

0.9 

0.7 

1.1 

14.1 

1.1 

1.3 

1.6 

0.9 

1.0 

0.9 

0.8 

0.9 

14.9 

1.2 

1.6 

1.3 

1.0 

0.9 

0.9 

0.9 

0.8 

Decay  Constants  for  Pn  at  Rice  Valley 


Hz 

El 

!  at  E2  at 

at  G2 

ai 

1.6 

31.4 

31.0 

44.3 

7.2 

3.3 

11.5 

5.9 

2.7 

9,3 

2.3 

19.2 

oo 

12.1 

5.5 

4.6 

5.8 

4.6 

3.2 

5,1 

3.1 

44.0 

oo 

29.6 

8.6 

7.3 

9.0 

7.0 

5.5 

7.5 

3.9 

45.8 

oo 

28.5 

8.7 

16.4 

8.2 

■a 

8.9 

6,9 

4.7 

18.4 

oo 

16.3 

5.3 

3.3 

6.6 

■a 

2.6 

5.6 

5.5 

6.2 

5.5 

9.3 

2.7 

1.2 

5.2 

2.5 

1.1 

4.4 

6.3 

5.1 

4.2 

8.1 

2.4 

1.1 

4.9 

2.2 

0.9 

4.1 

7.0 

4.1 

3.2 

6.6 

2.0 

0.9 

4.4 

1.9 

0.8 

3.8 

7.8 

3.0 

2.4 

4.9 

1.6 

0.7 

3.7 

1.5 

0.6 

3.3 

8.6 

2.0 

1.6 

3.3 

1.1 

0.5 

2.7 

1.0 

0.5 

2.7 

9.4 

1.8 

1.4 

3.1 

1.0 

0.5 

2.4 

0.9 

0.4 

2.3 

10.2 

1.7 

1.3 

3,1 

1.0 

0.5 

2.4 

0.8 

0.4 

2.3 

11.0 

1.8 

1.2 

3.6 

1.0 

1  0.5 

2.8 

0.8 

0.4 

2.6 

11.7 

2.0 

1.3 

4.0 

1.0 

0.5 

3.1 

0.9 

0.4 

2.9 

12.5 

2.1 

1.5 

3.7 

1.1 

1  0.5 

3.0 

0.9 

0.4 

2.7 

13.3 

2.3 

1.8 

3.8 

1.1 

0.5 

3.1 

1,0 

0.5 

2.7 

14.1 

2.3 

2.1 

3.5 

1.1 

0.5 

2.9 

1.0 

0.4 

2,6 

14.9 

2.4 

1.9 

4.0 

1.1 

0.5 

3.1 

0.9 

0.5 

2.7 

Table  Pn  decay  consianis  in  kilometers  for  each  of  the  six  models  at  each  array  site. 
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Decay  Constants  for  Pn  at  Savahia  Mountain 


Hz 

\~El 

ai.  E2 

.SI 

a£  S2 

ar 

G1 

aL  G2 

ar 

1.6 

5.4 

14.7 

5.7 

2.1 

44.4 

1.6 

2.0 

1.9 

2.1 

2.3 

5.0 

5.1 

8.5 

2.0 

4.6 

1.8 

1.9 

1.5 

2.3 

3.1 

5.5 

4.8 

2.1 

2.0 

2.4 

2.0 

1.4 

2.6 

3.9 

4.5 

4.8 

7.4 

1.8 

176.9 

1.5 

1.8 

1.3 

2.1 

4.7 

4.1 

7.4 

4.9 

1.7 

7.5 

1.3 

1.6 

1.2 

1.7 

5.5 

3.1 

DO 

2.6 

1.3 

6.4 

1.1 

1.2 

1.3 

1.2 

6.3 

2.4 

oo 

2.3 

1.1 

4.3 

0.9 

1.0 

2.3 

0.8 

7.0 

1.8 

OO 

0.8 

11.0 

0.6 

0.8 

3.3 

0.6 

7.8 

1.4 

OO 

0.8 

0.7 

3.8 

0.5 

0.6 

2.2 

0.5 

8.6 

1.3 

CO 

1.1 

0.6 

5.4 

0.4 

0.5 

1.6 

0.4 

9.4 

1.2 

OO 

0.5 

53.7 

0.4 

0.4 

2.4 

0.4 

1.0 

OO 

0.4 

2.7 

0.3 

0.4 

1.4 

0.3 

11.0 

0.9 

OO 

0.4 

4.5 

0.3 

0.3 

1.1 

0.3 

11.7 

14.9 

0.4 

7.8 

0.3 

0.3 

1.2 

0.3 

12.5 

3.2 

IBM 

0.4 

2.1 

0.3 

0.3 

1.0 

0.3 

13.3 

1.0 

1.9 

1.1 

0.4 

18.3 

0.3 

0.3 

0.7 

0.3 

14.1 

1.21 

1.7 

1.6 

0.42 

1.8 

0.3 

0.4 

14.9 

1.1 

1.9 

1.4 

0.4 

2.4 

0.3 

0.3 

IB 

Pn:  El  Decay  Constants 


Pn;  a  *  f 


^  Ruby  Valley 
•  Rice  Valley 
o  Savahia  Mountain 


Pg:  El  Decay  Constonts 


c 

10  1 

d 

/T 

10’  3 

10®- 

Frequency  (Hz) 


Frequency  (Hz) 


FIG.  5.11  (a)  Isotropic  exponential  decay  constants  for  Pn  at  each  array  site.  See  also  Table  5.3.  (b) 
Product  of  the  1-D  exponential  decay  constants  and  frequency.  A  constant  value  of  a*i  over  frequency 
would  indicate  that  the  decay  constant  is  a  linear  function  of  wavelength,  (c)  and  (d)  Same  as  (a)  and 
(b)  but  for  Pg. 


Frequency  (Hz) 


FIG.  5.12  Two-dimensional  exponential  decay  constant  ratios  for  Pn  at  each  array  site.  The  cross- 
hatched  bars  represent  values  of  at/at  and  the  blank  bars  values  of  Except  for  the  11.7  Hz  value 
at  Savahia,  ratios  f’’’  Rice  Valley  and  Savahia  that  reach  a  value  of  10  actually  extend  to  infinity  and 
are  insignificant  (see  text). 
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Standard  Deviations  of  Best  Fitting  Models 

Pg  at  Ruby  Valley 


0.0  5.0  10.0  15.0 


Pg  at  Savahia  Mountain 


Frequency  (Hz) 


FIG  5.13  Same  as  Figure  5.10  but  for  Pg. 
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Decay  Constants  for  Pg  at  Ruby  Valley 


El 

aL  E2  ar 

SI 

at  S2  aT 

G1 

2.4 

5.3 

2.0 

■El 

1.9 

1.2 

1.3 

2.1 

5.2 

1.7 

1.9 

1.1 

1.2 

1.8 

4.6 

1.4 

■a 

1.7 

1.0 

1.1 

1.5 

3.7 

1.2 

1.1 

1.5 

0.9 

1.0 

1.3 

3.3 

1.0 

1.0 

1.4 

0.8 

0.9 

1.0 

2.6 

0.8 

0.8 

1.2 

0.7 

0.8 

1.0 

2.6 

0.8 

0.8 

1.2 

0.6 

0.8 

0.8 

1.7 

0.7 

0.6 

0.9 

0.5 

0.7 

0.7 

1.2 

0.7 

0.6 

0.7 

0.5 

0.6 

0.7 

1.1 

0.7 

0.5 

0.7 

0.5 

0.5 

0.6 

0.9 

0.8 

0.5 

0.6 

0.6 

0.5 

0.6 

0.8 

0.8 

0.5 

0.6 

0.6 

0.5 

0.7 

0.9 

0.8 

0.5 

0.6 

0.6 

0.6 

0.7 

1.0 

0.7 

0.6 

0.6 

0.6 

0.6 

0.6 

0.8 

0.7 

0.5 

0.5 

0.5 

0.5 

0.6 

0.8 

0.7 

0.5 

0.5 

0.5 

0.5 

0.5 

0.8 

0.6 

0.5 

0.5 

0.5 

0.5 

0.5 

0.8 

0.6 

0.5 

0.5 

0.5 

0.5 

G2  aj 


1.6  1.1 

1.6  1.0 

1.5  0.9 

1.3  0.8 

1.2  0.7 

1.1  0.6 

1.1  0.6 

0.8  0.5 

0.7  0.5 

0.6  0.5 

0.5  0.5 

0.5  0.6 

0.5  0.6 

0.6  0.5 

0.5  0.5 

0.5  0.5 

0.5  0.4 

0.5  0.4 


Decay  Constants  for  Pg  at  Rice  Valley 


at  G2  aj 

3.4 

3.5 

2.2 

2.2 

1.3 

3.4 

1.4 

3.2 

1.0 

2.9 

0.9 

3.3 

0.9 

4.4 

1.1 

3.3 

0.5 

1.4 

0.6 

0.7 

Hz  I  El 


at 

S2  at 

9.6 

3.8 

3.4 

2.4 

8.9 

2.3 

1.6 

3.3 

1.2 

3.1 

1.1 

3.5 

1.1 

4.9 

1.4 

3.1 

0.6 

1.7 

0.6 

0.7 

0.6 

0.9 

0.4 

2.3 

0.3 

2.1 

0.3 

1.4 

0.3 

1.7 

0.3 

3.0 

0.2 

2.2 

0.2 

2.7 

10.0 

10.0 

3.5 

2.2 

2.4 

2.3 

1.9 

2.0 

2.2 

2.1 

0.9 


10.0  10.0 


Table  5.5  Same  as  Table  5.3  but  for  Pg. 
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Decay  Constants  for  Pg  at  Savahia  Mountain 


Hz 

El 

E2  Of 

SI 

ai,  S2  ax 

G1 

aL  G2 

a? 

2.6 

oo 

2.5 

1.4 

8.1 

l.l 

1.3 

2.1 

1.1 

2.2 

oo 

2.0 

1.1 

1.3 

1.0 

1.1 

1.3 

1.0 

2.1 

5.5 

2.2 

1.0 

6.5 

0.7 

0.9 

0.9 

0.9 

3.9 

2.1 

6.2 

2.2 

0.9 

4.5 

0.7 

0.8 

0.9 

0.8 

4.7 

1.5 

8.1 

1.5 

0.7 

5.8 

0.5 

0.6 

0.8 

0.5 

5.5 

1.2 

OO 

1.0 

0.5 

2.7 

0.4 

0.5 

1.3 

0.4 

6.3 

1.6 

oo 

1.4 

0.6 

3.3 

0.5 

0.6 

2.0 

0.5 

7.0 

1.6 

oo 

1.3 

0.6 

3.4 

0.4 

0.5 

1.8 

0.4 

7.8 

1.3 

oo 

1.2 

0.4 

6.6 

0.3 

0.3 

3.8 

0.3 

8.6 

1.0 

4.5 

1.0 

0.3 

3.8 

0.2 

0.3 

0.5 

0.2 

9.4 

0.6 

0.5 

1.3 

0.2 

0.1 

5.7 

0.2 

0.1 

0.3 

10.2 

0.7 

0.7 

1.1 

0.2 

0.1 

1.0 

0.2 

0.1 

0.3 

11.0 

0.8 

1.5 

1.0 

0.2 

1.2 

0.2 

0.2 

0.2 

0.2 

11.7 

0.8 

1.7 

0.9 

0.2 

0.3 

0.2 

0.2 

0.3 

0.2 

12.5 

0.8 

oo 

0.6 

0.2 

2.3 

0.2 

0.2 

0.4 

0.2 

13.3 

0.6 

1.3 

0.7 

0.2 

0.2 

0.1 

0.2 

0.2 

0.1 

14.1 

0.5 

0.7 

0.7 

0.1 

1.3 

0.1 

0.1 

0.1 

0.1 

14.9 

0.4 

0.4 

0.5 

0.1 

0.3 

0.1 

0.1 

0.1 

0.1 

Table  5.5  continued 
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Exponential  2-D  Decay  Constant  Ratios 


Pg  at  Ruby  Valley 


Pg  at  Savahia  Mountain 


Frequency  (Hz) 

FIG.  5.14  Two-dimensional  exponential  decay  constant  ratios  for  Pg  at  each  array  site.  Except  for  the 
1.6  Hz  value  at  Rice  valley,  ratios  for  Rice  Valley  and  Savahia  that  reach  a  value  of  10  actually  extend 
to  infinity  and  are  insignificant  (see  text). 
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Pn:  Ruby  Valley 


Array  Signal  Gain  Coherence  Gain 


Coherence  Gain  with  Phase  Model  Correction  Factors 


FIG.  5.15  (a)  Pn  model  array  gain  ai  Ruby  Valley  (solid  line)  and  Pn  array  gain  (solid  dots)  obtained 
from  the  cross-spectra  of  the  data  S,j(f).  (b)  Model  array  gain  in  (a)  and  the  observed  magnitude  coher¬ 
ence  gam  obtained  by  setting  the  cross-spectra  equal  to  the  observed  magnitude  coherence, 
S,j(0  =  (c)  Model  array  gain  in  (a)  and  the  observed  complex  coherence  gain  obtained  by  set- 

ung  the  cross-spectra  equal  to  the  complex  coherence,  S,j(0  =  Y,j(0-  (d)  Correction  factors  which  when 
multiplied  by  the  model  gain  produce  the  observed  array  gain  in  (a). 
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Pn:  Rice  Valley  - 1.5  km  Aperture 


Sub-array  Signal  Gains  Coherence  Gains 


Frequency  (Hz)  Frequency  (Hz) 


Model  Correction  Factors  Mean  Sub-array  Signal  Gain 


Frequency  (Hz)  Frequency  (Hz) 


FIG.  5.16  (a)  Observed  Pn  array  signal  gain  from  the  32  1.5  km  aperture  sub-arrays  at  Rice  Valley, 
(b)  Model  array  gains  for  each  sub-array  configuration.  Also  shown  are  the  observed  sub-array  magni¬ 
tude  coherence  gain  values  obtained  by  seuing  Si/O  =  lYijfOl-  (c)  Model  correction  factors  for  each 
sub-array  which  produce  the  observed  gains  in  (a).  The  mean  correction  factors  fall  along  the  solid 
line,  (d)  Mean  of  the  sub-array  signal  gains  in  (a).  This  is  equivalent  to  the  mean  of  the  model  gains 
in  (b)  scaled  by  the  mean  correction  factors  in  (c). 
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Pn:  Savahia  Mountain  -  1.5  km  Aperture 


Sub-array  Signal  Gains  Coherence  Gains 


Frequency  (Hz)  Frequency  (Hz) 


Model  Correction  Factors  Mean  Sub-array  Signal  Gain 


Frequency  (Hz)  Frequency  (Hz) 


FIG.  5.17  Same  as  Figure  5.16  but  for  the  97  1.5  km  aperture  sub-arrays  at  Savahia  Mountain. 
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Pn:  Rice  Valley  -  3.0  km  Aperture 


Sub— array  Signal  Gains  Coherence  Gains 


Frequency  (Hz)  Frequency  (Hz) 


fvlodel  Correction  Factors  Mean  Sub-array  Signal  Gain 

1.0 

0.8 

0.6 

0.4 

0.2 

0.0 

0  5  10  15  0  5  10  15 

Frequency  (Hz)  Frequency  (Hz) 


FIG.  5.18  Same  as  Figure  5.16  bul  for  the  15  3.0  km  aperture  sub-arrays  ai  Rice  Valley. 
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Pn:  Savahia  Mountain  -  3.0  km  Aperture 

Sub-array  Signal  Gains  Coherence  Gains 


Frequency  (Hz)  Frequency  (Hz) 


0  5  10  15  0  5  10  15 

Frequency  (Hz)  Frequency  (Hz) 


FIG.  5.19  Same  as  Figure  5.16  but  for  ihc  45  3.0  km  aperture  sub-arrays  at  Savahia  Mountain. 
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Pg:  Ruby  Valley 


Array  Signal  Gain  Coherence  Gain 


0  5  10  15  0  5  10  15 

Frequency  (Hz)  Frequency  (Hz) 


Coherence  Gain  with  Phase  Model  Correction  Factors 


0  5  10  15  0  5  10  15 

Frequency  (Hz)  Frequency  (Hz) 


FIG.  5.20  Same  as  for  Figure  5.15  bul  for  Pg. 
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Pg:  Rice  Valley  -  1.5  km  Aperture 


Sub-array  Signal  Gains 
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FIG.  5.21  Same  as  Figure  5.16  but  for  Pg. 
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Pg:  Savahia  Mountain  -  1.5  km  Aperture 
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HG.  5.22  Same  a.s  Figure  5.17  but  for  Pg, 


Pg:  Rice  Valley  -  3.0  km  Aperture 
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FK;.  5.23  Same  as  Figure  5.18  but  for  Pg. 


Pg:  Savahia  Mountain  -  3.0  km  Aperture 
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FIG.  5.24  Same  as  Figure  5.19  but  for  Pg. 
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Predicted  Array  Gains  -  NORRESS  Configuration 


Pn:  1.5  km  Aperture  Pg:  1.5  km  Aperture 
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Pn:  3.0  km  Aperture  Pg:  3.0  km  Aperture 
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FIG.  5.26  Predicted  Pn  and  Pg  array  gains  ba.sed  on  the  NORESS  configuration  and  the  coherence 
models  used  in  the  proceeding  figures  of  array  gain.  The  mean  model  correction  factors  derived  earlier 
have  been  applied  to  obtain  the  Rice  Valley  and  Savahia  Mountain  results.  The  3.0  km  aperture  gains 
for  Ruby  Valley  are  derived  from  the  1.5  km  aperture  correction  factors.  The  propagation  direction  is 
due  north  in  all  cases. 
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Ruby  Valley  Pn  Model  Gain 


NORESS  1.5  km  Aperture 


Frequency  (Hz) 


FIG.  5.27  F*redicied  Pn  array  gain  based  on  the  NORESS  1.5  km  configuration  and  the  Ruby  Valley 
exponential  models.  Plotted  is  the  gain  for  the  isotropic  model  and  gains  for  the  directional  model  using 
orthogonal  directions  of  wavefront  propagation.  The  directional  model  extends  only  to  9  Hz.  The  direc¬ 
tional  gains  essentially  overlap  and  are  only  slightly  lower  than  the  isotropic  gain. 
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Chapter  6 


Recommendations 

Regional  seismic  arrays  in  Scandanavia  have  proven  valuable  for  monitoring  small- 
magnitude  events  of  the  kind  that  would  be  of  concern  under  the  restrictions  of  a  low-yield 
Threshold  Test  Ban  Treaty.  However,  their  utility  in  a  place  like  the  tectonically  active 
western  United  States  has  been  less  well  established.  This  study  was  undertaken  to  begin  to 
assess  the  regional  array  monitoring  characteristics  and  potential  of  high-frequency  Pn  and  Pg 
wavefields  from  the  Nevada  Test  Site.  The  conclusions  of  the  work  have  been  described  fully 
at  the  end  of  the  preceding  chapters.  This  final  chapter  offers  a  number  of  recommendations 
on  the  deployment  and  processing  of  regional  arrays. 

(1)  Amplitude  Spectra 

The  advantages  of  usings  arrays  over  single-station  recordings  sites  was  described  at 
length  in  Chapter  1,  however  it  is  worth  reiterating  the  point  made  regarding  the  statistical 
improvement  of  spectral  amplitude  e.stimates.  In  Chapter  2  we  saw  very  significant  variation  in 
spectral  amplitude  over  fairly  short  distances.  A  factor  of  10  variation  at  a  given  frequency 
over  the  4  km  aperture  arrays  at  Rice  Valley  and  Savahia  Mountain  was  typical.  At  the  same 
time  we  observed  very  similar  mean  spectra  at  high  frequencies  over  the  much  larger  distances 
separating  the  array  sites  themselves.  These  observations  underscore  the  need  for  statistical 
averaging  of  spectral  estimates  to  improve  the  reliability  of  spectral-based  discrimination  tech¬ 
niques.  The  recommendation  here  is  that  (a)  permanent  arrays  be  used  to  achieve  this,  and  (b) 
if  a  single-station  monitoring  site  is  used  instead  that  a  temporary  array  be  installed  around  it 
to  estimate  the  statististical  properties  of  the  site. 
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(2)  Array  Design 

An  effective  array  design  is  critically  dependent  on  the  signal  and  noise  correlation 
characteristics  of  the  site.  These  characteristics  must  be  determined  before  a  permanent  array 
is  installed.  In  particular,  given  the  wide  variability  of  signal  correlation  that  we  observed  in 
Chapter  5,  this  should  be  carried  out  with  a  high  density  of  temporary  sensors  so  as  to  obtain 
statistically  meaningful  results.  In  geologic  regions  similar  to  those  upon  which  the  arrays  of 
this  study  were  located,  it’s  unlikely  that  significant  advantage  can  be  gained  by  installing  per¬ 
manent  arrays  having  apertures  much  beyond  1.5  km.  Exceeding  this  aperture  at  even  the 
NORESS  array,  where  signal  correlation  is  greater,  results  in  degraded  estimates  of  source 
azimuth  (Bame  et  al.,  1990).  The  gain  simulations  in  Chapter  5  showed  that  while  a  1.5  km 
aperture  may  be  appropriate  for  Pn  and  Pg  at  the  Rice  Valley  array  and  Pn  at  the  Ruby  Valley 
array,  it  is  too  large  for  Pg  at  Ruby  Valley  and,  even  more  so,  for  both  Pn  and  Pg  at  Savahia 
Mountain.  Decreasing  the  aperture  can  of  course  increase  array  gain,  but  then  one  needs  to  be 
concerned  about  increased  noise  correlation  and  degraded  resolving  power  of  the  array.  In 
these  cases  it  may  be  better  to  go  in  search  of  more  suitable  locations.  Of  the  three  sites 
examined  in  this  study.  Rice  Valley,  because  of  its  overall  greater  Pn  and  Pg  signal  correlation 
shows  the  most  promise  as  a  potential  monitoring  site.  An  additional  consideration  are  the 
array  characteristics  of  Lg,  which  have  not  been  examined  here.  Because  of  the  importance  of 
this  phase  in  discrimination  and  yield  estimation,  reconaissance  siting  of  regional  arrays  must 
also  take  its  properties  into  account. 

(3)  Using  Noise 

Other  than  compute  spectral  amplitude  levels  in  Chapter  3,  we  have  done  little  analysis  of 
noise  wavefields.  With  very  small  magnitude  events,  however,  the  noise  field  becomes  very 
important  and,  if  possible,  knowledge  of  its  correlation  structure  should  be  taken  advantage  of. 
Der  et  al.  (1988),  for  example,  prc.scnt  a  general  method  for  using  parameterized  models  of 
signal  and  noise  cross-spectra  to  obtain  stations  weights  which  optimize  signal  to  noi.se  ratios. 
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For  later  arriving  phases  like  Pg  and  Lg  one  can  model  the  noise  from  the  coda  of  the  previous 
arrival. 

(4)  Identification  of  Coherent  Scattering  Sources 

An  obvious  hindrance  to  processing  regional  data  is  its  complexity  due  to  scattering. 
However,  in  some  cases  the  scattering  is  coherent  and  its  source  may  be  understood.  For 
example,  careful  examination  of  the  SALUT  wavefields  recorded  at  the  Savahia  Mountain 
array  site  reveals  coherent  low-velocity  wavefronts  propagating  across  the  array  (see  Figures 
2.6  and  2.11).  The  direction  and  low  velocity  suggests  that  these  are  scattered  surface  waves 
from  Savahia  Mountain.  This  kind,  and  other  less  obvious  types,  of  coherent  scattering  is 
amenable  to  sophisticated  frequency-wavenumber  methods,  such  as  the  multiple  signal  charac¬ 
terization  (or  MUSIC)  method  of  Schmidt  (1986).  The  MUSIC  method  is  superior  to  the  con¬ 
ventional  method  in  its  ability  to  detect  multiple  signals,  and  its  utility  should  be  explored  for 
regional  data.  Velocity  filtering  techniques  can  also  be  used  to  isolate  (or  remove)  coherently 
scattered  wavefields.  Additional  analysis  of  the  Savahia  Mountain  and  Rice  Valley  data  sets 
for  coherent  sources  of  scattering  is  warranted.  Note,  however,  that  without  the  high  density 
of  sensors  at  Savahia,  spatial  aliasing  of  this  low  velocity  scattered  wavefield  would  prevent 
these  analysis  methods  from  being  useful.  This  again  stresses  the  need  for  the  temporary 
placement  of  high-density  arrays  at  potential  monitoring  sites. 

(5)  Wavefield  Simulations 

We  have  seen  throughout  this  study  that  regional  wavefields  are  significantly  affected  by 
very  local  variations  in  geologic  structure.  In  some  cases  the  effects  are  systematic,  such  as 
the  velocity-azimuth  bias  at  Ruby  Valley  and  the  coherent  scattering  at  Savahia  Mountain,  and 
in  some  cases  they  are  random,  as  characterized  by  the  coherence  functions  themselves  at  each 
array  site.  Effects  like  these  will  no  doubt  be  typical  of  array  monitoring  sites  in  the  western 
United  States.  This  clearly  presents  problems  for  routine  array  monitoring  of  regional  events. 
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However,  given  a  sufficiently  detaUed  knowledge  of  the  geologic  structure  surrounding  an 
array  site,  we  can  approach  this  problem  through  numerical  wavefield  simulations.  The  effects 
of  the  structure  on  array  processing  characteristics  can  be  determined  by  propagating  synthetic 
regional  wavefields  through  the  structure  using  finite  difference  methods.  Effects  such  as  the 
biasing  of  phase  velocity  and  source  azimuth  could  be  predicted  and  therefore  compensated. 
Coherent  scattering  sources,  once  identified,  can  be  routinely  removed  through  velocity  filtering 
or  simple  subtraction  in  the  frequency-wavenumber  domain  (Gupta  et  al.,  1990).  Because 
detailed  subsurface  'msging  is  so  costly,  a  first  approach  to  array  site  selection  is  to  consider 
regions  which  have  already  undergone  extensive  crustal  imaging.  The  areas  around  Rice  Valley 
and  Savahia  Mountain,  and  other  regions  explored  by  CALCRUST  are  examples.  The  effects 
of  random  crustal  variations  on  array  processing  can  be  also  simulated  as  well,  given  a  distri¬ 
bution  of  the  randomness.  Array  data  itself  can  be  used  to  estimate  these  distributions  (Flatte' 
and  Wu,  1988;  Wu  and  Hatte',  1990). 

In  general,  effective  routine  use  of  high-frequency  regional  monitoring  arrays  in  the 
western  U.S.  will  require  site-specific  studies  aimed  at  fully  understanding  the  distortional 
effects  caused  by  local  structural  irregularities. 
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Appendix  A 


Spectral  Amplitude  Estimation 


This  purpose  of  this  appendix  is  to  describe  a  few  of  the  details  of  the  multiple  taper 
spectral  estimation  procedure  of  Thomson  (1982),  and  to  examine  how  these  tapers  compare  in 
synthetic  tests  for  which  the  true  amplitude  spectrum  is  known.  Results  from  the  more  com¬ 
mon  method  of  simple  cosine  tapering,  described  immediately  below,  are  also  included.  The 
principal  advantage  of  the  multi-taper  method  is  its  ability  to  form  relatively  low  bias,  low 
variance  spectral  estimates.  A  more  complete  description  of  multi-taper  spectral  estimation  and 
its  statistical  advantages  over  other  spectral  estimation  procedures  can  be  found  in  Thomson 
(1982)  and  Park  et  a/.,  (1987). 

Cosine  Tapering 

Let  Uj(tn);  t^=nT;  n=l,...,N  be  a  time  series  of  sample  length  N  (time  length  T)  and  sam¬ 
pling  interval  x  recorded  at  the  station.  The  discrete  Fourier  transform  (DFT)  is  expressed 
by 

Uj(0  =  TSv(gu/ge"‘^’^\  (A.l) 

n=l 

where  v(tn)  is  the  data  taper  and  f  is  the  frequency.  The  estimate  of  the  amplitude  spectrum  is 
then  given  by  Aj(0  =  luj(0l. 

A  common  choice  of  data  taper  is  the  p%  split-cosine  bell  taper.  Shown  in  Figure  Ala  is 
the  20%  taper,  i.e.,  20%  of  the  time  series  is  modulated  by  the  cosine  function.  In  the  fre¬ 
quency  domain,  the  operation  in  (A.l)  is  equivalent  to  a  convolution  between  the  true  fourier 
transform  of  the  time  series  and  the  transform  of  the  data  taper.  In  effect,  the  spectral  estimate 
is  blurred  and  biased  by  the  effects  of  the  spectral  resolution  and  leakage  characteristics  of  the 
taper.  The  resolution  and  spectral  leakage  characteristics  of  this  taper  can  be  seen  in  its  DFT 
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amplitude  spectrum,  shown  in  Figure  A2a.  The  spectrum  is  plotted  as  a  function  of  frequency 
bin  number  J,  where  J/T  are  the  Rayleigh  frequencies,  which  are  the  sampled  frequencies 
returned  from  a  fast  fouricr  transform.  The  resolving  width,  i.e.,  the  width  of  the  main  lobe, 
for  this  taper  is  approximately  2/T,  and  the  amplitude  of  the  first  sidclobc  is  quite  large,  within 
one  order  of  magnitude  of  the  main  lobe. 

Multiple  Taper  Method 

An  alternative  to  conventional  tapering  is  the  method  of  Thomson  (1982),  which  employs 
the  use  of  multiple  data  tapers,  or  “cigentapers”.  These  tapers  v/(tn;T,W);  /=1,.,.,L  arc 
specifically  designed  to  minimize  spectral  leakage  outside  a  chosen  frequency  band  of  width 
2W.  For  each  time  series  one  computes  a  total  of  L  spectra,  or  ‘‘cigcnspccira”: 

Uj/(f)  =  'cXv/(t„)u/Ue"'’^‘^.  (A.2) 

n=l 

If  the  local  variations  in  tlic  spectrum  arc  not  too  extreme,  the  cigcnspcctra  will  be  nearly 
uncorrclated  and  the  estimate  of  the  amplitude  spectrum  Aj(f)  can  be  formed  from  a  weighted 
sum  of  the  cigcnspcctra  (Thomson,  1982), 

Xlw,(f)  Uj,(Ol" 

Af(0  =  '^^4 - •  (A.3) 

Zlw/(Ol" 

1=1 

The  amount  of  spectral  leakage  associated  with  each  cigentaper  is  reflected  in  its  corresponding 
eigenvalue  X,  where  the  amount  of  fractional  leakage  from  outside  the  2W  bandwidth  is  1  -  A.. 
A  set  of  tapers  having  WT  =  P  belongs  to  the  family  of  “Pit”  tapers.  Values  of  X  for  the 
2k,  3k,  and  4k  tapers  arc  given  in  Table  Al.  Here  we  will  consider  only  the  first  L=2P-1 
lowest  order  tapers.  From  Table  Al,  this  ensures  that  the  fractional  leakage  remains  below,  at 
most,  6%  for  each  taper.  This  amounts  to  using  the  first  3,  5,  and  7  of  the  2k,  3k,  and  4k 
tapers,  respectively. 
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The  first  Lp„  tapers  arc  shown  in  Figure  Al;  tapers  of  higher  order  cross  die  time  axis  a 
greater  number  of  times.  Any  single  taper  will  weight  the  time  series  very  unevenly.  However, 
when  all  tapers  are  combined,  the  time  series  is  much  more  evenly  sampled.  The  DFT  ampli¬ 
tude  spectra  of  each  of  the  tapers  are  shown  in  Figure  A2  and  rcficcl  the  differences  in  resolu¬ 
tion  width  and  leakage  between  the  families  of  tapers.  The  whole-width  of  the  main  lobes  can 

be  seen  to  be  4/T,  6/T,  and  8/T  for  27c,  3k,  4k  tapers,  respectively.  Within  a  family  of  tapers, 
the  higher  order  tapers  have  a  successively  greater  amount  of  spectral  leakage.  In  general,  the 
spectral  leakage  from  most  of  the  Pk  tapers  shown  is  less  than  that  from  the  cosine  taper. 
Between  families  of  tapers,  the  amount  of  spectral  leakage  incrca.scs  with  decreasing  resolution 
widih.  The  number  of  "igcnspcclra  L=2P-1  used  in  the  spectral  estimate  increases  going  from 

2k,  3k,  4k  .  Therefore,  while,  for  example,  the  27c  estimate  will  have  greater  rc.solution,  it  will 
also  have  greater  variance  and  leakage  relative  to  the  3jr  and  4k  estimates. 

The  frequency  dependent  weights  w/(0  arc  computed  through  an  adaptive  iteration  tech¬ 
nique  and  provide  a  means  of  reducing  the  bias  due  to  spectral  leakage  of  the  higher  order 
tapers.  The  details  of  this  weighting  scheme  can  be  found  in  Thomson  (1982)  and  Park  et  al. , 
(1987).  In  regions  where  the  spectrum  is  relatively  constant  w;(0  is  approximately  equal  to 
one.  However,  in  regions  where  spectral  leakage  is  mere  signiriL-aiii,  i.»,.,  wiicrc  the  spcclia  is 
changing  more  rapidly,  the  higher  order  weights  (which  correspond  to  higher  order  tapers  hav¬ 
ing  relatively  greater  spectral  leakage)  arc  reduced.  The  method  therefore  provides  an  objec¬ 
tive  means  of  trading  off  variance,  which  is  increased  by  the  weighting,  and  bias,  which  is 
decreased. 

Synthetic  Tests 

The  characteristics  of  multi-tapering  can  be  further  illuminated  by  analyzing  a  synthetic 
seismic  time  series  for  which  the  true  amplitude  spectrum  is  known.  Here  we  can  use  the  dis¬ 
placement  time  .series  from  a  spherically  symmetric  comprcssional  point  .source  having  a  step 
source  time  function,  as  given  by  Sharpe  (1942'*,  and  shown  in  Figure  A3a  tor  a  chosen  set  of 
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physical  parameters.  The  20%  cosine,  2n,  3k,  and  47t  tapered  DFT  spectra  are  compared  in 
the  three  synthetic  tests  described  below.  A  time  window  of  approximately  0.5  seconds  is 
used,  which  amounts  to  spectral  resolution  widths  of  4  Hz,  8  Hz,  12  Hz,  and  16  Hz  for  the 
20%  cosine,  2jc,  3k,  and  47C  estimates,  respectively.  In  all  tests,  to  further  reduce  spectral  leak¬ 
age,  the  tapered  time  series  were  prewhitened  using  a  sixth  order  auto-regressive  filter  deter¬ 
mined  from  the  tapered  time  series. 

Test  1:  Displacement  pulse  only 

Here  amplitude  spectral  estimates  of  the  time  series  as  shown  in  Figure  A3a  arc  com¬ 
puted.  The  theoretical  amplitude  spectrum  is  shown  at  the  top  of  Figures  A3b  and  A3c.  For 
clarity,  both  log-log  and  log-linear  plots  of  the  spectra  are  shown.  The  true  spectrum  is  quite 
simple,  characterized  by  a  constant  low  frequency  level,  a  comer  frequency  near  4  Hz,  and  a 
constant  high-frequency  log-log  slope  of  -2.  All  of  the  tapering  methods  estimate  the  high  fre¬ 
quency  slope  and  comer  frequency  about  equally  well.  The  20%  cosine  taper,  however,  does  a 
relatively  poor  job  at  low  frequencies,  underestimating  the  true  amplitude  level.  This  is  a 
result  of  its  relatively  heavy  tapering  at  the  ends  of  the  time  series,  effectively  shortening  the 
length  of  the  time  scries  and  therefore  removing  some  of  the  low  frequency  energy.  A  number 
of  the  multiple  tapers,  however,  still  apply  relatively  large  weighting  at  the  ends  of  the  time 
scries  and  so  do  not  lose  so  much  low  frequency  energy  (see  Figure  Al).  Of  the  multi-tapers, 
the  27t  tapers  downweight  the  ends  of  the  time  scries  the  most,  and  so  some  underestimation  of 
the  low  frequency  levels  is  also  evident  in  its  spectral  estimate. 

TTic  eigenspcctra  and  corresponding  computed  weights  for  each  taper  arc  shown  in  Figure 
A4.  The  differences  in  eigenspcctra  reflect  the  differences  in  the  data  tapers  used.  For  exam¬ 
ple,  the  second  order  27c  taper  goes  through  a  node  near  where  the  amplitude  of  the  time  scries 
is  greatest,  resulting  in  the  relatively  low  spectral  amplitude  esii.naic  in  Figure  A4a  The 
weighting  for  any  taper  basically  depends  on  the  amount  of  leakage  associated  with  it  and  the 
variation  in  the  tme  spectrum.  The  tapers  with  the  least  amount  of  leakage,  the  zeroth  order 
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371  and  4k  tapers,  have  essentially  constant  weighting  of  one.  Tapers  with  increasing  leakage 
ate  weighted  progressively  less. 

Test  2'  Displacement  pulse  with  added  noise 

Here  spectra  are  computed  from  the  same  time  series  with  a  slight  amount  of  white  noise 
added,  as  shown  in  Figure  A5a.  These  spectral  estimates  differ  from  those  in  Figure  A3  only 
beyond  about  20  Hz,  where  the  signal  to  noise  level  becomes  relatively  low.  The  20%  cosine 
taper  displays  the  most  variation  here,  while  the  471  tapers,  because  of  their  large  resolution 
width,  produce  the  smoothest  estimate. 

The  eigenspectra  and  weights  for  this  case  are  shown  in  Figure  A6.  The  weights,  in  gen¬ 
eral,  are  reduced  relative  to  Test  1  due  to  the  increased  spectral  variation  caused  by  the  added 
noise.  The  relative  reduction  in  weights  becomes  greater  for  the  higher  order  tapers,  even  at 
the  low  frequencies  where  the  signal  to  noise  ratio  is  relatively  large. 

Test  3:  Displacement  pulse  with  two  added  harmonics 

To  contrast  resolution  and  leakage  characteristics  among  the  tapers,  harmonics  at  20  Hz 
and  55  Hz,  differing  in  amplitude  by  a  factor  of  ten,  are  added  to  the  time  series,  as  shown  in 
Figure  A7.  No  noise  has  been  added  to  the  time  series.  The  20%  cosine  estimate  produces 
the  best  resolution  of  the  harmonics,  as  well  as  the  greatest  amount  of  spectral  leakage.  Were 
the  relative  amplitude  difference  between  the  harmonics  larger,  the  presence  of  the  smaller 
amplitude  harmonic  could  easily  be  lost  in  the  leakage  of  the  larger  harmonic.  The  multi-taper 
estimates  are  all  quite  smooth,  displaying  very  little  spectral  leakage  effects.  The  resolution 
differences  between  the  multitaper  estimates  are  clearly  evident  at  the  harmonic  frequencies. 
All  estimates  retain  the  high-frequency  slope  about  equally  well,  however,  there  is  a  progres¬ 
sive  blurring  of  the  comer  frequency  proceeding  from  the  27i  to  the  47t  estimate  due  to 
increased  smoothing  of  the  20  Hz  harmonic.  This  pffeef  v’ould  not  be  as  prc..ounccd  were  rhe 
time  series  of  longer  duration,  which  would  narrow  the  resolved  width  of  the  harmonics. 
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The  eigenspectra  and  weights  for  this  case  are  shown  in  Figure  A8.  Rather  than  mono- 
tonically  decay  as  in  the  first  lest,  here  the  weights  increase  around  the  frequency  of  the  har¬ 
monics.  Note  that  the  weighting  of  the  higher  frequency  harmonic  is  scaled  by  the  background 
weighting,  (i.e.  the  weighting  shown  in  Figure  A4b)  and  so  for  any  taper  it  is  weighted  rela¬ 
tively  less  than  the  lower  frequency  harmonic. 

Comments 

These  synthetic  tests  have  displayed  some  of  the  features  of  the  multi-taper  spectral  esti¬ 
mation  method.  For  each  test  case  we  saw  how  reduced  spectral  leakage  from  higher  order 
tapers  trades  off  with  increased  variance  through  the  adaptive  weighting  of  eigenspectra.  We 
also  examined  variance  and  resolution  characteristics.  In  Tests  1  and  2,  the  multi-taper  method 
reproduced  the  spectrum  for  a  simple  explosive  source  better  than  that  from  simple  cosine 
tapering,  particularly  at  low  frequencies.  The  accuracy  of  the  estimate  improved  proceeding 
from  the  27t  to  the  47t  estimate,  as  more  and  more  tapers  were  added,  which  reduced  the  vari¬ 
ance.  However  in  Test  3,  where  small  bandwidth  structure  was  added  to  the  true  spectrum,  the 
low-variance  4k  estimate  was  no  longer  the  obvious  choice  as  the  best  estimate.  The  higher- 
resolution  2k  tapers  appeared  to  provide  a  better  overall  estimate  for  this  case.  And  it  was  the 
20%  cosine  estimate  that  best  reproduced  the  harmonics  by  virtue  of  its  narrow  resolution 
width. 

In  general,  the  choice  of  spectral  estimation  method  should  depend  on  what  is  known 
about  the  data,  what  one  is  looking  for  in  the  spectrum,  and  the  statistical  requiremenLs  of  the 
estimate.  For  example,  if  one  is  looking  for  harmonics,  then  high-resolution  methods,  such  as 
the  split  cosine  and  maximum  entropy  estimates,  would  be  appropriate.  However,  if  details  in 
the  spectrum  on  scales  smaller  than  the  relatively  larger  resolution  widths  of  the  multi-tapers  is 
not  terribly  significant  and  not  of  great  interest,  then  the  multi-taper  estimate  may  be  favored. 
Such  could  be  the  case,  for  example,  when  analyzing  time  scries  from  explosive  sources,  as 
suggested  by  the  above  test  ca.scs.  In  any  case,  it  is  instructive  to  experiment  with  all  tapers 
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for  a  particular  data  set  to  get  a  feel  for  the  differences  and  for  which  estimate  may  be 
appropriate.  This  is  done  below  for  data  recorded  at  the  Ruby  Valley  array. 

Ruby  Valley  Example 

In  this  section  we  compare  the  differences  in  the  spectral  estimates  for  seismic  data 
recorded  at  th'*  Ruby  Valley  array  and  the  LLNL  seismic  station.  The  time  series  here  consists 
of  a  2.6-second  window  of  the  radial  Pn  wave  from  the  explosion  HARDIN,  the  same  window 
analyzed  in  section  3.3.1.  Figure  A9  shows  the  instrument  corrected  spectral  estimates  for 
each  type  of  taper.  The  Ruby  Valley  spectrum  shown  is  actually  the  mean  of  the  spectra  com¬ 
puted  for  each  of  the  twelve  stations  at  the  array,  assuming  a  log-normal  spectral  distribution. 
Also  shown  is  the  recording  system  noise  level  at  the  Ruby  Valley  array. 

The  differences  here  are  similar  to  the  differences  found  in  the  synthetic  cases,  namely 
the  decreasing  spectral  leakage  and  decreasing  resolution  proceeding  from  the  20%  taper  to  the 
4tc  tapers.  All  estimates  produce  essentially  the  same  high-frequency  slope  of  approximately  -7 
beyond  10  Hz.  The  spectra,  even  though  from  an  explosive  source,  do  show  prominent  peak¬ 
ing,  possibly  due  to  site  response  effects.  These  features,  near  6  Hz  at  the  array  and  near  3  Hz 
and  8  Hz  at  the  LLNL  site,  remain  visible  even  up  to  the  low- resolution  4k  estimate.  Here 
one  might  choose  to  favor  the  37t  estimate  as  a  compromise  between  the  2k  and  47t  estimates, 
retaining  relatively  low  variance,  characteristic  of  the  47t  estimate,  and  acceptable  resolution 
characteristic  of  the  37t  estimate. 

For  compari.son  to  the  test  cases,  the  eigenspectra  and  multi-tajxr  weights  for  station  #1 
of  the  Ruby  Valley  array  are  shown  in  Figure  AlO.  The  eigen.spectra  shown  here  have  not 
been  instrument  corrected.  Notice  the  significant  drop  in  weighting  for  the  lower  order  tapers 
near  30  Hz.  Referring  to  Figure  A9,  this  is  just  where  the  signal  runs  into  the  noise  and  the 
spectral  variation  greatly  increases,  resulting  in  the  low  weighting. 
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2n 


3n 


FRACTIONAL  TIME 


471 


split  cosine  taper  and  the  first  L  =  2P  -  1  multiple  tapers,  for  P  =  2,  3.  and  4. 
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TAPER  TRANSFORMS 


FREQUENCY  BIN  FREQUENCY  BIN 


FKJ.  A2  Discrete  Fourier  Transform  amplitude  spectra  of  the  tapers  shown  in  f-i^ture  Al,  plo;- 
ted  as  a  function  of  frequency  bin  number  J,  where  J/T  are  the  sampled  frequencies  relumed 
from  a  fast  fourier  transform.  T  being  the  length  of  the  time  series  in  seconds 
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Multi-taper  Eigenvalues 

order 

27C 

3n 

471 

0 

0.99994 

0.99999987 

0.9999999997 

1 

0.998 

0.999991 

0.99999997 

2 

0.96 

0.9997 

0.9999988 

3 

0.73 

0.995 

0.99997 

4 

0.95 

0.9994 

5 

0.72 

0.993 

6 

0.94 

7 

_ 

0.72 

Table  A1  Eigenvalues  X  for  ihe  lower  order  lapers  of  uie  2k,  Sir,  and  4n  fami¬ 
lies  of  tapers.  The  amount  of  fractional  specu-al  leakage  for  each  taper  is  I  - 


ly? 


FIG.  A3  (a)  Discrete  time  series  of  the  thcorcticai  displacement  from  a  spherically  symmetric 
explosive  source,  (b)  Corresponding  theoretical  displacement  amplitude  spectrum  and  esti¬ 
mates  based  on  the  tapers  indicated,  (c)  Log-linear  plot  of  (b). 
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FU;.  A4  Eigenspcctra  for  the  27i,  37i,  and  4rt  families  of  multiple  tapers  computed  from  the 
time  series  in  Figure  A3.  The  corresponding  frequency  dependent  eigenspcctra  weighting 
functions  lw;(OI^  are  shown  to  the  right 
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FIG.  A5  Same  as  Figure  A3  with  white  noise  added  to  the  time  series. 
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FKf.  A6  Same  as  Figure  A4  with  white  noise  added  to  the  lime  series. 


FIG.  A7  Same  as  Figure  A3  with  harmonics  added  to  the  time  series  at  20  Hz  and  55  Hz.  No 
noise  added. 
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FIG.  A8  Same  as  Figure  A4  with  harmonics  added  to  ihe  lime  series 
noise  added. 
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HARDIN  Pn 


FIG.  A9  Instrument  corrected  displacement  amplitude  spectral  estimates  of  the  Pn  wave  from 
the  explosion  HARDIN  recorded  at  the  Ruby  Valley  array  and  at  the  LLNL  site  for  the  tapers 
indicated.  The  Ruby  Valley  spectrum  shown  represents  an  average  of  the  spectra  computed  at 
each  of  the  12  array  recording  stations.  The  Ruby  Valley  system  noise  is  also  shown. 
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FIG.  A 10  Eigenspiectra  and  corresponding  weights  at  station  #  i  of  the  Ruby  Vailcy  array  for 
the  Pn  wave  of  the  cx^-losion  HARDIN. 
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Appendix  B 


Bandpass  Filtering  Example  of  a  Regional  Waveform 


An  important  task  of  regional  seismic  monitoring  is  the  detection  and  isolation  of  regional 
phases.  This  is  often  done  in  the  time  domain  with  the  use  of  STA^LTA  algorithms,  which 
may  be  made  more  effective  by  first  bandpass  filtering  the  data.  For  example,  regional  data 
recorded  at  the  NORESS  array  often  reveal  otherwise  concealed  Pg  and  Sn  phases  when 
bandpassed  filtered  toward  the  higher  frequency  end  (Mykkeltveit  et  al. ,  1990) 

As  an  example  of  our  western  U.S.  data,  bandpass  filtered  waveforms  from  one  of  the 
Ruby  Valley  array  recordings  of  the  explosion  HARDIN  are  displayed  in  Figure  B 1 .  The  raw 
three-component  data  are  shown  at  the  top  of  the  figure.  A  12  pole  zero-phase  Butterworth 
filter  was  used  to  compute  the  filtered  traces.  The  comer  frequencies,  or  bandpass,  of  the 
filters  increase  in  octave  steps,  as  indicated  in  the  figure.  Each  waveform  is  plotted  to  74 
seconds  and  is  normalized  by  its  maximum  amplitude.  On  the  raw  traces  can  be  seen  the 
onsets  of  Pn  near  2  seconds,  Pg  near  10  seconds,  and,  on  the  horizontals,  Lg  near  60  seconds. 
In  the  0.5-4  Hz  range  the  waveforms  are  dominated  by  Pg  and  Lg.  In  the  4-8  Hz  range  Pn  has 
surpassed  Pg  in  amplitude  on  the  vertical  component  and  Lg  is  greatly  attenuated.  By  the  8-16 
Hz  range,  Pn  has  reached  equal  or  greater  amplitude  relative  to  Pg  on  all  components  and  Lg 
appears  undetectable.  Here  the  peak  amplitude  of  Pn  is  about  the  same  on  all  three  com¬ 
ponents.  Finally,  in  the  16-32  Hz  range  Pn  on  the  vertical  component  sustains  the  largest 
amplitudes,  while  on  the  horizontals  Pn  falls  back  below  Pg.  In  this  frequency  band  Pg  is 
detectable  only  on  the  horizontals,  and  here  the  detectability  is  not  great.  In  no  ir‘'^tance  was  a 
clear  new  arrival  revealed  in  the  filtering  process.  However  we  do  see  that  the  detectability  of 
Pn,  Pg,  and  Lg,  measured  here  simply  by  the  signal  amplitude  relative  to  the  preceding  ground 
motion,  is  dependent  on  frequency  and  component  of  ground  motion.  Overall  the  detectability 
of  Pn  appears  greatest  on  the  vertical  component.  The  detectability  of  Pg  appears  comparable 
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on  all  components,  while  Lg  is  greatest  on  the  horizontals  at  the  lowest  frequencies. 
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Appendix  C 


Coherence  Estimation  and  Statistics 

The  purpose  of  this  appendix  is  to  describe  the  procedure  used  in  the  estimation  of  coher¬ 
ence  and  to  describe  the  corresponding  statistics  of  the  estimate.  Both  the  multiple  taper  spec¬ 
tral  estimation  procedure  of  Thompson  (1982),  and  the  more  conventional  frequency  averaging 
method  (e.g.  Beauchamp  and  Yuen,  1979)  are  considered. 

First,  recall  that  the  estimate  of  coherence  between  two  time  series  Uj(t)  and  Ui^(t)  at  fre¬ 
quency  f  is  given  by 

ISjkCOl 

- - - nr  •  (C.i) 

[S^(f)  Sy,(f)]>^ 

where  Sji^  is  the  cross-spectral  estimate  given  by 

Sjk(f)  =  £[Uj(Ouk(f)*]  ,  (C.2) 

where  E  is  the  expectation  operator,  Uj(0  and  U|t(f)  are  the  Fourier  U'ansforms  of  the  time 
series,  and  *  denotes  the  complex  conjugate.  The  value  of  depends  fundamentally  on  the 
means  by  which  the  expectation  in  Sji^CO  is  computed.  The  non-stationarity  typical  of  our 
seismic  data  precludes  averaging  cross-spectra  over  multiple  time  windows,  as  can  be  done 
with  stationary  time  series  (e.g.  Beauchamp  and  Yuen,  1979).  Therefore  some  other  means  of 
spectral  averaging  must  be  used.  Two  alternatives  are  described  below. 

Frequency  Averaging  Method 

In  this  method  the  assumption  is  made  that  the  cross-spectrum  varies  slowly  with  fre¬ 
quency,  so  that  we  can  average  the  cross-spcctrum  over  a  frequency  band  centered  about  the 
frequency  of  interest.  A  single  data  taper,  such  as  a  split  cosine  taper,  is  used  on  each  of  the 
two  time  series.  An  estimate  of  the  cross-spcctrum  between  the  time  series  at  a  “center 
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frequency”  ^  is  obtained  by  a  weighted  average  over  a  frequency  band  of  width  2M/T  cen¬ 
tered  at  fc,  i.e., 

Sjk(fc)=  E  a(fj  Uj(f  J  u,(fj*,  (C.3) 

m=-M 

where  2M+1  discrete  Rayleigh  frequencies  are  used  in  the  averaging  (f^  =  fc  +  ni/T)  and  aff^) 
is  a  weighting  function.  Typical  weighting  functions  are  the  boxcar,  triangle,  and  Hamming 
functions.  The  weights  are  normalized  such  that  J^a^  =  1  to  prevent  biasing  the  estimate. 
Again,  this  estimate  of  SjjtffJ  assumes  that  the  cross-spectral  values  at  each  frequency  within 
the  averaging  bandwidth  are  independent  and  identically  normally  distributed. 

Multiple  Taper  Method 

An  alternative  to  the  above  frequency  averaging  is  the  method  of  Thomson  (1982),  which 
employs  the  use  of  multiple  data  tapers.  As  described  in  Appendix  A,  these  tapers 
v/(tn;T,W);  /=!,.. .,L  are  designed  to  minimize  spectral  leakage  outside  a  chosen  frequency 
band  and  are  used  to  construct  eigenspectra  Uj/(0.  If  the  local  variations  in  the  spectrum  are 
not  too  extreme,  the  eigenspectra  will  be  nearly  uncorrelated  and  the  estimate  of  the  cross¬ 
spectrum  can  be  constructed  from  an  ensemble  average  at  a  single  frequency,  rather  than  over 
a  bandwidth  as  before: 

L 

Sjk(0  =  Z  a,  Uj;(0  UkKO*.  (C.4) 

/=i 

The  coherence  is  then  given  as  before  in  (C.l).  Here  the  weights  a/  are  a  function  of 
taper  order  number  rather  than  frequency  in  (C.3).  The  cross-spectrum  estimated  in  this  way 
was  used  by  Thomson  (1982)  in  synthetic  coherency  calculations,  and  later  by  Park  etal., 
(1987)  in  a  polarization  analysis  of  seismic  data.  Because  the  coherence  is  essentially  a  corre¬ 
lation  coefficient,  if  only  one  realization  of  Uj(f)Uk(0*  is  used  in  the  averaging  for  Sj^(0  the 
correlation  will  be  perfect  and  will  always  equal  1.  Therefore,  to  obtain  meaningful 
results  we  must  have  M  significantly  greater  than  0  in  (C.3)  and  L  significantly  greater  than  1 
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in  (C.4). 


Notice  that  in  (C.4)  frequency  dependent  weights  ha’/e  not  been  applied  to  the  eigenspec- 
tra  as  they  were  in  the  multi-taper  estimate  of  the  amplitude  spectrum  described  in  Appendix 
A.  This  is  because  the  low  weighting  which  would  be  applied  to  the  higher  order  tapers  is 
effectively  equivalent  to  decreasing  L  in  (C.4).  This  would  result  in  a  high-variance  estimate 
of  biased  towards  a  value  of  one.  However,  without  weighting,  spectral  leakage  of 
higher  order  tapers  can  be  significant.  We  can  reduce  this  effect  by  using  fewer  higher-order 
tapers  in  (C.4).  If  instead  of  considering  the  first  3,  5,  and  7  of  the  2k  ,  Stc,  and  47c  tapers, 
respectively,  as  we  did  in  the  spectral  amplitude  estimate,  and  iastead  consider  the  first  1,  3, 
and  5  tapers,  the  maximum  amount  of  leakage  from  any  taper  will  be  reduced  to  below  1% 
(see  Table  1  in  Appendix  A).  In  doing  so,  we  have  reduced  the  bias  from  spectral  leakage,  but 
now  wc  may  be  in  danger  of  unacceptably  large  variance  due  to  the  smaller  number  of  tapers 
u.sed.  To  further  quantify  this,  and  to  compare  the  frequency  averaging  and  multi-taper  coher¬ 
ence  estimation  methods,  we  now  look  at  the  st?tistics  of  the  estimates. 


Coherence  vStatistics 

The  suitistics  of  coherence  arc  quite  complicated.  It  is  not  normally  distributed  and  does 
not  have  a  constant  variance,  i.e.,  the  variance  of  lYj|(l  depends  on  the  value  of  lyj^l.  This 
makes  lea.st  squares  fitting  of  coherence  data  non-trivial.  However,  under  a  tanh"’  variance 
stabilizing  transformation,  lyji(l  can  be  transformed  to  a  variable  which  is  normally  distributed 
and  has  a  constant  variance  (Enochson  and  Goodman,  1965),  thus  making  it  more  suitable  to 
least  .squares  procedures  (e.g.  Abraham.son,  1988). 

Let  lyj|(l  denote  the  true  coherence.  If  the  value  of  lyji^l  is  not  too  small  and,  again, 
a.s.suming  that  the  Fourier  coefficients  in  (C.3)  or  (C.4)  arc  independent  and  normally  distri¬ 
buted,  then  the  distribution  of  tanh^My^^l  is  approximately  normal  with  mean  and  variance 
given  by 


=  tanh'lyj  +  — £1^ 
'  2(1 -g‘) 


:i  1 


(C.5) 


and 


var  tanh  ’lYjkij  = 


(C.6) 


(BriJlinger,1981;  Bloom  field,  1976).  The  quantity  is  given  by  Bloomfield  n976)  to  approxi¬ 
mate  the  statistical  effect  of  the  data  taper  v(tn)  and  frequency  averaging,  where 

N 


n=l 

- - - 

N 


M 


Tn=^M 


(C.7) 


Analogously,  with  multi-taper  ensemble  averaging,  one  can  show  that 

\4 


g  =  N  X 
/=! 


n=l 


Zv/(tn)' 

n=\ 


2 


(C.8) 


In  general,  a  taper  which  downweights  the  time  series  more  will  have  a  larger  g^  factor 
and  therefore  produce  a  larger  bias  g^/2(I-g^)  and  variance  g^/2.  However,  at  the  same  time, 
such  a  taper  typically  results  in  relatively  less  spectral  leakage.  Such  is  the  case,  for  example, 
in  contrasting  the  20%  cosine  taper  with  the  lower  order  multiple  tapers  (refer  to  Figures  A1 
and  A2  in  Appendix  A).  Therefore,  we  should  expect  the  multi-tapers  to  produce  somewhat 
greater  bias  and  variance  as  defined  in  (C.5)  and  (C.6).  Keep  in  mind,  however,  the  significant 
reduction  in  bias  due  to  spectral  leakage  for  the  multi-lapcrs,  not  accounted  for  in  (C.5)  and 
(C.6). 

In  regard  to  frequency  averaging,  by  increasing  the  value  of  M  in  (C.7),  i.e.,  the  number 
of  frequencies  over  which  Sj^  is  averaged,  one  gains  a  decrease  in  the  above  bias  (C.5)  and 
variance  (C.6).  However,  in  doing  so,  resolution  is  lost  as  the  averaging  bandwidth  increases 
with  increasing  M.  Additionally,  incrca.sing  M  increases  the  bias  which  will  result  when  the 
spccira  are  not  identically  distributed  over  the  averaging  bandwidth.  Concentrating  the  fre¬ 
quency  averaging  weights  towards  the  center  frequency,  such  as  with  a  Hamming  or  triangle 
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function,  will  reduce  this  bias  but  at  the  same  time  will  tend  to  increase  the  bias  in  (C.5). 

The  multi-taper  estimate  in  (C,4),  in  addition  to  its  resistance  to  spectral  leakage,  has  the 
added  advantage  of  not  losing  resolution  or  gaining  additional  bias  from  frequency  averaging. 
However,  as  mentioned  above,  the  ensemble  average  over  only  L  tapers  may  produce  unac¬ 
ceptable  variance,  where  we  have  now  reduced  L  to  1,  3,  and  5  for  the  2n,  3k,  and  esti¬ 
mates,  respectively.  To  reduce  the  variance,  one  can  supplement  the  multi-taper  ensemble 
average  with  a  small  amount  of  frequency  averaging.  Having  multiple  cross-spectral  estimates 
at  each  frequency,  one  can  still  achieve  relatively  better  resolution  than  that  from  the  single¬ 
taper  frequency  averaging  method  without  suffering  greatly  increased  bias  and  variance. 

To  demonstrate  this  more  clearly,  the  bias  from  (C.6)  and  the  standard  deviation  a  =  g/V2 
from  (C.7)  are  plotted  in  Figures  Cl  and  C.2,  respectively,  as  a  functions  of  Nx,  the  number  of 
the  cro.ss-spcctra  averaged  in  the  estimate  of  lyj^l.  For  frequency  averaging,  Nx  =  2M-t-l, 
where  a^  is  a  Hamming  window  function  and  v(tn)  is  a  20%  cosine  data  taper.  For  the  multi¬ 
tapers,  a  boxcar  weighting  function  is  used,  3/  =  1/L.  The  multi-taper  curves  begin  as  a  func¬ 
tion  of  increasing  taper  order  up  to  Nx  =  L.  This  taper  sequence  is  then  continuously  repeated 
to  demonstrate  the  effect  of  including  frequency  averaging  over  additional  frequencies  to  cither 
side  of  a  center  frequency.  For  example,  Nx  =  5  on  the  371  curve  indicates  that  in  addition  to 
the  L3„  =  3  tapers  used  at  the  center  frequency  f^,  the  lowest  order  37i  tapper  was  also  used  on 
the  two  adjacent  frequencies,  f  =  fc  ±  1/T.  Similarly,  Nx  =  9  implies  all  three  37t  tapers  used 
on  the  two  adjacent  frequencies. 

Below  about  Nx  =  5,  the  bias  and  variance  are  exceedingly  large  for  each  of  the  four 
methods.  Beyond  about  Nx  =  5,  the  tapers  which  downweight  the  time  series  relatively  less, 
such  as  the  20%  cosine  and  the  27i  tapers,  do  produce  the  lowest  bias  and  variance  as 
expected,  however  the  differences  between  the  four  methods  arc  not  great.  The  principal 
differences  here  lie  in  the  resolution  of  the  estimates.  Consider  an  example  where  one  decides 
that  acceptable  values  of  the  bias  and  o  occur  near  Nx  =  15  (  all  four  estimates  give  about  the 
same  values  here).  To  achieve  Nx  =  15,  the  20%  cosine  single-taper  estimate  would  require 
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averaging  over  the  seven  adjacent  Rayleigh  frequencies  to  either  side  of  the  center  frequency. 
With  its  2/T  spectral  resolution  width,  where  T  is  the  record  length,  this  results  in  a  coherence 
estimate  averaged,  primarily,  over  a  total  of  16/r  Hz.  The  27t  estimate,  with  =  1  and  hav¬ 
ing  a  greater  4/T  resolution  width,  will  average  over  18/T  Hz.  However,  the  37t  and  47c  esti¬ 
mates,  with  L3„  =  3  and  L4„  =  5,  and  having  resolution  widths  of  6/r  and  8^,  respectively, 
will  both  average  over  a  total  of  only  lO/T  Hz.  Thus  we  see  that,  in  addition  to  their  resis¬ 
tance  to  spectral  leakage,  the  3re  and  4iz  coherence  estimates  will  produce  essentially  the  same 
statistics,  as  defined  by  (C.5)  and  (C.6),  as  the  20%  cosine  single-taper  estimate  but  with 
significant  improvement  in  spectral  resolution,  lO/T  Hz  as  opposed  to  16/T  Hz  for  =  15. 
Also,  with  this  better  resolution,  the  estimates  are  not  as  susceptible  to  bias  due  statistical  vari¬ 
ations  in  cross-spectra  over  wide  averaging  bandwidths. 

The  other  important  statistic  is  the  null  distribution  of  coherence.  This  is  the  value  of 
coherence  above  which  one  can  be  100p%  confident  that  the  estimate  was  not  produced  from 
random  noise.  The  lOOp  %  confidence  point  of  the  null  distribution  of  lyjkl^  can  be  approxi¬ 
mated  by 

p\p)  =  1  -  (C.9) 

(Bloomfield,  1976).  The  50%  and  90%  confidence  levels  for  both  tanh^Myl  and  lyl  are  shown 
in  Figure  C3  as  functions  of  Nx  for  the  taper  and  frequency-averaging  parameters  described 
above.  These  curves,  like  those  in  Figures  Cl  and  C2,  can  be  used  to  determine  how  much 
frequency  averaging  is  required  for  acceptable  statistics  for  the  four  different  estimation  pro¬ 
cedures.  Here  again,  the  20%  cosine  frequency-averaged  estimate  produces  somewhat  more 
optimum  statistics,  in  this  case  lower  values  for  the  null  distribution.  However,  as  before,  the 
differences  here  are  not  terribly  great  and,  as  discussed  above,  the  multi-taper  method  offers 
significant  advantages  in  resolution  and  decreased  bias  from  spectral  leakage,  making  it  a  more 
attractive  estimation  procedure  for  coherence,  and  cross-spectra  in  general,  for  non-stationary 
time  series. 


214 


References  Cited 


Abrahamson,  N.  (1988),  Empirical  models  of  spatial  coherency  of  strong  ground  motion,  per¬ 
sonal  manuscript. 

Beauchamp,  K.  and  C.  Yuen  (1979),  Digital  Methods  for  Signal  Analysis,  George  Allen  & 
Unwin,  London. 

Bloomfield,  P.  (1976),  Fourier  Analysis  of  Time  Series:  An  Introduction,  John  Wiley  &  Sons, 
New  York. 

Brillinger,  D.  R.  (1981),  Time  Series,  Data  Analysis  and  Theory,  Expanded  edition,  Holden- 
Day,  San  Francisco. 

Enochson,  L.  D.  and  N.  R.  Goodman  (1965),  Gaussian  approximations  to  the  distribution  of 
sample  coherence.  Tech.  Rep.  AFFDL-TR-65-57,  Wright-Patterson  Air  Force  Base. 

Park,  J.,  F,  L.  Vernon  III,  and  C.R.  Linberg  (1987),  Frequency  dependent  polarization  analysis 
of  high-frequency  seismograms,  7.  Geophys.  Res.,  92,  12664-12674. 

Thomson,  D.J.  (1982),  Spectrum  estimation  and  harmonic  analysis,  IEEE  Proc.,  70,  1055- 
1096. 


215 


Nx 


FIG.  Cl  Bias  of  lanh'Myl  as  a  function  of  N^.  the  number  of  cross-spcctra  used  in  the  esti 
mate  of  ly',  for  each  of  the  four  estimates  described  in  the  text.  This  does  not  include  the  bia,' 
which  will  result  from  frequency  averaging  when  the  cross-spectra  are  not  identically  statist! 
cally  distributed,  nor  the  bias  due  to  spectral  leakage. 


2 


tanh  ^  lyl  standard  deviation 


FIG.  C2  Standard  deviation  of  tanh  ‘  lyl  as  a  function  of  the  number  of  cross-spectra  used  in 
each  of  the  four  estimates. 
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of  cross-spectra  used  in  each  of  the  four  cslimaies.  The  curve  symbols  have  the  same  meaning 
as  in  Figures  Cl  and  C2.  (b)  Same  as  (a)  but  for  coherence  I7I  . 


Appendix  D 


Estimates  of  tanh  ^Yl  Used  in  Chapter  5  Inversions 


This  appendix  contains  the  estimates  of  tanh'Myl  used  in  the  inversions  for  coherence 
model  parameters  in  Chapter  5.  Estimates  for  Pn  and  Pg  at  each  array  site  up  to  14.9  Hz  are 
displayed.  The  best-fitting  isotropic  models  are  superimposed  on  the  estimates. 
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